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ABSTRACT 


The present report is mainly concerned with a careful analysis of the fundamental 
ideas, methods, and postulates which underlie current approaches to the cosmological 
problem. In particular, the methods of general relativity and the author’s kinematical 
relativity are contrasted. In the later method all the fundamental notions (‘“‘uniform 
time,”’ “‘rigid length scale’’) are carefully scrutinized and emphasis is laid on the methods 
by which two arbitrary observers can set up ‘“‘congruent”’ clocks and distance measures 
by dispatching and receiving light-signals. A linear equivalence, defined as a class of 
observers possessing congruent clocks, is then analyzed. This leads to the considera- 
tion of the two scales of time, ¢ and 7 = ¢, log t/t) + to, such that a uniform motion 
equivalence in the scale of ¢ becomes one at relative rest in the scale of r. The physical 
significance of the two scales of time are then further discussed. The methods by 
which the dynamical equations of motion are to be set up by projecting a “free” 
test particle in the presence of a density distribution and specified equivalence are then 
considered. The fundamental problem here is the determination of the motion of a free 
particle for an equivalence satisfying the cosmological principle, i.e., the case when the 
equivalence is such that the statistical description of the whole equivalence from any 
one member of the equivalence coincides with that from any other member. The solu- 
tion of the cosmological problem which the consideration of this problem leads to is 
then discussed; particular attention is here given to the degree of arbitrariness which is 
left at each stage of the solution of the problem and how the arbitrariness is then later 
removed. The kinematic arguments are shown to lead to (1) the law of inertia for a 
substratum (i.e., a density distribution satisfying the cosmological principle), (2) the 
existence of ‘‘gravitational’’ mass, and (3) the inverse-square character of ‘‘gravita- 
tional force.’’ Further problems which can be treated by the methods outlined in this 
paper are also indicated. The importance of the two time scales, that of ¢ and that of 7, 
in answering questions about the universe (e.g., whether the universe is expanding or 
is nonexpanding) is pointed out. 


Cosmological theory is concerned with the distribution of matter 
and motion in the universe as a whole. Any particular cosmological 
« This lecture was delivered at the Symposium on Galactic and Extragalactic Struc- 
ture, held in connection with the dedication of the McDonald Observatory on May 5-8, 


1939. 
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theory attempts to indicate the relationship between the distribution 
of matter (the density distribution) and the distribution of motion 
(the velocity distribution relative to a specified observer). Or, to 
put it in another way, a cosmological theory may attempt to predict 
the main outlines of motion in the universe from a knowledge of the 
main outlines of the density distribution; or, more ambitiously, to 
show why a certain main outline of motion-distribution should even- 
tually lead to some particular—say, the observed—distribution of 
matter. 

There are, fundamentally, two different possible starting points 
for an investigation of this kind-two possible modes of attack. One 
is to suppose the investigator to be equipped with all the knowledge, 
empirical (that is, experimental and observational) and theoretical 
(that is, deductive) that has been acquired in the course of the age- 
long history of science up to the present day; to be in possession of 
all the relevant laws of nature; and to attempt to combine these 
various items of knowledge just as in attacking any other scientific 
problem, by making initial hypotheses, by applying to these the ac- 
cepted laws of nature, and by determining whether the deductions 
agree with observation. If they do so, the initial hypotheses may be 
presumed to be in some measure confirmed; or, at least, the way is 
open to an investigation as to whether any other initial hypothesis 
would give an equally successful agreement with observations. The 
relevant laws of nature for such an investigation will include the 
laws of motion (of dynamics) and of gravitation, bound together by 
the principle of relativity. I use the latter phrase out of deference to 
the reader’s possible predilections, although I confess that I do not 
know what the principle of relativity really is, as I do not, in my own 
work, make any appeal to, or application of, it. Still, though now- 
adays the principle of relativity tends to be lost in wider and more 
general modes of thought, ‘‘relativity,”’ the word itself, stands for a 
technique, a calculus; and the calculus will most probably be utilized 
by the investigator who looks at it from the point of view which I 
have described. This calculus, or technique, includes the attribution 
of intrinsic properties to ‘‘mere space,” the combination of the con- 
cept of space with the concept of time to yield the concept “‘space- 
time,’ and the linking of the intrinsic properties of this con- 
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cept ‘‘space-time’’ with the occupation of space by matter and with 
the motions of this matter that ensue within it. Accepting, then, all 
these doctrines and points of view, the investigator will proceed to 
attack the problem of the universe as a whole as if it were any other 
dynamical problem. He simply broadens the field of investigation 
until it includes the totality of existent matter. He synthesizes his 
solution out of locally acquired knowledge of nature; and he is en- 
titled, in this mode of solution, to make use of all his sources of 
knowledge, however obtained. 

It will be observed that in such a solution three types of knowledge 
are involved: intuitive, deductive, and empirical. Among items 
of intuitive knowledge is the awareness of the passage of time for the 
individual observer, about which I shall have more to say later. 
Among items of deductive knowledge are the consequences of the 
notions of number, of point, line, and plane, and, in short, the the- 
orems of arithmetic, algebra, and geometry. Among items of empiri- 
cal knowledge is the empirical fact that the number of spatial di- 
mensions or degrees of freedom of a point is three, together with the 
empirical laws of dynamics and gravitation, as already mentioned. 
The empirical laws of dynamics and gravitation have, of course, a 
different status from such a brute fact as the number of spatial di- 
mensions. Their ‘‘sanction,”’ 1.e., the authority behind them, is their 
verification in a large number of concrete instances, or, rather, the 
verification of their consequences when they have been combined 
with deductive elements. That is to say, they are inductive generali- 
zations based on the supposed ‘‘principle of induction,” which is so 
difficult to state and still more difficult to establish. 

There is, of course, nothing illegitimate in this, the first, method of 
attacking the cosmological problem. It is the method adopted for 
all smaller problems. Why, then, should we not be content with it 
for the grandest of all problems? The answer is that, even if it can 
be applied completely, it leaves our curiosity largely unsatisfied, for 
why should an observer situated on a body, a member of the universe, 
have to become aware of the results of a large number of small-scale 
experiments before he can understand the spectacle presented to 
him? The spectacle itself constitutes a grand experiment in being. 
Should not this experiment tell its own tale and point its own moral? 
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Moreover, to say that a large-scale system can be understood be- 
cause it exemplifies principles known to be true in small-scale sys- 
tems gives us no insight into why these principles hold good at all. 
To take the principles as given is to admit an element of magic, of 
the irrational, in the working of the world, for if they have no fur- 
ther exaplanation, they are just magical. It is the claim of the gen- 
eral theory of relativity that it rests on observation and experiment 

for example, on the fact that inertial mass and gravitational mass 
are equal. In so doing, it necessarily admits, great as are its achieve- 
ments, that it is not fundamental, for it is always possible to ask 
deeper questions, such as ‘‘Why is gravitational mass equal to in- 
ertial mass?”’ and ‘‘Why should the occupation of space by matter 
cause what relativists call ‘curvature of the space’?”” We are en- 
titled, it appears to me, to ask these and many other deeper ques- 
tions, such as ‘‘Why does such and such a law of nature, supposedly 
irreducible to any simple law, hold good?” ‘‘What is meant by the 
‘rigid body,’ or the ‘uniform time-keeper,’ of relativistic metrology?” 
‘‘Why should there be an invariant measure of the separation of two 
events?”’ And so on. 

It was a view held by Mach, and since assented to by many think- 
ers, that some, at least, of the laws of nature, such as the law of in- 
ertia, are consequences of the distribution of matter in the universe; 
that is to say, it is not that the motions of particles in the universe 
ensue in virtue of (a) the distribution of matter in the universe, to- 
gether with () certain laws of nature; but that (0), the laws of na- 
ture themselves, hold good simply in virtue of (a), of the particles be- 
ing in the presence of the rest of the universe. If there is any validity 
in such a view, it can be ascertained or tested only by actually apply- 
ing this view, that is to say, by refraining from assuming at the out- 
set those laws of nature whose origin is in question. This leads us to 
the second main approach to a solution of the cosmological problem, 
which is to study the universe as it actually presents itself to us, with- 
out appealing to any small-scale experiments or to knowledge de- 
rived from small-scale experiments. That is, in short, to use only 
such brute facts, such irreducible facts, as are of the intuitive sort or 
do not rest on the questionable principle of induction, and thus to 
appeal to no empirical “‘laws of nature” of a quantitative kind. The 
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difference between the two methods can be appreciated by a com- 
parison with the theory of numbers or with geometry. A theorem in 
the theory of numbers may be discovered first empirically, and the 
weight of evidence for it may be such as to compel a belief in its gen- 
eral validity; but no mathematician accepts this as the last word on 
the subject. He essays to prove it, that is, to establish it deductively 
from axioms which, in effect, are definitions of the subjects of dis- 
course, without appealing to any empirical evidence. It is not that 
he despises or rejects the empirical evidence; the empirical evidence 
is most welcome as a verification of the absence of logical flaws from 
the details of the process of inference, besides its initial role in dis- 
covering the theorem. But the theorem cannot be held to be con- 
nected with the discipline which is the theory of numbers so long as 
the evidence for it is wholly or partly empirical. The mathematician 
must, at some stage, perform an act of renunciation and proceed 
without reference to empirical facts. Some empirical facts, like Fer- 
mat’s ‘‘last theorem” and Riemann’s hypothesis, have so far defied 
proof; but this does not mean that we are to regard them as inaccessi- 
ble to proof. 

Similarly, no one now tests a geometrical theorem empirically— 
no one, for example, publishes memoirs of experiments designed to 
ascertain how exactly the three perpendiculars to the sides of a tri- 
angle from the vertices meet in a point. Yet, presumably in its pre- 
Greek stage, geometry passed through such a phase in which the- 
orems were discovered or known empirically. It remained for the 
Greeks to invent new and more powerful combinations of processes 
of inference which were capable of establishing the empirically 
known theorems deductively, from axioms which (in their modern 
usage) are concealed definitions. 

There are, then, these two separate ways of attacking the struc- 
ture of the universe: the one to make use of every available piece 
of empirical knowledge known to be valid on the small scale; the 
other to begin with the situation actually presented to us by the 
totality of things without supposing ourselves to know anything, to 
start with, of the facts used in the first method. In the second meth- 
od we attempt a complete reconstruction of physics from the bottom 
upward, on an axiomatic basis. It might be supposed that the sec- 
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ond method, starting with a poorer armory of resources, could never 
get as far as the first method. Actually, it proves to be more power- 
ful, for, as it proceeds, it frees itself from the presuppositions tacitly 
made in the first method. 

I will now say something about the results of the two methods. 
In principle, the natural way of applying the first method would be 
to assume some distribution of density in the universe, apply the 
laws of dynamics and gravitation and calculate the resulting motion. 
It is natural, in the first instance, to treat each extragalactic nebula 
as a particle and to mean by “‘density at any locality”’ the aver- 
age number of such particles per unit volume; the system is then 
compared to a dust cloud, in which the separate particle-members are 
in free motions under another’s gravitational influence. The motion 
will then depend on the initially assumed density distribution. The 
method followed by general relativity is, however, indirect: it first 
assumes a metric (described picturesquely as ‘‘some kind of expand- 
ing space’), determines the motions implied, and then by means of 
dynamical equations (field equations) calculates what uniform density 
(spatially constant but a function of the time) is implied by this set 
of equations. 

It is a remarkable fact, pointed out in 1934 by McCrea and the 
writer, that the set of equations thus obtained connecting density 
and motion, obtained from the tield equations of general relativity, 
are identical with those deduced from a crude nonrelativistic New- 
tonian theory in which the main role is played by the Newtonian 
equation of motion of a particle under the gravitational influence of 
the matter inclosed within the sphere through the particle centered 
at the observer and a hydrodynamical equation of continuity. 
It then appears that there are three types of motion— motion with 
a velocity equal to, exceeding, or less than the ‘‘velocity of escape”’ 
from the sphere of matter just mentioned—-and that these three 
types of motion correspond exactly to the three possible types of ‘‘ex- 
panding space,’ namely, of zero, negative, or positive curvature. 
The interpretation of the equations in terms of observations is dif- 
ferent in the Newtonian and relativistic schemes, for they employ 
space in different senses. Nevertheless, the general dynamical con- 
tent of the two schemes is essentially the same. And the astounding 
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fact emerges that neither of these schemes, which both follow the 
first method of our general discussion, contains any criterion by 
which one can decide which to expect in nature, out of the three sup- 
posed possibilities. General relativity is powerless to tell us whether 
the space in terms of which it analyzes the universe is expected to be 
of zero, negative, or possible curvature, whether the motions to be 
expected are of parabolic, hyperbolic, or elliptic type; whether the 
universe is expected to be expanding forever, to be oscillating be- 
tween finite nonzero radii, or to be expanding and then contracting 
to zero. No theoretical criterion has been found capable of distin- 
guishing between what are, to general relativity, fundamentally dif- 
ferent possibilities. 

Science has long been acquainted with problems of the type that 
cannot be solved. Such have in the past been, for example, the prob- 
lem of perpetual motion, the problem of squaring the circle, of tri- 
secting an angle. It has usually been proved in such cases not that 
man’s intellectual resources were incompletely applied but that some 
erroneous tacit assumption has been made in the formulation of the 
problem, or, again, that the problem was a no-problem, that it did 
not exist. After men had striven for years, for example, to prove Eu- 
clid’s parallel postulate, it occurred to the founders of non-Euclidean 
geometry to do without it altogether—in short, to deny it. Is it not 
probable that this continued failure of general relativity to solve a 
fundamental problem proposed by itself means that there is some- 
thing wrong with the formulation of the problem, or with its embodi- 
ment in analysis, or with the assumptions tacitly made in stating the 
problem? I invite you to consider the most attractive alternative: 
that all these differing possibilities are different descriptions of the 
same entity-—descriptions which originate from the tacit adoption 
of different scales of time; that there is one universe but many differ- 
ent possible descriptions of it; that there is actually no diversity of 
possibility of the thing in itself, but only a diversity of possibility of 
description. This will lead us directly to the second mode of attack 
on the cosmological problem. 

One of the first questions that can be asked about a system of 
moving particles is whether it is in a steady state or not. If it isina 
steady state, its statistical description at any one instant is the same 
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as its statistical description at any other instant; the variable ¢ can- 
not enter into the description—one moment is as good as another. 
But if it is not in a steady state, the state at one instant distinguishes 
itself from that at another. It should, therefore, be possible to indi- 
cate to an intelligent person which instant is being the subject of dis- 
cussion or examination. When the same problem occurs on a less 
grand scale, when the person is interested in a system which is some 
small part of the universe, the same problem will occur; but, then, 
the person will have some experience of other small systems, into 
which experience he can fit his experience of the given small system. 
In other words, once his scale of time is fixed by some one of his in- 
dividually experienced systems, he can correlate his experience of other 
systems with that one. That one system constitutes his clock, and 
he can use it to record the evolution of other systems open to his ob- 
servation. But if his field of observation is the whole universe, and, 
following our second method, he renounces appeals to smaller 
systems, that whole universe must itself constitute his clock. And 
as one aspect of the spectacle succeeds another and so distinguishes 
itself from another, so the observer must be able to differentiate his 
description of it and refer them to different instants in his clock 
scale. To have a meaning, these must be communicable to any other 
intelligence, and, therefore, to describe them he must be able to refer 
to some natural time-zero inherent in that description of the universe 
which he is trying to communicate. The universe, if not in a steady 
state, must thus contain within its own description some identifiable 
time-zero. On the other hand, if the universe is in a steady state, any 
epoch can be taken as time-zero; past and future are illimitable; and 
the “earliest possible’? moment recedes to minus infinity on this 
time scale. 

I confess that this situation has always created for me a difficulty 
in the well-known paper which Hubble and Tolman published in the 
Astrophysical Journal in 1935.7 In that paper a certain function, 
g(t), was expanded in powers of ¢ in such a way that the epoch of 
observation was taken as ¢ = o. But this places all epochs of ob- 
servation on the same footing; it bars the development of the system 
of the universe in time, for it contains no specification of what mo- 


282, 302. 
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ment is the moment / = o. In spite of the red shift, then, the pro- 
cedure is tantamount to discussing only a stationary universe, in the 
sense that the state contemplated at any one epoch of observation is 
the same as that contemplated at any other epoch. To remove this 
restriction it is necessary to incorporate into the analysis some param- 
eter corresponding to a zero of time, i.e., to some recognizable 
event in the evolutionary history of the subject; and then to relate 
any present epoch with the introduced zero, by reckoning it from 
that zero. In the investigations contained in my book World Struc- 
ture, which were parallel to those of Hubble and Tolman, I always 
left in a time variable ¢, which could measure the epoch of observa- 
tion reckoned from some well-defined zero. It is the absence of a ¢ in 
Hubble and Tolman’s final formulae which makes me suspect that 
their analysis describes only a stationary system. Not that it is any 
the less apropos for that, as I hope to show. 

This forcing on our notice of a zero of time, this recognition that 
the universe itself must ultimately be the clock to be used in the 
second main method of approach to the cosmological problem, leads 
us to consider in more detail the problem of time as it occurs in sci- 
ence. 

In relativity the notion of uniform time is adopted without dis- 
cussion. It is said that a material clock is constituted by any “regu- 
larly periodic phenomenon,” though how one is to recognize when 
a phenomenon is strictly periodic and when not is not stated. Sup- 
pose I have two clocks, side by side. One contains a mechanism 
which makes it gain on the other by an amount which, starting from 
the epoch when the clocks agreed, varies as the square of the epoch 
read by the other. Which measures uniform time? No one can say. No 
one can catch the fleeting minute and set it side by side with another 
supposed minute. Again, we can, if we like, say that uniform time 
is that which makes the laws of dynamics hold good. But this is 
equivalent to graduating one clock—say that given by the rotating 
earth—in a particular way and calling that way “uniform.” Ob- 
viously, it is no more entitled to the word “uniform” than many 
other possible ways of graduating the clock. We have to recognize 
that the word “‘uniform” in the context of “uniform time”’ is one of 
those emotion-soliciting words which have no place in science, and 
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we must attempt to proceed without it. All modes of graduating a 
“clock”’ are then on the same footing, to begin with. 

The word “uniform” in the phrase ‘uniform time” has the same 
dangerous plausibility as the word “rigid” in the phrase “rigid body.”’ 
Most people feel that they have an intuitive sense of what is meant 
by a “rigid length scale,” just as they believe themselves to have an 
intuitive sense of what is meant by “uniform time.’ But no one can 
say offhand what is meant by either phrase. I will attempt to illus- 
trate by an example. It is a commonplace that motion is relative. It 
is perhaps less of a commonplace that position in the universe is 
relative. It is often asked whether the distant nebulae are ‘‘really”’ 
receding from us or not. I shall show later that this question has no 
meaning. For the moment, let us see what the propounder of the 
question usually means. He means to inquire: If we could put up a 
fixed “‘post,”’ say a “bench mark” in space, at a fixed distance from 
ourselves, in the vicinity of a nebular nucleus, would the nucleus be 
passing the “bench mark” or not? To make sure that the bench 
mark is at a fixed distance from ourselves, we must be continually 
measuring its distance, ideally with “rigid’’ meter scales. Suppose 
these are put end to end, so as to join ourselves rigidly to the bench 
mark. Now suppose a uniform extensible rope joining ourselves to 
the nebular nucleus, with the marks on it coinciding, at a certain 
instant ¢, with the meter markings of the rigid scale. Then the prob- 
lem of whether the nebular nucleus is “really”? moving or not is the 
same as the problem of whether the first meter marking on the rope 
is or is not passing the first meter marking on the “‘rigid” scale. If 
the nebular nucleus is supposed to be moving at the rate suggested by 
Hubble’s determination of the red shift—say 500 km/sec, 10° parsecs 
—then the first meter marking on the rope will be moving past the 


first meter marking on the “‘rigid”’ scale at the rate of 


500 X 103 XK 3.16 X 107 
“ meters per year , 
EG” 3.00 X 10°" ms 


or by one part in 2X10’ per year. It is sometimes supposed that the 
nebulae are receding (if they are receding at all) very fast. It appears 
from this illustration that, considering their distances, they are re- 
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ceding excessively slowly. In fact, we could rarely distinguish be- 
tween the supposed fixed length of the meter scale and the supposed 
expanding first segment of the expanding rope. But if so, how can we 
be sure that the supposed “‘rigid”’ meter is not “‘naturally”’ expanding 
with time? And again, if so, we could join the nebular nucleus to our 
own galaxy with a rigid rod without the latter’s snapping. But that 
would be the same thing as the nebula’s not being in receding motion 
at all. It is just when we come to very small rates of motion like that 
of the nebulae that it is of supreme importance to be clear what we 
mean by a “rigid length scale.” It is as imperative to define the lat- 
ter as to define ‘‘uniform time.” But both are indefinable, conven- 
tional expressions. The only thing left to do is to link the two. 

We can do this by means of the to-and-fro transmission of light. 
Having chosen a mode of graduating our clock, we shall say that a 
given particle is distant a “constant distance” from ourselves by this 
clock, when the interval between the transmission of a light-signal and 
the reception of the reflection signal is constant, as recorded on this 
clock. It is to be observed that we assert nothing about the velocity 
of light. The velocity of light is, in fact, quite indefinable until we 
have specified a clock and given a definition of a distance. Thus, we 
are not assuming that the ‘“‘velocity of light is constant.’’ Nor do 
we assert that the forward and backward velocities of light are equal 

another old bogey in the supposed assumptions of relativity, for 
we have, as yet, no means of specifying the epoch of arrival of the 
signal at the distant particle—-we have so far constructed no “clock 
at a distance.”” We have, however, linked our distance measure with 
our clock graduation. 

We put this in symbols. A “clock” is, for us, just some arbi- 
trary correlation of the events of our experience with the real num- 
bers. Such a correlation is possible because the events of an observ- 
er’s experience form a well-ordered sequence—ordered by the ‘‘be- 
fore-and-after”’ relation. As Shakespeare says, ‘““We see which way 
the stream of time doth run.’’s (I will remark at this stage that there 
is no actual need to restrict the correlation to the real numbers. Any 
linear continuum of numbers, real or complex, will serve to order the 
well-ordered sequence of events in the observer's consciousness, and 


3 Henry IV, Part IT, Act IV, scene 1. 
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it actually occurs in certain cases of regraduation of clocks that it is 
convenient to adopt a complex time scale.) Our clock, then, is at 
ourselves and records only the epochs of events at ourselves. If, 
therefore, /, is the epoch of the transmission of a light-signal to a dis- 
tant particle, and ¢; the time of return, we put 3c (¢; — ¢,), where c 
is an arbitrary constant, as the distance of the distant particle by 
this clock, i.e., by this mode of clock graduation. We do not, of 
course, know or care whether this assessment of ‘‘distance”’ agrees 
with any other measure of distance that may be made by some other 
method; it suffices that we have a measure of the separation of the 
distant particle from ourselves at a certain definite instant. What 
instant? The instant of reflection. How are we to give a measure to 
this? Out of our two independent observations /, and ¢; we can con- 
struct two equivalent data in an infinite variety of ways. One we 
have already constructed, namely, 3c(t; — ¢,). Let us adopt for the 
other 3(¢; + ¢,) and call this the epoch, by the same clock, of the 
event of reflection. We put 7, = 3c (t; — 4), t. = 3(t; +h). 

It is clear that the measure of epoch-at-a-distance, ¢,, equal to 
> (t; + ¢,), is largely arbitrary. We might equally well have chosen 
(t,t,)"/? or [5 (8 + #)|"”’, etc. This is a familiar topic in relativity 

the conventionality of assignment of an epoch-at-a-distance. But 
this shows up at the same time, strikingly, that the assignment of 
distance-at-an-epoch is equally conventional. Instead of taking 
+c (t; — t,), we could take $c t, log (t,/t,) or 3c (2 — @)'”?, etc. Thisis 
a less familiar topic in relativity. General relativity is accustomed to 
regard all co-ordinates and combinations of co-ordinates as on the 
same footing. It is then prepared to take any function of /, and 4, 
as a co-ordinate defining ‘‘epoch-at-a-distance.”’ In particular, we 
may replace ¢,; and ?, themselves by 7; and 7,, where 7, = x(¢;), 
7, = x(t,), and the 7’s denote clock-readings regraduated. How 
then, if all co-ordinates are on the same footing, does it succeed in 
picking out one mode of clock graduation as unique? It does so by 
just the same process as it uses when it claims to pick out one mode 
of distance measure as unique—in both cases by refusing to define, 
owing to the impossibility of defining, a uniform “clock” or a “rigid”’ 
scale. The result is that it cannot be sure that one of the co-ordinates 
it employs is identifiable with the clock time of physics, or that an- 
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other is identifiable with the supposed rigid “length measure”’ of 
physics. 

I do not propose to contest ground, inch by inch, with general 
relativity. (It makes so many disconnected assumptions that it can 
usually extricate itself from the inadequacy of any one of these by 
relying on some other.) Instead, I hope I have succeeded in showing 
how fundamental in a discussion of the cosmological problem is the 
fixation of a distant particle’s state of rest or motion relative to our- 
selves, and that the description—1.e., “as a state of relative motion”’ 
or ‘‘as a state of relative rest’’—depends partly on the way our clocks 
are graduated and partly on conventions. 

Closely connected with the problems of setting up a clock at our- 
selves and of choosing our conventions for epoch-at-a-distance and 
distance-at-an-epoch is the problem of setting up a clock B at a dis- 
tance which shall be a copy of the clock A set up at ourselves. (We 
shall see later why this has an important bearing on the cosmological 
question.) The clock at a distance B is supposed to be represented 
by the again arbitrary graduation of the temporal experience of the 
observer at B. This observer is identifiable with some particle which 
may be “moving” in any manner, its motion being assessed by our 
clock A, by specifying the relation of its distance-at-an-epoch to the 
epoch concerned. How can A dictate instructions to B, enabling B 
to graduate his clock (or, if already graduated, to regraduate it) so 
that B’s clock may be said to keep the same time as A’s? Clearly, on 
the possibility of doing this rests the whole significance of speaking 
of ‘‘time”’ at one place in relation to “time”’ at another. 

If A’s and B’s experiences are completely uncorrelated, no mean- 
ing can be attached to asking whether they both “keep the same 
time,”’ for their experiences would be independent. Their experiences 
may be correlated, in principle, merely by allowing them to view one 
another—in addition, that is, to allowing them to assess one anoth- 
er’s motion by to-and-fro light-signals. Now, if there is any asym- 
metry in the relation holding between A and B, we could not deny 
the possibility of an asymmetry in their description of one another, 
each using his own clock; and we should not know how much of this 
asymmetry in description was due to an asymmetry in their respective 
methods of time-keeping. If, on the other hand, the relation between 
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A and B is a symmetrical one, it must be possible for them to gradu- 
ate their clocks so that the symmetry extends to their descriptions of 
one another when made by each using his own clock and viewing the 
other’s. When such graduation is achieved, they may be said to be 
“keeping the same time.”’ More precisely, it can be shown that A, 
having arbitrarily graduated his own clock, can dictate instructions 
to B (derived by A from his observations of B) so that B can re- 
graduate his clock to keep the same time as A’s. The clocks are then 
said to be congruent, the observer’s equivalent. It can be shown 
further that, if at any time during their relative motion, A and B 
coincided momentarily in position, then B’s regraduation of his clock 
can be uniquely determined. (This is connected with the circum 
stance that, when the two observers coincide, their clocks must read 
the same epoch. But it is of a much more subtle character than this 
circumstance alone would appear to indicate.)* If A and B never 
coincide, there remains an arbitrariness in B’s regraduation. 

Now that A and B keep the same time, we may inquire under what 
circumstances a third observer, C, keeps the same time as both A and 
B. It is convenient to confine attention first to observer C collinear 
with A and B. Then, to state the problem precisely, we may ask 
whether, when A has dictated instructions to B and to C, so as to 
put them in possession of clocks congruent to A’s, B and C will find 
themselves in possession of clocks congruent to one another. The 
answer is, in general, that they will not. C’s clock, in these circum- 
stances, will be congruent to B’s only if there is a certain restriction 
placed on C’s motion—his motion, that is, relative to A or B. It can 
be shown further that when C satisfies this condition, if a fourth 
collinear observer D is considered, and it is asked that he shall have a 
clock congruent to A’s, B’s, and C’s, then the motion of D, and 
a fortiori of all other collinear observers FE, F,...., is completely 
determinate. Such a collinearity of observers possessing congruent 
clocks is called a “linear equivalence.’ Its study is fundamental in 
all questions of time-keeping and in the general problem of the motion 
structure of the universe. 

A property of a linear equivalence is that, if two members ever 

4 See a paper by M. Crum communicated to the Quarterly Journal of Mathematics 
(Oxford). (10, 155, 1939.) 
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coincide in position at some epoch, then all members of the same 
equivalence coincide in the same position at the same epoch. But 
there need be no epoch of coincidence. Examples of linear equiva- 
lences are: (1) aset of stationary particles, (2) a set of particles in uni- 
form relative motion which all coincided at one epoch; (3) all kinds 
of accelerated motions of particles starting from coincidence, among 
which is the equivalence of particles in so-called “uniform’”’ accelera- 
tion considered by Leigh Page; (4) all kinds of accelerated motions of 
particles which never coincide.* 

But, however much these types of motion may appear to differ, they 
are essentially one and the same entity, in that by appropriately re- 
graduating any one observer’s clock in any one linear equivalence 
according to some law T = x(t), and then regraduating all other 
observers’ clocks according to the same law, this linear equivalence 
may be reduced to any other selected linear equivalence. That is to 
say, the descriptions of the motions of the members of the first 
equivalence as seen from any one member and determined by means 
of his clock can be brought into coincidence with the description by 
any other member of any other equivalence, by suitable clock re- 
graduation. It may occur that in certain cases the regraduation is 
complex, but this does not affect the principle. 

By suitable means the linear equivalences may now be generalized 
to give three-dimensional equivalences, which fill three-dimensional 
space. We then have, as examples, the relatively stationary equiva- 
lence, consisting of a three-dimensional set of particles at relative rest 
and a set of particles in uniform relative motion expanding outward 
from any one of them, considered as origin. We can also have any 
variety of accelerated motion. In accordance with the general prin- 
ciple, these can all be reduced to one another by clock regraduation. 
For example, if we regraduate the clocks of the uniform motion 
equivalence, reading time ¢, by a logarithmic regraduation to time 
tT, where rt = log ¢, or, better, r = ¢, log (t/t.) + t., where ¢, is an 
arbitrary constant, then the uniform motion equivalence turns into 
the one at relative rest. 

There is, in accordance with expectations mentioned earlier and 

‘It has been shown by Whitrow that the general linear equivalence is character- 


ized by a generating function y, such that ¢ Oan (t:) = Waap Y(t), where #2 is B’s 
time of reception of the signal leaving A at A’s time ¢,, and agg is a number. 
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based on general principles, a natural time-zero for the uniform mo- 
tion equivalence, which may be taken to be the epoch of coincidence 
of all the members. In accordance with the same general principle, 
the relatively stationary equivalence, being in a steady state, can 
have no zero of time. We notice at once, in confirmation, that t = o 
corresponds to r = —. By adjustment of the arbitrary constant 
t., 7 can be made to take any value for any given value of ¢ and, in 
particular, to agree with ¢ when ¢ takes the value fy. 

An equivalence is simply a distribution of motions, the particles 
executing which have a common measure of time and carry congru- 
ent clocks. For a world-wide system of congruent time-keeping it ap- 
pears necessary that the clocks shall form members of an equiva- 
lence. The fundamental problem therefore propounds itself: Are 
the extragalactic nebulae, considered as particles, members of an 
equivalence? 

If they are, several other questions answer themselves at once. 
For example, though an equivalence may be described in different 
ways, according to its mode of clock graduation, there is essentially 
only one equivalence. Hence, if the universe is of the motion type of 
an equivalence, there is only one type of universe possible. This is at 
once an immense improvement on the answer given by the older, 
relativistic (first-method) theories, which, as we saw, led to different 
possible types of universes, among which they could not discrimi- 
nate. We have, instead, different possible descriptions of one and the 
same underlying kinematic entity. The allegedly different motion 
types of the older theories must be due either to their failure to dis- 
criminate between the effects of different types of clock graduation 
in altering the description of a system or in their admission of possi- 
bilities, such as that a light-signal may be received back at two dis- 
tinct epochs, which we have not admitted. Such possibilities, which 
lead to such profound difficulties on further analysis, may be dis- 
carded; they have led to such fantastic results as that a particle may 
affect the motion of another particle the more the farther away it is 

results repugnant to common sense and easily seen to derive them- 
selves as mathematical rather than physical possibilities. It is highly 
satisfactory that, when we approach the problem of the structure 
of the universe by the second method, we are led to only one eventual 
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possibility, different as are the appearances which this possibility 
may take. It has been sometimes claimed as an accusation of lack of 
generality in the present treatment that it does not provide these 
tiresome choices of k = 0, +1, —1, which relativists so often dwell 
on. But we have been led to a single possibility inevitably and not 
by ignoring choices of possibilities presented to us. After all, there 
is only one universe; and a satisfactory analysis of the possibilities 
of universe structures should lead to only one type of such structure. 

Again, the objection has been raised to the present solution that it 
is “‘special,’’ in the sense in which the restricted theory of relativity 
is special compared with the so-called ‘“‘general”’ theory; that it con- 
templates only uniform motion, and handles it only by Lorentz for- 
mulae. But this does not impugn the generality of the system dis- 
cussed, for, according to the general theorem, we can reduce any 
motion type to uniform relative motion, any equivalence to the 
equivalence in uniform relative motion. We can, therefore, handle 
the most general type of motion, involving accelerated relative mo- 
tion of any character, by suitable regraduation of clocks; and the 
formulae of transformation among members of any equivalence, 
easily obtained in their general form, pass, by suitable regraduation 
of clocks, into Lorentz formulae. Actually the present analysis of 
the cosmological problem is perfectly general; it is simply that I 
choose a particular idiom in which to describe it, instead of using 
a characterless general language. 

If the universe is representable as an equivalence, we see at once 
that the question often asked, “Is the universe expanding or not?”’ 
is a no-question. For an equivalence can equally well be regarded as 
relatively stationary or as in any kind of relative motion; in the for- 
mer case there is no epoch of coincidence of all the members, in the 
latter case there may be. Simply by regraduating clocks, the same 
equivalence may be represented as of the one character or the other. 
Thus, whether the nebulae all coincided at one epoch or not is an- 
other no-question. Certain modes of clock graduation lead to special 
features in the resulting description. For example, the only mode of 
description which contains no parameter /, or its equivalent, in itself, 
is that corresponding to uniform relative motion. The only one con- 
taining no natural zero of time is the relatively stationary one; and 
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in this the epoch assigned by any observer to an event is the same 
whatever the observer. (It should be particularly noted that, when 
we are considering equivalences, the relatively stationary equiva- 
lence is not a particular case of the uniform relative-motion equiva- 
lence.) 

It thus becomes necessary to ascertain observationally whether 
the system of the extragalactic nebulae constitutes an equivalence 
whether, that is, there is a world-wide possibility of the construction 
of congruent clocks. This is a question logically quite anterior to the 
question which has so preoccupied the minds and time of astrono- 
mers, namely, whether the universe is homogeneous or not. For it 
should be carefully noticed that in the specification of an equivalence 
no density distribution is concerned; density of members is entirely 
irrelevant. The specifications of an equivalence involve simply the 
motion of each individual member relative to some given member 
the relation, that is, of distance-at-an-epoch to epoch for that mem- 
ber. We want to know, for the extragalactic nebulae, whether they 
belong to the equivalence defined by any three of them. 

If we could dispatch to and receive light-signals from them, these 
would be, in principle, a simple but slow matter. Actually we can only 
receive signals. And the building-up of these signals into the conven- 
tional numbers called “distances” and “‘velocities” is a very difficult 
business if we are to relate them to fundamental measures made with 
a given clock. The methods of determining nebular distances in- 
volve as an intermediary the construction of a good deal of geometry 
and a good deal of physics; and the further discussion of motion in 
volves the differentiation of these distances, when ascertained, with 
time. The observation of the red shift by itself tells us nothing save 
that, to a first approximation, atoms keep kinematic time the time 
of a clock graduated to represent the equivalence as being in uniform 
relative motion; and any equivalence obeys to a first approximation 
a linear relation between distance and velocity. One of the greatest 
of theoretical needs is the determination of what a Mount Wilson or 
Harvard distance means on the fundamental time scale it must be 
tacitly based on. 

And one of the greatest of observational needs is the detection of 
absolute change in the geometrical appearances of extragalactic ob- 
jects with time. It seems to me that here is one of the principal tasks 
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of large telescopes. The most distant visible objects, increasingly 
dimmed by recession, but more and more numerous, will have enor- 
mous red shifts. But they will be at a finite distance. They will 
therefore give the maximum chance of yielding changes of angular 
diameter with time, for example. I have not the slightest doubt that 
no such changes will ever be detected, even if observationally de- 
tectable in principle. For, if the universe can be assimilated to a 
three-dimensional equivalence, it can be regarded as a stationary 
universe, in which an angular diameter remains constant. The inter- 
pretation of this on the kinematic scale, i.e., on the scale correspond- 
ing to uniform relative motion, is that the linear dimensions of the 
nebulae, on that scale, are expanding in proportion to their dis- 
tances. The same consideration applies to a so-called “rigid” length 
scale. For, if one is supposed erected at a nebular nucleus, perpendic- 
ular to the line of sight from ourselves, it must subtend a constant 
angle in a stationary universe on the r-scale, and therefore also on 
the /-scale; hence, since its distance away is increasing on the ¢-scale, 
its measure of length on the /-scale is increasing. This again shows 
that, if a nebular nucleus is supposed joined to ourselves by a chain 
of “rigid” rods, it will not part company from the end of the chain on 
the r-scale, and so not on the ¢-scale; and hence the rods must be ex- 
panding uniformly on the ¢-scale. 

The fact that the red shift is proportional to distance is, as re- 
marked, evidence that the nebulae compose an equivalence and that 
atoms keep approximately kinematic time as regards the frequencies 
they emit. It does not tell us whether the orbital frequencies of the 
atoms are constant in any particular scale. That question takes us to 
dynamics, and with it to the distribution of density in the universe 
an important question, but one coming second in logical order. Be- 
fore, however, we can discuss atomic dynamics, we must have a 
dynamics capable of general applications. Such a dynamics must 
be based on the motion of a “‘free’’ particle, for a force can be meas- 
ured only in terms of departure from ‘“‘free’”’ motion. Examples of 
‘free’ particles are afforded by the nebular nuclei themselves. So 
far we have considered these “‘motions”’ or states of rest kinemati- 
cally, i.e., as something given. We can now go on to ask whether an 
equivalence represents a natural motion, i.e., whether the members 
of an equivalence can of themselves, without an “external” agency, 
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pursue the trajectories provided for them in the specifications of an 
equivalence. To solve this problem it is first necessary to consider 
the motion of a free particle in general in the presence of an equiva- 
lence. And we choose, for the scale of time, i.e., the mode of clock 
graduation in which to conduct our analysis, the ¢-scale, the one in 
which the members of the equivalence are in uniform relative motion. 

Can we assert anything about the motion of such a free particle? 
Clearly, all we can do is to pass from one member of the equivalence, 
considered as point of observation, to another—from one “funda- 
mental” particle to another. We can transform the equation of mo- 
tion from one origin to another. The new form which this equation 
of motion takes will depend on the distribution of the totality of 
particles constituting the equivalence. For example, the equation of 
motion of a free particle in a dense part of the equivalence might be 
of a different form from that of a free particle in a less dense part of 
the equivalence. For we must remember that, when we speak of a 
“free” particle, we mean a particle in the presence of the totality of 
matter in the universe, or, as we should normally say, a particle ex- 
posed to the “gravitational pull” of the rest of the universe. We can 
assert nothing whatever about the motion of a particle “ignoring” 
this effect of the rest of the universe, for the rest of the universe is 
necessarily present—without it there would be no means even of de- 
scribing or locating the “‘test’”’ particle. Thus, the dynamics of a 
“free” test-particle must depend on the density distribution in the 
equivalence. Until the latter is specified, the “dynamics of a free 
particle’ has no meaning. This density distribution of the equiva- 
lence is something quite different from the motion pattern which de- 
fines the nature of the equivalence in the chosen scale of time. For 
the same motion pattern we can have kinematically any kind of 
density distribution. The specifications are of a quite different char- 
acter, and the specification of the motion pattern comes logically 
before that of the density distribution. The question so often dis- 
cussed, ‘‘Is the universe homogeneous?”’ is logically posterior to the 
question which ought to be asked first, “Is the universe an equiva- 
lence?” 

Let us then suppose we have specified our equivalence and speci- 
fied a density distribution in it. Two questions now arise: First, can 
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we determine from this the equation of motion of a free test-particle 
in the presence of this density distribution? Second, do the given 
particles of the equivalence obey this law of motion; i.e., is the equiv- 
alence, when kinematically specified, also a dynamical system itself? 

These questions have not yet been answered save in one important 
case, namely, that in which the equivalence is a substratum. A sub- 
stratum is an equivalence whose density distribution satisfies the 
cosmological principle; that is, if A and B are any two members of 
the equivalence, then the statistical description of the whole equiva- 
lence, from A, in its temporal unfolding, coincides with that from B. 
This we write A = B, a statement to be carefully distinguished from 
A = B, which means simply that A and B have congruent clocks. 
When A = B, for any two members A and B, the system constitutes 
an equivalence; when, in addition, A = B, it constitutes a sub- 
stratum. The ideas contained in the assertion A = B are a generali- 
zation of the ideas contained in the notion of homogeneity. Whether 
a substratum appears homogeneous or not depends on the scale of 
time adopted; for example, on the z-scale, or the scale in which the 
equivalence appears as at relative rest, a substratum is a homogene- 
ous distribution of fundamental particles; but on the ¢scale the den- 
sity distribution will appear to increase outward from any arbitrary 
observer, although no observer is in a favorable position—a remark 
which, I think, removes some of the difficulties stated in the 1936 
Rhodes lectures by Hubble. 

It is clearly necessary to solve this, the fundamental problem of 
dynamics—the determination of the motion of a free particle—for an 
equivalence satisfying the cosmological principle before attempting 
it for any other. This must afford a norm of knowledge, a standard 
of expectation. When we attack this problem analytically, when we 
attempt to find the equation of motion of a free particle in the pres- 
ence of a substratum, we find two highly significant results. 

The first is that the answer appears, at first sight, to be indeter- 
minate. The second is that, in spite of the indeterminacy, the kine- 
matically prescribed motions of the particle members of the equiva- 
lence satisfy the dynamical equation of motion. 

Choose the scale of time in which the equivalence (A = B) is in 
uniform relative motion, that is, the /-scale. Then the motion pat- 
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tern is V = P/t, where V is the velocity of the fundamental particle 
at P at time ¢. The corresponding substratum has a density dis- 


tribution 
Bt dx dy dz 


le- 3) 
Aah = = 
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provided that, in addition to the cosmological principle (A = B), a 
certain other principle, called the ‘‘principle of symmetry,” is sat- 


ndx dy dz = (1) 


isfied. 
The equation of motion of the particle passing through P at time / 
with velocity V is 


dV 
-=(P— Vi 
dt ¥ 


G (&), (2) 


where Y =1 — V2/?, X = # — P?/¢?, Z=it-—P-V/c?, andé = 
Z?, XY. The function G(&) is not determined by the two conditions 
A = BandA = B.° But whatever the form of G(£), the fundamen- 
tal particles moving according to V = P/t¢ satisfy the equation of 
motion (2). 

This is a result of fundamental importance, whose interpretation 
has, I think, not so far been given. It means, I think, that the 
kinematical system called an ‘‘equivalence” can only be also a dy- 
namical system when it is also a substratum satisfying A = B. For 
suppose that we now alter the density distribution so that it no 
longer obeys the cosmological principle, but that we leave the motion 
pattern unchanged. Then we shall certainly alter the equation of 
motion of a free particle. This alteration will not be merely a change 
in the function G(é), as we shall see ina moment. It must, therefore, 
be a more drastic change, and the law V = P ¢ will no longer satis- 
fy (2). Hence, the equivalence, as specified kinematically, will not 
conserve itself dynamically and must alter in its state of motion, 
not necessarily even remaining an equivalence. The substratum, or, 
as we may Call it loosely, the ‘homogeneous equivalence,” is there- 
fore a fundamentally important type of equivalence in that it alone 
conserves its character of being an equivalence. 

6 The reader who is unfamiliar with the mathematical background may be informed 
that A = B determines the transformation from one fundamental particle to another, 


while A = B determines the conservation of form of the equation of motion and the 
density distribution. 
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There are, of course, many difficult, unsolved problems to be in- 
vestigated in this region. We should like to know, for example, the 
conditions of stability of an equivalence arranged as a substratum; 
this would require the knowledge of the equation of motion of a free 
test-particle in the presence of an equivalence which was not a sub- 
stratum—again an unsolved problem. But one remark may be made 
before we pass on. It is commonly supposed that there is such a 
thing as a “dynamics,” that is, a law of nature giving the equation of 
motion of a free particle in a form independent of epoch and location 

a law of nature satisfying the principle of the uniformity of nature. 
My point is that those who believe this are tacitly using the cosmo- 
logical principle. If the universe is not homogeneous (in our extend- 
ed usage of the word), then the form of the dynamics of a free particle 
will depend on locality, and there will be no world-wide dynamics 
possible. There is not sufficient observational evidence to say wheth- 
er a world-wide dynamics is necessary or not; but it is inconsistent 
to appeal to, or to assume, the principle of the uniformity of nature 
and yet to deny the cosmological principle as soon as it has been 
shown that the forms of the laws of nature depend on the distribu- 
tion of matter in the universe and that these forms can be the same 
everywhere only when the distribution of matter is describable in 
the same terms from every fundamental particle as observing-point. 

The next thing to notice is the apparent indeterminacy contained 
in the function G. It may be remarked that the function G in the 
form G(£) has been reduced from the more general form G(X, &), 
from the consideration that the substratum, as specified by (1), con- 
tains no constant of the dimension of a time in its description; and 
as Y could occur only in the form X/#2, where ¢, was a constant, it 
cannot occur at all.?7 Thus G = G(§). 

Now, why is G thus indeterminate? The substratum is indetermi- 
nate only to a constant B, and this is equivalent only to the arbi- 
trariness in reckoning what is a fundamental particle: it has no 
dynamical significance. Hence, the equation of motion must be 
completely determinate—a free particle must move in an absolutely 
foreseeable way. 


On the contrary, when we use dynamical time 7, the density distribution (1) be- 
comes replaced by B cit}; and a constant fo appears and may reappear (as it does) in 


the equation of motion of a free particle 
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NO 


The contrary has been urged by H. P. Robertson in a passage 
which expresses the view that the universe is, at bottom, magica! and 
irrational. Thus he remarks: ‘The cosmological principle reduces 
the problem of motion to that of determining, by the imposition of 
physical law [italics mine], a single scalar function.”’ The view that 
we have to impose a physical law, i.e., something extraneous to the 
considerations before us, is entirely foreign to the whole set of ideas 
here being developed. Robertson, in fact, misconstrued the force of 
the kinematic argument. 

The reason that G is indeterminate from the considerations so 
far advanced is not that G requires some extraneous “‘physical”’ law 
to fix it but that the system defined by specifying that it shall satisfy 
the cosmological principle is still itself indeterminate. A substratum, 
i.e., a homogeneous equivalence, is a system of hydrodynamical char- 
acter; but far more general systems will have an equation of motion 
of a free particle of precisely the same form, simply because they 
satisfy the cosmological principle A = B. The indeterminacy is 
G(é) corresponds to the possibility of defining vastly more of these 
general systems. Since there is essentially but a single hydrodynami- 
cal system, more general systems cannot be of hydrodynamical char- 
acter but must be of a statistical character. And as soon as we in- 
vestigate the possible diversity of such statistical systems, we find 
that there is just the degree of freedom corresponding to the inde- 
terminacy of G; namely, a single function (&), itself arbitrary, de- 
scribes each statistical system. Its density-velocity distribution is 
found to be, in fact, of the form 

_ Wl&) dx dy dz du dv dw 


n dx dy dz du dv dw = —— OX Ps ; (3) 


If this exhausts the possibility of obtaining systems satisfying the 
cosmological principle, then we should expect that there should be a 
relation between ¥(£) and G(é). And we find, in fact, from argu- 
ments equivalent to the mere counting of particles, that the relation 
is just 
+ 
OS) =F aay (4) 
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where C is again arbitrary. This relation is quite unlike any other 
relation in the whole range of mathematical physics, for without any 
“imposition of physical laws”’ it derives by purely a priori methods a 
relation between density and motion—in other words, a relation of 
purely gravitational character is derived without recourse to experi- 
ment. 

We have now not only arrived at the reason why G(§) is, at first 
sight, indeterminate; we have obtained an actual expression for 
G(£) for the general case of a statistical or nonhydronomical distribu- 
tion. By limiting processes, of which the details have been given 
elsewhere,* it can be shown that the first term, — 1, in (4) represents 
the contribution to G(é) of the hydrodynamical system by itself. 
This, therefore, determines G(£) for a pure substratum as 


G(é) = —I, 


which is equivalent to the purely logical deduction of the first law of 
motion. It is the fundamental theorem of dynamics. For, on trans- 
formation from the scale of ¢ to the scale of 7, the equation 


dV _ 


V 
i —(P — Vi) > 


becomes an equation stating that v, the r-velocity of a free particle 
in a certain public hyperbolic space, remains constant. It remains 
constant in spite of the filled-up-ness of the substratum or universe 
with matter, and allows fully for the gravitational effect of this mat- 
ter on the particle. This shows that in a dynamics based on the as- 
sumption of Newton’s first law of motion, the time variable must be 
what we have called 7. Our procedure is something quite different 
from the procedure followed in general relativity. There, the metric 
depends on the filled-up-ness of the space with matter, and the paths 
are assumed to be geodesies; in our procedure, the space chosen is 
independent of the density of the matter, and, in the measures, the 
free paths are not geodesies.? 

As soon as the motion of a free particle is known, a dynamics can 
be constructed. This has been done in papers published elsewhere. 

§ Ouart. J. Math. (Oxford), 8, 22, 1937. 


In ¢-measure each observer has his own private space. 
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The principal result of interest is that mass, inertial mass, is repre- 
sented in the ¢-scale by mé"/?; energy, by mc’, where 


V 
m( i — —- a 
rm 


= = (1 - y? ( ‘oe =) he 
ie 2 
"ig Cc’, 


When P = o, i.e., for an observer momentarily coinciding with 


the particle, this reduces to m/(1 — V*/c?)'?, which is Ein- 
stein’s expression, save that here V denotes specifically not the 
velocity relative to any standard of rest, as in the special theory of 
relativity, but the velocity relative to the local standard of rest. It 
is quite clear that in all dynamical calculations, such as, for example, 
the advance of perihelion of a planet, the local standard of rest—the 
standard fixed by the extragalactic nebulae as a whole—is destined 
to play a much more important part than that attributed to it in the 
theory of relativity. 
If, on the other hand, an extragalactic observer O is chosen so that 
V = o, then the expression reduces to m, (1 — P?/c’l?)'’? and meas- 
ures the “‘potential” energy of the now ex-centric particle. 
So much for the —1 in formula (4). What of the other term? This 
has an equally exciting interpretation. Since 
Eman tle ey - FA 


om) Cc? 
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the acceleration term in (2) corresponding to the second term in (4) 
reduces, under well-detined circumstances, to 


CBC Bshaa Vt 
¥(1) (P — Vt)3’ 


which is an inverse-square attraction toward the point Vé, which is 
the center of the neighboring condensation in the statistical system. 
The essential point is that, without at any point making any as- 
sumption corresponding to the inverse-square law, an inverse-square 
law has emerged. It is, of course, the law of local gravitation. The 
curious index 3/2 leads precisely to the negative exponent 2 in the 
distance of the particle from the condensation. It must be empha- 
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sized that the arguments employed in deducing (4) are purely 
kinematic; yet, an inverse-square law has appeared. 

By appropriate analysis, this formula can be modified to give the 
gravitational potential x for any pair of particles in the form 


ee { eee 
ee M, P, . P, 2 P? é P? /2? 
Ps i C2 


where M, is the mass of the equivalent homogeneous universe and /, 
and t, are the epochs at the two particles. Here we have introduced 
the masses m, and m, of the two particles, and this leads us to refer 


(5) 


for a moment to the arbitrary constant C in formula (4). 

Why should this formula contain an arbitrary constant? Given 
the function ¥(é), determining the population of the statistical sys- 
tem, it appears that there yet remains a degree of arbitrariness in the 
acceleration of a free particle. Why is this? The answer is that C 
corresponds to the mass of the condensation; C is thrown up by 
purely kinematic arguments, which are thus capable of predicting 
deductively the existence of gravitational mass. The actual con- 
nections are given by 

M, = 4xnm.B, 
C 


M,, 
y(1) ; 


where m, is the arbitrary mass number attached to a fundamental 
particle. The details are of less interest than the unexpected fact 
that purely kinematic arguments insist on the appearance of a new 
parameter which can be identified with what is known empirically as 
“sravitational mass.’ This gives a kinematic basis to gravitation. 
We are ultimately led to the relation y = c*t/M,., which allows a 
determination of M/, from ¢, the age of the universe, and y, the ‘‘con- 
stant” of gravitation. 

Thus, the kinematic arguments leading to (4) establish: (1) the 
law of inertia for a substratum; (2) the existence of ‘‘gravitational”’ 
mass; (3) the inverse-square character of “‘gravitational’’ force. 
These are thus removed from the domain of empirical knowledge to 
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the domain of deductive knowledge, though, of course, the appeal to 
observation remains in the act of identifying in nature entities cor- 
responding to those thrown up deductively. 

I have not time here to show in detail how these ideas can be en- 
larged so as to include electromagnetic forces. It suffices to remark 
that a “gravitational force” is derived by the differentiation of the 
potential (5), which is the simplest scalar potential which can be 
built up possessing a singularity at P, = P., t, = t,. That is to say, 
a gravitational force is the simplest force, arising in a substratum, 
that can ‘‘act”’ on a free particle. If we inquire as to the next sim- 
plest type of ‘‘force,”’ it appears that it arises from the double differ- 
entiation of super potentials, scalars which also possess singularities 
at P, = P,and have certain other properties. These ‘‘next simplest” 
types of forces then prove to be identifiable with electromagnetic 
forces. The further study of these, as between neighboring ‘‘charged”’ 
particles, leads to results of some interest in the theory of nuclear 
dynamics. It may be remarked that, just as “gravitational mass’”’ 
is simply the numerical label distinguishing the simplest type of sin- 
gularity possible, so ‘‘charge”’ is the numerical label distinguishing 
the next simplest type of singularity. 

Much of what I have been saying is not in accordance with certain 
current ways of thinking. But it appears to me that the experimen- 
tail and observational method—marvelous, as it is, for discovery—~is 
powerless to indicate the answers to the deeper philosophical ques 
tions of the origin of laws of nature—to tell us why they hold good. 
This requires a purely logical analysis. Actually, even a logical anal- 
ysis cannot tell us why they hold good, any more than a geometry 
can tell us why theorems hold good. But just as a geometry can tell 
us what geometrical theorems hold good, and give us a degree of be- 
lief in their validity deeper than can be gained from the principle of 
induction, so kinematical analysis (geometry and time) can tell us 
what kinematical, dynamical, gravitational, and electromagnetic 
theorems must hold good. For example, an appeal to the principle of 
induction cannot distinguish between an inverse-square law 1/7? and 
a law 1/r°*” ’; and the so-called ‘“‘appeals to simplicity’? urged by 
Jeffreys carry little conviction. Kinematic analysis leads to the deti- 
nite law 1/7’, in its relativistic form, for local gravitation, thus de- 
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riving a so-called “‘law of nature’’ purely from the self-consistency of 
all the elements of observation with which that law deals. The de- 
tails of this analysis are unlike the usual type of physical argument, 
and their very unfamiliarity tends to make them suspect. I am con- 
tent to let the march of time bring that conviction of their reliability 
which has come to all who have studied them without preconceived 
prejudices. 

In conclusion, I may mention that kinematic analysis gives defi- 
nite answers to questions which can be asked about the universe but 
to which empirical relativity gives ambiguous answers. Among these 
are the questions: 


1. Is the universe expanding or nonexpanding? 


iS) 


Is it finite in spatial extent or infinite? 

3. Is the number of particles in it finite or infinite? 
4. Is the space appropriate to it curved or flat? 

Is its past history finite? 


oa 


6. Is the universe homogeneous or not? 


Kinematic analysis says that the answers to these questions de- 
pend on the scale of time adopted and that the adoption of a scale of 
time is an arbitrary act on the part of the thinker. Thus, if the scale 
of time is that in which the laws of dynamics take their classical 
form (r7-scale), then the universe is not expanding, is infinite in spa- 
tial extent with an infinite past history, and the appropriate public 
space is hyperbolic. If the scale of time is that in which an atomic 
frequency is taken as constant (é-scale), then the universe is ex- 
panding uniformly, is finite in spatial extent with a finite past his- 
tory, and the appropriate private spaces are flat. In the former case 
the density distribution is homogeneous; in the latter case the den- 
sity increases outward from the observer—and, indeed, from every 
observer. In both cases the number of particles is infinite, but there 
is a parameter M, which plays the part of a finite mass for the uni- 
verse. 

These answers depend fundamentally on only one axiom, namely, 
that a meaning can be attached to saying that different observers in 
different parts of the (expanding or nonexpanding) universe can set 
up the same measures of time when a precise arithmetical meaning 
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is given to this phrase. The element of experience incorporated into 
the analysis is the individual observer’s consciousness of the passage 
of time, i.e., his awareness of a before-and-after relation between any 
pair of events he experiences. Starting from this awareness, and an 
initial arbitrary graduation of temporal experience, it proves possi- 
ble to isolate that measure of time, or graduation of temporal ex- 
perience, in which a free particle, projected among a “homogeneous” 
set of equivalent time-keepers, conserves its velocity. In this meas- 
ure of time the system is stationary, nonexpanding, and of infinite 
extent. The red shift is, therefore, evidence that atomic vibrations 
are not conserved in period in this particular measure of time. Fur- 
ther analysis shows that they must be conserved in a (in some ways) 
more fundamental measure of time, that in which to any particular 
observer the universe is expanding, finite in extent but still infinite 
in population, and nonhomogeneous, with density increasing out- 
ward from each individual observer. The physics appropriate to this 
measure of time contains no arbitrary constant but only convention- 
al constants like c and M,. The theorems of the analysis are recog- 
nizable as laws of nature, but from their logical origin they have a 
sanction exceeding that of any law of nature depending on the prin- 
ciple of induction and on the supposed principle of the uniformity 
of nature. The latter principle stands exploded, for in a universe 
which is a seat of change, where epoch has a meaning, the laws of 
nature themselves must be functions of epoch, reckoned from some 
natural zero. The spectator of the universe can take his choice. If 
he regards it as changing, epoch is significant; and the laws of nature 
will, in general, be expected to involve this epoch. If it is stationary, 
then dependence on epoch disappears; but the existence of the red 
shift then shows that atomic emission and absorption frequencies 
are not constant in the time scale chosen. The existence of the two 
time scales is closely connected with the disjunction between matter 
and light. The r-scale is appropriate to material particles and their 
dynamics and electrodynamics; the ¢-scale, to radiation. It is tempt- 
ing to say, in conclusion, that it is the interplay of these two time 
scales which gives the universe its interest. 
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ABSTRACT 


1. Apparatus.—The Mount Wilson observations were made with the double quartz 
monochromator (slightly modified) that was used in the Tucson work. Light was sup- 
plied by an aluminized siderostat mirror. Both lens-projection and mirror-projection 
systems were used. 

2. Measurements.—The measurements were made during the years 1934, 1937, and 
1939, With the thermopile, by the method of constant dispersion. The results, reduced 
to no atmosphere, yield spectral energy-curves for integrated sunlight and for light from 
the center of the disk which verify the principal features obtained at Tucson, namely, 
a sudden drop of 48 per cent in intensity from \ 0.40 to A 0.38 w, a nearly constant in- 
tensity from A 0.38 to \ 0.325 mw, and a nearly linear fall of intensity to a very low 
value at \ 0.292 py. 

3. Atmospheric transmission coe fficients—The means of determinations on 35 days, 
which include both integrated light and light from the center of the disk, do not differ 
more than 2 per cent from those of the Smithsonian observers until the ozone band is 
reached. Corrections due to the fall of energy within the second slit have been applied 
to the results for \ 0.295 and 4 0.292 yu. 

4. Atmospheric ozone content.—A plot of atmospheric transmission against wave 
length shows that the ozone content is about 0.10 cm at N.T.P., only slightly more than 
half that found at Tucson. 

5. Energy between the lines —To eliminate as far as possible the effects of Fraunhofer 
line absorption, energy-curves were obtained with a 21-foot concave-grating monochro- 
mator and photoelectric cell over 100 A units of the spectrum. The areas below contour 
lines drawn over the highest ordinates, divided by the areas below the record of the line 
spectrum, furnish coefficients by which the thermoelectric energy-curves may be re- 
duced to an energy-curve for the continuous spectrum. Although this procedure 
smooths out the ultraviolet energy-curve to some extent, most of the major depressions 
remain, and the curve in the ultraviolet is still far removed from that of a black body. 

6. Band spectrum of ozone.—An earlier laboratory investigation of the ozone spec- 
trum is reported. A calibration of the plates shows that the line intensities are so low 
that they will not be distinguishable in the solar spectrum and cannot alone account for 
depressions found in the photoelectric energy-curves. 

7. Line-group intensities —From the photoelectric records of the 21-foot monochro- 
mator the line-group intensities for the regions AA 3200-3300 A and AA 3900-4000 A are 
given. A comparison with Mulders’ calibration of Rowland intensities in the latter 
region shows an average ratio of 0.81, with considerable scattering. 


A determination of the distribution of energy in the ultraviolet 
solar spectrum made at Tucson, Arizona (elevation 2600 feet), in 
1931, has been described in a previous paper.’ As most of this work 
concerned the spectrum of integrated sunlight, it was thought useful 
to repeat the observations on Mount Wilson with special emphasis 

* Contributions from the Mount Wilson Observatory, Carnegie Institution of Wash- 
ington, No. 622. 

t Mt. W. Contr., No. 445; Ap. J., 75, 185, 1932. 
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on the radiation from the center of the disk, which, after all, best 
represents what the sun emits from the photosphere. 

The observations at Tucson were made with a double Bausch 
and Lomb quartz monochromator and a vacuum thermopile. Sun- 
light was supplied with a Foucault underdrive siderostat with a 
12-inch stellite mirror. The image-forming lens was of fused quartz, 
of 12 inches aperture and 6 feet focal length. 

In arranging the equipment for operation on Mount Wilson, the 
same monochromator and thermopile were used, but either a 6-inch 
fused quartz lens of 6 feet focal length or a mirror system replaced 
the lens used at Tucson. 

CONCERNING SIDEROSTATS 

Siderostats are of two types. That used at Tucson (Fig. 1, A) has 
the mirror set perpendicular to the driving rod R; the declination 
arm D connecting it swings in an arc beneath the mirror. This ar- 
rangement has the disadvantage that in the original adjustment of 
the instrument it is not easy to make the horizontal beam of light 
MN, projected back through the mirror, intersect the polar axis 
at the same point O with the declination arm. In the overdrive type 
(Fig. 1, B) adopted at Mount Wilson this difficulty is overcome by 
drilling the axis a of the mirror fork f, which now points in the beam 
direction mn. The ends of the drilled axis a are closed with glass 
plates pp’ having the centers marked. These act as sights. In 
setting up the siderostat for the initial run, the mirror is removed, 
and the axis sighted upon the optical part to receive the light. The 
mirror is then replaced, and the declination arm d, provided with a 
lens finder ss’, having been set upon the sun, is changed slightly in 
length until the beam from the mirror falls centrally upon the optical 
part desired. This procedure secures the most accurate following. 
The same roller slip-connection k between the declination arm d 
and the mirror driving rod r (now in the plane of the mirror) that 
was devised for Tucson has been adopted in the new instrument 
with equal success. An elevating screw e and an azimuth screw h 
facilitate the adjustment of the projected beam mm in both azimuth 
and altitude; the latter adjustment is mechanically difficult to effect 
with type A. The attachment z is an air-mass meter, 
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This siderostat, provided with a 12-inch aluminized mirror, was 
made part of the remounted 20-inch telescope and has electrical con- 
trols in both co-ordinates. The telescope, mounted on rails, can be 
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Fic. 1.—A, Underdrive siderostat as used at Tucson in the earlier measurement of the 


solar energy-curve in the ultraviolet. B, Overdrive siderostat as used on Mount Wilson. 


moved in the meridian a distance of 8 feet. The beam from the sid- 
erostat enters the building to the south where the double quartz 
monochromator and a 21-foot concave-grating monochromator are 


placed. 
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THE DOUBLE QUARTZ MONOCHROMATOR 

This consists of two monochromators with a common straight 
slit, the first and third slits being curved. The width of the second 
slit, which determines the wave-length range of the radiation meas- 
ured in a given setting, was read on a 5-inch dial graduated to 100 
parts. The four focusing pinions are now operated by a common 
knurled head. 

The procedure in making the measurements, with some modifica- 
tions, is the same as that already described.’ In half the work of 
1934 and in all that of 1937 an all-mirror projection system replaced 
the quartz lens, which eliminated the focusing of this part of the 
equipment for each wave length. This system involved two or three 
reflections from aluminized surfaces. In reducing the 1937 observa- 
tions it was found that, possibly owing to prolonged exposure to sun- 
light, the reflecting power of the mirrors to the violet of \ 0.37 u 
had fallen sharply. In 1939 this mirror system was therefore aban- 
doned in favor of the single mirror and fused quartz lens, an arrange- 
ment used at Tucson in 1931 and on Mount Wilson during a part of 
1934. The principle of eliminating all unnecessary reflecting surfaces 
in ultraviolet radiometric determinations seems to be well estab- 
lished by this experience. 

In beginning the observations on Mount Wilson it was found that, 
owing to the peculiar construction of the slits, the changes of tem- 
perature in the early summer seriously affected the slit widths. 
This was especially troublesome in the case of the first slit, since the 
zero of the second was easily tested. At Tucson the temperature 
range was small during the observing season and therefore did not 
produce much interference. In 1934 and 1937 the temperature of the 
observing room on Mount Wilson was in part controlled, and cor- 
recting formulae were devised to eliminate, in so far as possible, the 
outstanding effects. In 1939 the first slit was replaced by a fixed 
invar slit. The third slit was also removed, and the thermopile 
placed directly in the spectrum. The fused quartz transfer lens 
formerly used was thus eliminated. No photoelectric observations 
were made as in previous years, but the system, more sensitive than 


2 Mt. W. Contr., No. 445, p. 26; Ap. J., 75, 210, 1932. 
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that employed at Tucson, required a 5 per cent sector instead of 
the 10 per cent sector previously used, thereby making the deflec- 
tions in the extreme ultraviolet, when the sector was stopped, double 
those formerly observed. 

The effective length of the first slit was fixed by the length of the 
thermopile receivers, 4 mm. In 1934 and 1938 the diameter of the 
solar image was 16.8 mm so that the slit integrated the radiation 
over a region extending out from the center on either side to a dis- 
tance 0.24 of the solar radius. Since the darkening along this fraction 
of a radius is small for all wave lengths, the measurements of spectral 
energy refer efiectively to the center of the disk. In 1937 a solar 
diameter of 50 mm was used. 

THE OBSERVATIONS 

The program of observations was the same as that followed at 
Tucson. Since there were no considerable sky obstructions, the 
measures began early, when sec z=4 to 7. To enable the observer 
to read the air mass at any time, the air-mass meter’ was removed 
from the lunar equipment and attached to the declination arm of the 
siderostat (s, Fig. 1, B). This arrangement assists in securing the 
best distribution of points in the plots of the logarithms of deflec- 
tions against sec z 

The method of constant dispersion was used throughout. This 
procedure involved the following adjustments of the double mono- 
chromator: ‘The first slit was fixed at 1 division opening (0.1 mm, 
equal to 10 divisions on the second slit head); it is important that 
this opening remain undisturbed. The second slit was set for an 
opening which allowed 100 A of the spectrum of the first monochro- 
mator to pass through. The third slit was set to the same reading as 
the second plus one division. In the earlier work its image was trans- 
ferred to the thermopile with a quartz lens, but in 1939, as already 
noted, both slit and lens were eliminated, and the pile placed directly 
in the image formed by the second monochromator. The method of 
constant dispersion has the advantage that the deflections of the gal- 
vanometer, after correction for transmission, reflectivity, etc., are di- 
rectly a measure of the spectral energy. The method of constant slit 


3 E. Pettit, Mt. W. Contr., No. 504; Ap. J., 81, 17, 1935. 
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widths, on the other hand, involves dispersion corrections, which in 
the extreme ultraviolet become very large, more than 600 per cent 
at \ 0.30 w with a quartz monochromator. 

The wave lengths in which intensities were measured were the 
same as those used at Tucson and are indicated in Table 1. Each 
measurement consisted in the following operations: (1) The two 
wave-length dials of the monochromator were set. (2) With the aid 
of a table (second col., Table 1) on the observing sheet, the second 
slit was set to transmit 100 A. (3) Until its removal in 1939, the 
third slit was set to this value plus 1 division. (4) The monochroma- 
tors were focused. (5) The quartz projecting lens was focused. (6) 
A guiding image from the siderostat mirror was checked. (7) After 
the shutter was opened, the wave-length dial of the first monochro- 
mator was moved slightly by a lever slow-motion until maximum 
deflection was observed. (8) The hour angle was noted. All these 
operations were repeated for each of the 22 wave lengths. Further, 
the temperature and second-slit zero reading were noted at the be- 
ginning and end of each set. 

One such single set of measurements required about a half-hour 
to complete. The first two sets were run in immediate succession as 
soon as the observations began after sunrise. An interval of a half- 
hour was allowed between the second and third and between the 
third and fourth sets. A fifth and sometimes a sixth set was taken 
during the hour just preceding meridian passage. Usually only the 
wave lengths in the extreme ultraviolet were observed in the last set, 
since the intensities of the higher wave lengths showed little change. 

In all these measurements a rotating sector transmitting 5 per 
cent was employed until the deflections dropped to about 25 mm; 
then the sector was stopped, with the opening in front of the slit, 
thus increasing the deflections to 500 mm. This usually occurred 
near \ 0.33 uw for the first set, but near \ 0.30 uw in the last. No de- 
flections could be found at \ 0.30 in the first and at 0.295 and 
0.292 mw in the first two sets. This is illustrated in Table 1. It was 
customary therefore to add a few measures in this region as the sun 
approached meridian. In the work on integrated sunlight done in 
1934 this region was measured with the sodium photoelectric cell, 
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with an overlap to connect with the thermoelectric measures in a 
longer wave length. 
As an illustration, Table 1 gives the observing record for June 8, 
1939. Besides \ and the second-slit setting S,, the table gives the 
TABLE 1 


SOLAR ENERGY-CURVE OBSERVATIONS, JUNE 8, 1939 











(Morning) 
| 
| Set 1 Set 2 | Set 3 | Set 4 Set 5 
N | Sz | S, Zero= S, Zero= | S, Zero= | S2 Zero= S, Zero= 
| +0.05 +0.05 +0.06 +0.06 | +0.06 
T=22°5C | T=23°5C T=25°0C | T2623 C 
m | mm | mm | mm | mm | 
0.70 | 0.43] 613" 207 | sh40o™ 250) sho3™ 267) 3"00™ ses COE 
60 ©o.70} 6 II 281 5 38 358] 5 02 408] 2 59 GBB ees 
55 lo 89] 6 10 271 5 36 371] 5 00 428] 2 57 Ol ee 
50 | 1.17}609 273] 535 3721 4 58 4481 256 513] ... mm 
45.. 1.61} 6 08 245 | 5 34 366 4 56 445) 2 54 546] ob23™ 562 
42 | 2.00} 606 154] 5 32 2501 4 55 ana :2 63 407| © 22 429 
40.. 2.34) 0 05 125 | 5 31 224] 4 53 305] 2 §2 423} O 20 426 
39 2.54) 6 03 90 | 5 28 174) 4 siaT 240} 2 50 323] oO 20 332 
38 | 2.76] 6 02 57 |5 27 112] 4 08) 188] 2 49 )§6©224, 0 18¢ 245 
37 | 3.01] 6 00 50 | 5 25 105] 4 05 189! 2.47 2281017 249 
30 2.301 & co 4i | S24 QI] 4 04 172| 2 46 210) O15 228 
35 3.63] 5 58 x6oo*| 5 23 63] 4 03 139] 2 45 174) O I4 195 
34 | 4 02| 5 50 %*451 | 5 22 54] 4 02 I21} 2 44 155) oO 13 175 
33 | 4.45] 5 55 315 | 5 20 41} 4 OI IOI] 2 42 138] o II 156 
325 | 4 69] 5 54 *239 | 5 19 xXO70} 4 00 85} 2 40 118] 0 10d 132 
32 | 4 9o| 5 52 X141 | 5 17 430] 3 58 59] 2 39 86) 1 28e 104 
315 | 5.201 5 50 x 60] 5 15 4227] 3 57 37| 2 38 59| 1 26 70 
31 | 5.471 § 40 © 2215 13 £0713 55 xara) 2.36 35| I 25 42 
305 15.7515 47 * 4] 5 12 & 25] 3 54 XI151| 2 35 XZ0I| I 24 4395 
20 | 6.05 S511 * 31353 * 32] 2 34 # 87 1 22 2525 
205 | 6.37 13 50 42.31 232 % oO| 1 20 % 14 
0.292 | 6 57| 1348 «O:5] 2 30 “3.5| 3.28 £3.06 
| 
| T = 2375 ( T=24°8¢ T=26°0 C | T=27°4C | T=27°9 C 
Range in | 
O60'S5....43 ..| 5.76-3.86 3.52-2.64 2.47-1 . 58 t.32-1.22 | 1:02-5:08 


| 
| 
| 





* The symbol x indicates sector was stopped, increasing deflections twenty fold. 


t Additional temperature observations: a 25°3, b 25°5, c = 28°5, d = 28°98, and e = 27°6. In 
set 5 the measures to the violet of \ 0.325 w preceded the others 

data for five sets of measures, consisting of the hour angle followed 
by the observed galvanometer deflection in millimeters. The nota- 
tion x before a deflection means that the 5 per cent sector was 
stopped with its aperture in front of the slit, thus multiplying the 
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deflections by 20. At the top of each column is noted the second-slit 
zero and below it the temperature. The temperature is also noted 
where a discontinuity occurs in the hour-angle column and at the 
end of each set. Below is given the range in sec z (air mass) for that 
column. The negative sign of the hour angle has been omitted, as it 
applies throughout the table. 
TEMPERATURE EFFECTS 

There are several possible effects of temperature changes on the 
thermocouple and the galvanometer circuit. The resistance of most 
metals, particularly those used in making the thermopile and the 
copper of the galvanometer circuit, increases about 0.4 per cent per 
degree centigrade. On this account alone we would expect the de- 
flections to diminish by a like amount with increasing temperature. 
There is also the slightly decreasing torsion of the galvanometer 
suspension, which would increase the deflections, and the slight 
change in thermoelectric power of the thermopile and its decreasing 
heat engine efficiency, 7~'AT. Against these considerations must be 
placed the uncertain effect of temperature on the vacuum in the 
thermopile cell with its calcium tube. This is a very important item, 
for, if the vacuum is near the limit corresponding to the steep siope 
of the pressure-sensitivity curve,’ the deflections may be consider- 
ably decreased if the temperature change slightly raises the pressure. 

To test the question, the thermopile (freshly exhausted) and the 
galvanometer used in the present measurements were inclosed in a 
paper box, provided with windows to admit light to the thermopile 
(one receiver of which was shaded) and to expose the galvanometer 
mirror. The temperature was controlled by a 75-watt lamp within 
the box. The deflections were obtained by switching on a 1oo-watt 
lamp at 3 meters’ distance. This lamp, with a voltmeter across the 
circuit, was fed from a D.C. generator. The voltage was maintained 
at 115 by means of a microvoltage control. The readings made 
are given in the following table. For the first to of tempera- 
ture increase the deflection decreased only 1 per cent, and as this 
interval includes the range usually occurring during observations 
of the solar energy-curve, we may conclude that the effect of ordi- 


4E. Pettit and S. B. Nicholson, Mt. W. Contr., No. 246, p. 7; Ap. J., 56, 301, 1922. 
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nary temperature changes on radiation sensitivity may be neg- 








lected. 

io hl . ° | ° | ° ° | ° 
Temperature (C) ...... 33:8: | 27.0 | 3054, ) “agus 39.0 
Deflection (mm)....... | 330 338 | 338 | 335 | 320 





The slit mechanisms were such that increasing temperature de- 
creased the slit widths. While the effect on the first slit differed from 
that on the second and the third, we may for convenience suppose 
that the effect on all three is involved in the second slit alone. Let us 
suppose that the normal width of the slit is 0.043 mm for the first 
morning observation at \ 0.7 w and the corresponding deflection D,; 
succeeding settings give D,..... Since the atmospheric transmis- 
sion on Mount Wilson for \ 0.7 u is 0.95 (Abbot), the deflection D 


for any air mass is 
log D = (log D, — 9.978 sec 2) + 9.978 sec 2. (1) 


If the fundamental atmospheric transparency is constant and if 
the sensibility of the galvanometer circuit remains unchanged or can 
be corrected, we may attribute any deviation of the observed de- 
flections from those calculated by equation (1) to change of slit 
width. Suppose the observed deflection is D’, then the change in 
slit width for any sec 3’ is 

bw = 0.043 (D’D™ — 1). (2) 

In practice the values of dw were plotted against hour angle and 
the corrections to slit width were taken from the curve. These cor- 
rections apply to all wave lengths. The correct deflection D, for any 
wave length at any air mass is then 

D. = Daw(w + bw), (3) 
where w is the tabular slit width and D, the observed deflection. The 
value of é6w did not usually exceed +0.008 mm. 

PROPORTIONALITY OF DEFLECTIONS TO SLIT WIDTH 


Since the second slit is opened in proportion to the dispersion for 
each wave-length setting, it is important to know if the deflections 
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are proportional to this opening. This was investigated for \ 0.32 u 
with both thermopile and photocell, and for \ 0.6 » with the thermo- 
pile. Deflections were read for slit openings increasing by intervals 
of 0.05 of a screw revolution over the range used in observation. 
When plotted against slit width, these deflections form a straight line 
which does not quite pass through the origin, showing that the pro- 
portionality rule holds save for slits narrower than 0.05 revolutions. 
In these data the reading for zero slit width was carefully deter- 
mined; beginning at 0.05 revolutions, the line representing the plot 
is therefore curved as it approaches the origin. To obtain the correc- 
tions a line is drawn parallel to the plot through the origin; the 
ratio of any ordinate on this line to that on the plot will be the cor- 
rection to an observed deflection to make it proportional to slit 
opening. These corrections appear in Table 2, column S. 


REDUCTION OF THE MEASUREMENTS 

The logarithms of the observed deflections (log d), corrected if 
necessary by equation (3), were plotted against air mass (sec 3), 
and the straight lines best representing the observations were ex- 
tended to the zenith (air mass 1) and to no atmosphere (air mass 0). 
Figure 2 shows such a plot for a satisfactory day’s run. From these 
plots the values of log d, and log d, were read for air masses 0 and 1, 
giving d,, the deflection for no atmosphere, and log t, = log d,; —log do, 
the logarithm of atmospheric transmission. Inspection of Fig- 
ure 2 shows that the values of d, between 0.38 and 0.325 yu 
cluster about the same ordinate on air mass o. Since the corrections 
to be applied to d, are not large, it is evident that the ultraviolet 
solar intensity does not vary greatly over this range in wave length. 

The corrections involved depend somewhat on the arrangement of 
the optical system. There is always the reflecting power of the sider- 
ostat mirror. During part of 1934 there was one, and throughout 
1937 a second, additional mirror. For the remainder of 1934 a fused 
quartz lens replaced the added mirror; and this lens was also used 
in 1939. Then there are two quartz monochromators and, except in 
1939, a fused quartz transfer lens. The window on the vacuum 
thermopile is so thin (1.5 mm) that the losses are sensibly those by 
reflection from the two surfaces. According to Fresnel’s formula, the 
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Fic. 2.—Plot of logarithm of deflections against air mass, for June 8, 1939. Note 


that between d 0.38 and 0.325 uw the deflections at air mass o have approximately 


the same value. 
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loss for one surface at \0.3 uw is 0.0384 and at X0.7 w, 0.0344—a 
range of less than o.5 per cent, which can be neglected. When the 
photoelectric cell was used, a fused quartz prism was thrown in front 
of the transfer lens to send the light to the cell. Since the optical 
paths through the prism and lens are about the same, this alteration 
of conditions need not be considered. 

If Rk. be the reflecting power of an aluminized surface,’ Tp the 
transmission of fused quartz’ (both the image-forming and the 
transfer lenses were of approximately the same thickness), T, the 
transmission of one monochromator,’ F the focal length (or aperture 
ratio) of the fused quartz image-forming lens, and S the reciprocal 
of the slit efficiency, we have the following correction coefficients K 
to be applied to the mean value of d, for each wave length: 


1934t K,=SR,*T,;? TR’, (4) 
1934 A K,=S F? Rz' Tz? Tp’, (5) 
1934 BK, = S Rx? Tz? Tz, (6) 
1937 K, = SRz3T;z? Tz, (7) 
1939 K, = SF? Rz' Tz? TR. (8) 


Formula (4) refers to integrated light, (5) to a lens projector, (6) 
to a single-mirror projector, (7) to a two-mirror projector; (8) is the 
same as (5) with the transfer lens and third slit eliminated. The 
numerical values of A are given in Table 2. It will be noted that 
the ranges in AK, and A, are much less than for the other coefficients. 
This is due to the compensating effect of the change in focal length 
F of the quartz projecting lens with wave length. In fact, if a tele- 
photo arrangement were used, a nearly complete compensation 
would be possible. 

After the averages were taken and before AK was applied, the 
measurements were reduced to a scale of watts per square meter per 
100 A by making the values at \ 0.45 u equal to 22.4 for integrated 
light and to 26.0 for the center of the disk. The method of arriving 
at the former value is as follows: The energy in the integrated 


Pettit, Pub. A.S.P., 46, 27, 1934. 
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solar spectrum® between 0.45 and \o.70 uw is 0.71 cal cm~? min, 
or 4.95 X 10° erg cm~? sec~*. If we plot the solar energy-curve for 
no atmosphere as given in arbitrary units in the Smithsonian 
Physical Tables’ over this spectral region, we find that the area 
within an interval of 100 A at \ 0.45 u is 0.0453 of the whole. The 
corresponding energy is 0.0453 X 4.95 X 10°, or 22.4 X 103 erg 


cm~? sec~' (100 A)~", or 22.4 watts per square meter per 100 A. In 
the previous calibration used for the work at Tucson, the approxi- 
mate value 20 watts was used; but it is believed that the value here 
found will give a more accurate scale, and it has therefore been 
adopted. A calibration from Planck’s radiation formula for a solar 
temperature of 6000° K, reduced to mean solar distance, yields 21.4 
watts per square meter. 

To obtain the calibration for the center of the disk the area of the 
solar energy-curve was increased by the factor 1.16, the ratio’ of 
the average intensity of total radiation for the whole solar disk to 
the intensity at the center. The slight shift of the maximum to the 
violet, which the Smithsonian observations seem to show in the 
energy spectrum for the center of the disk, has been neglected; it is 
too small to be of importance here. 

The results in watts per square meter per 100 A are shown in 
Table 2, where 7 following the date indicates integrated light, and 
c light from the center of the disk. ‘‘Av.” is the average value, in- 
cluding the 1931 work at Tucson reduced to the new scale. “‘Tra.” 
is the average atmospheric transmission for all the observations on 
Mount Wilson, a total of 35 determinations for each wave length. 

Figure 3 shows a plot of the results for each year’s work. In com- 
paring them it should be remembered that they have unequal weight, 
the number of days of observations being different, as shown in 
the following table. But since all these different series of measures 
apply to different years, no weighting of normal places was made in 
taking the averages in Table 2 or in drawing the mean curves of 
Figure 3. The same general features are seen in this plot as were 
found in the 1931 plot, viz., a sudden drop of 48 per cent in the solar 

© Smithsonian Physical Tables (8th rev. ed.), p. 608, 1933. 

E. Pettit and S. B. Nicholson, Mt. W. Contr., No. 397; Ap. J., 71, 153, 1930. 
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energy between \ 0.40 and 0.38 uw and a nearly constant value of 
the solar energy from \ 0.38 to 40.325 wy. 


| 1931 1934 | 19037 1939 
Center ....... : el 2 | 10 | 9 10 
Integrated 4 7 

| | 





The ordinates of the curve for integrated light are nearly every- 
where proportional to those for the center of the disk. This is par- 
ticularly interesting in the ultraviolet, where the rise in intensity 
between 0.38 and 0.35 yu is revealed in both curves. Below 
Xo.32 u, the integrated-light curve runs even a little higher than 
that for the center of the disk. This is hard to understand, since it 
is supposed that darkening at the limb would tend to make it fall 
even farther below than for higher wave lengths. The same phe- 
nomenon was noted’ in 1931. There is some indication here that 
darkening at the limb does not vary greatly throughout the ultra- 
violet region. The high values for integrated light found in 1931 for 
AA 0.30, 0.295, and 0.292 w were not verified in 1934 and may be 
erroneous. The curve through the averages for these wave lengths 
has been shown here as a dotted line, the full line including only the 
Mount Wilson values. The fall in intensity from \ 0.325 to 0.292 
is practically linear. This interval is all within the ozone band, and, 
while the effect of the band should be eliminated by the method of 
reduction, the steepness of the transmission plots (Fig. 2) leads to 
considerable uncertainty. This is all the more patent when we con- 
sider the scattering of points below \ 0.30 uw in Figure 3. 

There is, however, an effect which would make the measures be- 
low \. 0.30 uw too high. The second slit opening, being always too A, 
includes the region AX 0.30-0.29 uw for the setting \ 0.295 uw, and 
AA 0.297—-0.287 w for the setting \ 0.292 w, and the thermopile inte- 
grates the radiation over these wave lengths. If the decrease of 
radiation across the 100-A slit decelerates toward the violet or if it 
falls to zero within the second slit, the observed deflections will be 
too high or will refer to a higher wave length than that correspond- 
ing to the center of the slit. For wave lengths greater than 0.30 u 
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this criticism does not apply. This effect would also tend to make the 
atmospheric transmission coefficients too large. 


ATMOSPHERIC TRANSMISSION AND OZONE 


The atmospheric transmission coefficients have been found for both 
integrated light and light from the center of the disk, the data for 


TABLE 3 


ATMOSPHERIC TRANSMISSION ON MOUNT WILSON 


| | 





Au 19034 1037 1939 Av. | Smithsonian 

0.7 | 0.950 0.938 0.941 | 0.943 | 0.950 

6 all Qos 893 943 | oI! . goo 

55 | 884 875 888 } Gee Eesscevuswes 

cO 869 866 871 | 869 862 

As 827 840 827 | 832 ( .825)* 

42 790 784 793 | TE lea bea wero ae 

40 794 740 743 752 729 

39 735 740 731 742 win are 

28 716 60909 700 705 | 676 

37 682 OSS 670 | 6090 © wie ae oe 

30 6062 640 640 050 | 6035 

35 O17 007 592 005 Pr 

34 585 65 45 | 565 | 580 

22 541 520 501 52 } 5 hacbanes Sooke 

325 502 195 470 49I hed 

32 $70 462 431 450 } 520 

315 421 403 350 401 

31 349 325 | 311 329 Er or 

305 272 238 227 | QQO foseecasives 

30 195 167 137 166 | (0.460) 

205 114 o80 o71 | >| aE Poe eee ee 
0.292 0.078 0.026 0.037 0.047T 


* Interpolated 
t When corrected for slit width, these become 0.071 and 32, respectively. 


which may be combined by taking simple means for any one season. 
These coeflicients were obtained for 15 days in 1934, 9 days in 1937, 
and 11 days in 1939. The averages for each year, together with the 
corresponding coefiicients from the Smithsonian work on Mount 
Wilson,° are shown in Table 3. 

We may determine the atmospheric ozone in the same manner as 
was done for the Tucson observations. From Rayleigh’s law of molec- 


ular scattering, 


log t, = —3273 (3N)“H log e (u — 1)?A4, (9) 
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LOG 
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we see that for zenith conditions, where the atmospheric depth H 
and the molecular density N are constant, the logarithm of the 
transmission, log ¢,, depends on the factor — (u — 1)*A~4. Hence if we 
plot the observed log ¢, against a wave-length scale which expands 
according to the function —(u — 1)?A~4, we should expect a straight 
line as long as the atmospheric interference is of pure molecular 
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Fic. 4.—Atmospheric transmissions plotted on a wave-length scale that expands 


according to the function —(u—1)?\—4. Heavy lines indicate ozone transmission for 


layers 0.08 and 0.10 cm equivalent thickness at N.T.P. 
é, | 


Figure 4 is such a plot of the observations for 1934, 1937, and 
1939 from the data in Table 3. It is seen that the portion to the red 
of 0.33 uw is represented by a straight line, but for shorter wave 
lengths the deviation increases continuously. If we extend the 
straight portion of the line to the violet, the difference in ordinates 
for any wave length, in the sense of extended line minus observed 
curve, is the logarithm of the absorption due to atmospheric ozone. 

The coefiicient of absorption of ozone for a layer 1 cm thick at 


N.T.P. given by Fabry and Buisson is 


log a = 17.58 — 0.00564), (10) 
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where \ is in angstroms. The thickness of the equivalent ozone 
layer due to any observed atmospheric transmission fp is 


h = a" (log t, — log tp) , (11) 


where / is in centimeters and log ¢, is taken from the projected 
line for the same wave length. 

The more direct way of estimating ozone is to compute (log #, 
— log tg) from equation (11) for assumed values of ¢ and compare 
the curves with the plotted points. This has been done for h = 0.08 
and o.10 cm, whence it is seen that the Mount Wilson observations 
indicate an atmospheric ozone content of about o.1 cm. The 1937 
series gives 0.08, and the 1939 series about 0.11 cm. These values 
are considerably lower than that found at Tucson (0.18 cm) at the 
same season of the year. 

It will be noted that the values of tp for \ 0.295 and \ 0.292 w are 
greater than those for the computed ozone-curves. The values in 
note t of Table 3 will bring these points between the two ozone- 
curves. This is again evidence that the measures in these wave 
lengths should be made with narrower slits. 

On account of the low value of the atmospheric ozone content 
above Mount Wilson, it was thought that some error might have 
occurred in the calibration of the wave-length dials in 1934, which 
was done with a mercury arc before and after the series. This cali- 
bration was checked in the following winter, after the observations 
had been reduced, and no error was discovered. In 1937 the wave 
lengths were checked in the middle of the series, as well as before 
and after; and in 1939 the checks were made before and after each 
day’s run. It seems then that the observed low value of the ozone 
content of the atmosphere above Mount Wilson cannot have its 


origin in wave-length error. 
ULTRAVIOLET ENERGY BETWEEN THE LINES 
The energy-curves in Figure 3 represent the total solar energy 
received from each interval of 100 A as modified by the absorption 
of the Fraunhofer lines. If these lines could be eliminated by measur- 
ing the radiation between them and if the wings of the lines did not 
seriously overlap, the result would be the energy-curve of the con- 
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tinuous spectrum of the sun. In the ultraviolet, however, there is 
no space between the lines in which the overlap of the wings can be 
ignored, and any direct measurement must give something less than 
the energy of the continuous spectrum. 

To investigate the matter the 21-foot concave-grating mono- 
chromator, arranged in Eagle mounting, was provided with an 8-inch 
grating giving a bright first-order spectrum. The turntable sup- 
porting the grating is mechanically connected with the gear train 
that drives the recording drum of the galvanometer. Light entering 
the first slit and falling upon the grating returns as a spectrum along 
approximately the same path and passes through a second slit be- 
hind which is placed a vacuum chamber containing the quartz 
photoelectric cell and amplifying (“‘pliotron’’) tube. A fused quartz 
window admits light to the vacuum chamber. The control box and 
batteries are placed near the galvanometer (Leeds and Northrup 
goo ohms). A 2-inch concave mirror of 26 feet radius is mounted 
on the grating vertical axis. This mirror casts an image of a straight- 
filament lamp upon a scale graduated in 1oo-A units for setting pur- 
poses. A fused quartz convex cylindrical lens can be put behind the 
first slit to bring the focal plane to the slit jaws for any wave length.® 
A filter holder is inserted in the light path near this lens. 

A concave-mirror projecting system of 4o feet focal length was 
used to form the solar image upon the first slit. Ordinarily the two 
slits had the same width, about 0.08 mm. For registration to the 
violet of 4.0.40 w a Shott UGr filter was used to reduce scattered 
light. Because of the sluggishness of the electrical system and the 
3-second period of the galvanometer, about 15 minutes were required 
to record each r1oo-A interval. For this reason records were not 
made when sec z was greater than 2. 

For the region \\ 0.49-0.37 uw the gas-filled potassium cell was 


The color-sensitivity function of the sodium cell has already been 
given.’ Tests made with a photographic step wedge showed that the 
galvanometer deflections were proportional to the radiation in the 
21-foot monochromator when used with these cells. 


SE. Pettit, Pub. A.S.P., 43, 75, 1931. 
9 Mt. W. Contr., No. 445, p. 25; Ap. J., 75, 209, Table VII, 1932 
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Usually the recordings were made over the integral hundreds of 
angstroms, for example, AX 0.35-0.36 u, AX 0.36-0.37 w, etc. Con- 
tour lines were drawn through the apparently highest points of the 
records and subdivided into 50-A units (Pl. I). The areas below 
the contour line and below the curve itself were then measured with 
a planimeter. The ratio of the former to the latter, averaged for two 
contiguous regions, is the ratio of the energy in the 1oo-A interval, 
as it would be if measured between the lines most widely separated, 
to that measured with the thermopile and quartz monochromator 
shown in Figure 3. The data for the 50-A units were then combined 
by simple averages to give the ratios for AXA 0.355-0.365 wu, etc., 
which correspond to settings of the quartz monochromator at 
0.360 uw, etc. This procedure was necessary since the planimeter 
was not large enough to measure an entire 100-A record. 

Plate I shows galvanometer records for the familiar H and K 
region and for AX 0.32-0.33 uw. In the former the contour has been 
drawn over the H and K lines for comparison with the contours 
found by A. D. Thackeray,'® who used the method of photographic 
photometry with a 21-foot concave grating. Thackeray’s observa- 
tional points for his contour of the K line are shown by the heavy 
dots in the figure. Since he gives no data for the H line, some of his 
observational points for K have been transposed upon the H con- 
tour to show the agreement between the two contours; they should 
not be mistaken for observations of the H line itself. This compari- 
son shows that the agreement for K is as close as could be expected 
from two such different methods of procedure and that the con- 
tours of H and K are nearly alike. The region AX 0.32-0.33 pw is at 
the beginning of the ozone band, and the positions of several of the 
broad ozone lines are indicated on the curve by the braces under the 
symbol O,. 

Table 4 shows the ratio, e, of the energy between the lines to that 
in the integrated spectrum averaged over too A, found by the fore- 
going method. The column eM is obtained from estimates made by 
G. F. W. Mulders™ based on a calibration-curve of Rowland’s 
scale of intensities. The most serious disagreement is for \ 0.39 
and \ 0.40 yw; but this at least cannot be due to scattered light in the 


M.N., 95, 293, 1935. 11 Zs. f. Ap., 11, 132, Table 2, 1935. 








180 EDISON PETTIT 


21-foot monochromator, since it is near the red cutoff of the filter 
and the contour of the K line agrees with that of Thackeray, as 
shown in Plate I. 

The quantities e Av. i and e Av. c are plotted in Figure 3 as 
broken lines. It is seen that the effect of correcting the thermo- 
electric measures for line absorption in this manner is to smooth out 
the violet slope of the energy-curve of the sun. Although the major 
features are still present, the energy in the ultraviolet no longer ap- 
proximates a constant value. The black-body radiation-curves for 


TABLE 4 


RATIO e OF THE ENERGY BETWEEN THE LINES TO 
THAT IN THE INTEGRATED SPECTRUM 


Au | e eM @Av.t | ¢€Av.¢ 

6.320. . I.39 E.22 12.8 | 12.6 
330. i .33 E32 r2.0 14.3 
34° 1.43 I. 33 14.1 15.7 
350 1.35 | I Ra 3 '$..2 
ac. I.47 1.54 10.1 15.1 
37 1.43 | 1.61 16.0 17.8 
380 1.50 1.59 10.1 15.5 
390 1.50 1.55 20.9 24.1 
305 1.61 ee 29.0 
400 i .32 1.54 a3.2 27.3 
420 cas 23..7 26.0 
445 I.14 I.r2 25.5 290.6 
485 I.12 Pee 24.5 28.6 
500 L316 24.4 28.1 
550. 1.09 22:7 25.8 
600 I.02 20.2 42-9 
0.700 I.O1 £3.77 15.6 


6000°, 6200°, and 6400° K have been drawn for comparison. These 
were obtained by use of Planck’s radiation law, reduced to mean 
solar distance. The application of the coefficient e to the measures 
most nearly approximates 6200° for integrated light and 6400° for 
the center of the disk. In the ultraviolet, however, the corrected 
energy-curve of the sun still falls far below that of a black body. 
Since the thermoelectric measures at \ 0.39 uw include the K line 
and those at 0.40 uw the H line, the effect of the absorption of the 
two together is smoothed out. To estimate the intensity at \ 0.395 y, 
we take the quotient of the ordinates of the smoothed curves of 
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e Av.i and e Av.c at that wave length and e for \ 0.395 wu. These, 
plotted upon the original energy-curves, give the depressions at 
0.395 w seen in Figure 3. It is significant that the H and K lines 
seem to play little part in the great drop of intensity between 
dh o.40 and \.0.38 uw. There are no lines in the ultraviolet solar spec- 
trum attributable to terrestrial origin until we reach the Huggins 
band of ozone superposed upon the great absorption band at the 
violet end of the solar spectrum. The Huggins band consists prin- 
cipally of a cluster of pairs of broad lines shading to the red, spaced 
about 28 A apart, with the individual lines of each pair separated 
about 7 A. The width of the lines themselves is about 8 A, and of a 
pair about 13 A, the shading of the violet component overlapping 
the red. Since their intensities are rather low, the lines are entirely 
masked by lines of solar origin, even with low dispersions. 

In 1926 the author made a laboratory study of the ozone spec- 
trum, using the first order of a 1-m concave-grating spectrograph, 
mounted Eagle fashion (dispersion 17.4 A/mm). The ozone tube, 
67 cm long, closed by crystal quartz windows, was charged from a 
Berthelot generator operated from a 50,000-volt transformer. With 
this equipment an ozone density could be obtained equivalent to a 
layer 1.5 cm thick at N.T.P. The analysis of the ozone content of 
the tube was made by the iodine method. The exposures were on 
Cramer Contrast plates, with light from the crater of the carbon arc. 
An iron-arc comparison spectrum was added to each exposure by 
using a focal-plane occulting bar. Several exposures were impressed 
on each plate together with exposures through a calibration wedge 
for intensities. 

From this material four spectra referring to a variety of ozone 
densities were selected as having the best exposures. These were 
measured direct and reverse and referred to the comparison lines in 
the usual way. The results are given in Table 5. The first column 
gives the wave lengths, those in parentheses being measured by 
Fowler and Strutt’? only. Three lines not given by Fowler and 
Strutt were found at AXA 3148.6, 3213.9, and 3440.2. Since the ac- 
curacy of measurement cannot be very great, only the nearest 0.1 A 
is recorded. The second column gives the mean residuals for the 


12 Proc. R. Soc., 93, 577, 1917. 
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four spectra, and the third the differences between the present meas- 
ures and those of Fowler and Strutt. The strong lines bracketed in 
the first column are members of pairs, and their positions are 
marked in Plate I by their intensity numbers. The braces in this 
plate show the region of ozone absorption for each pair of lines. It 
appears that the Huggins band only partly accounts for the general 


TABLE 5 
WAVE LENGTHS AND LINE INTENSITIES IN THE 
ULTRAVIOLET OZONE BAND 





| | a To ee 
Mean | Mt. W. | | | Mean | Mt. W 





| 
| | | 
rN | Residu- | minus | I d | Residu- | minus | I 

| als F and S als |FandS| 

| | 
(3089.5) | 8 3242.0 | O88 | —Xe0:) a 
(3096.5) | 4 3248.5....| I Lie es 
(3105.0) L -¥ 3255.0 2 | —0.5 6 
3114.7 +o.4| 8 | {3271.9 1 | —o.1 2\( 
> 7” — > | ie) 
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absorption at AXA 3193, 3213, 3227, 3250, and 3281 and that these 
depressions are mostly due to groups of the stronger solar lines. 
The intensities in the fourth column are on a scale of to. For 
wave lengths shorter than \ 3187 they are taken from Fowler and 
Strutt (save \ 3148.6, which is not in their list), who used greater 
dilutions of ozone. Values in parentheses, each referring to a pair of 
lines, were obtained by microphotometry, the calibration wedge 
being used to convert the photometric curves into intensities. The 
intensity ratio with ozone in the tube to that with oxygen was 
formed for each wave length. Plotted against wave length, they 
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formed a curve with depressions for each ozone line. The ratio of the 
depth of each depression to the ordinate of the smoothed curve is 
then the absorption of the line. 

The sums of the estimated intensities for a pair of lines approxi- 
mate the measured intensities for 1 cm of ozone, if we disregard 
the decimal. Since the maximum in absorption is 18 per cent and 
averages 10, we would expect only about 1 per cent absorption in 
the bands in the region Ad 0.32-0.33 uw indicated by the record in 
Plate I, which was made at air mass 1.3 on Mount Wilson. From 
this discussion we see that the Huggins band of lines will be scarcely 
distinguishable in the clusters of Fraunhofer lines in a G-type spec- 
trum. 

PRELIMINARY LINE-GROUP INTENSITIES 

The individual lines are not sufficiently separated in the photo- 
electric energy-curves from the 21-foot monochromator (PI. I) 
to permit accurate measures of their intensities. Such curves, how- 
ever, furnish material for a study of the intensities of groups of lines 
included between wave lengths at which the continuous spectrum 
approaches the general contour of the energy-curve. The 100-A 
strips for AA 3900-4000 and AA 3200-3300 A were accordingly sub- 
divided as shown in Plate I, the adopted wave-length limits being 
those given in Table 6. 

Mulders'} has defined the absorption of a line in terms of the 
width of an equivalent totally absorbing rectangular line in thou- 
sandths of an angstrom (EW in mA). We may determine EW, the 
line-group absorptions (intensities) according to this definition, by 


the formula 


EW = 1000 A\ A; AZ", (12) 


where Ad is in angstroms, A; is the area of the line and A, the area 
below the contour, both within the wave-length limits assigned. The 
results for EW from planimeter measurements are given in Table 6. 

The corresponding line absorptions given under ‘“‘Mulders’ cali- 
bration”? were obtained by summing for the lines within each inter- 
val the values obtained from Mulders’ calibration-curves of the 


13Z5.f. Ap., 10, 297, 1935. 
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Rowland intensities, which overlap the region AX 3900-4000 A. 
Mulders’ measured intensities were, however, used whenever given 
in his table. 

TABLE 6 





| | 


| 


Mul- | Mul 
Spectral . ders P Spectral : ders’ P 
Limits ; a) Limits = 7 oo) OS 
| bra- | bra 
tion | tion 
mA mA | mA mA | 
3894 .5-97.0 775 656 {1.39 || 39609.0-70.8 54 802 |0.07* 
97.0-98.8 390 501 |O.79 10 771 7 30 223 |O.10 
98 .8-o1 .2 816 908 |o.go* 71.7-72.7 210 216 10.97 
3901. 2-04.4 1056 | 1206 |o.88* 72.97-75.4 459 826 |o.56* 
04.4-07.4 1380 | 1520 |o.g1* 75.4-77.2 | 360 466 |0.77 
07 .4-09. 5 | 504] 533 |0.95 77.2-79.0 | 288 | 382 jo. 76* 
09Q.5-11.0 525 044 |O.5I 79.0-50.0 100 227 |0.44 
[lt .O-12.9 195 294 |0.66 50 .0-55.5 1045 1244 |0.54 
I2.Q-15.2. 460 687 |0.67 85.5-88.0 550 570 |0.96* 
15.2-18.2 690 832 |0.83* 88 .o-9g0.8 560 546 |1.03 
18.c—20.0 513 | 615 |o.84 go .8-93 .6 392 505 |0.78 
20 .0-22.5 1303 968 {1.41 3093 .0-00 1564 1516 |1.03 
22.5-24.0 494 | 459 |1.07 || 3189.5-98.9 4251 
24.0-25.0 |} I50 | 205 |0.70 Q5.9-04.5 1410 
25.0-260.7 | 204] 522 |0.56 3204. 8-06. 5 255 
26.7-29.0 | 630] 582 |1.08* 06.5-19.0 3750 
29 .0-32.0 | 506 581 |0.87* IQ.0-21.0 320 
K  15.0-53.0 11820 |15900 |0.75 21.0-24.5 980 
35.8-38.0. 28 172 |0.16 24..5§-32.5 2880 | 
38 .0-40.0 216 292 |0.74 $2:.5-35:5 2160 
40.0-40.8 2 68 |o.35 38.5-40.7 550 
40.8-42.8 360 | 629 |0.57 40.7-46.5 925 
42.8-46.0 102 1189 |0.86* 46.5-55.5 2790 
49 .0-49.5 600 955 |0.63 55.5~-50:5 440 
49.5-51.0 126 224 |0.56 59.5-62.2 297 
51.0-54.9 | 720 978 |0.74* 62.2-65.09 666 
54.9-57-5 672 | 810 |o.83* 65.9-75.0 924 
57 «5-50: 255 315 |0.92 75 -0-59.9 3954 
59.1-64.9 1500 | 1617 |0.93* 89.9-93.4 s16 
04.9-07.0 1003 450 |2.09 93.4797 .5 405 
H 3953.0-82.0 10050 3297.5-00.0 510 
* Region includes lines of Rowland intensity greater than 5 or more than one of intensity 5 


The column P/M gives the ratios for the two series of group 





intensities. The average ratio is 0.84, but the scatter is very large. 
If we examine the 15 ratios (indicated by asterisks) for which the 
spectral intervals include lines of Rowland intensities greater than 5, 
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or more than one line of intensity 5, the scatter is much reduced 
save for the interval \X 3969.0-3970.8 A near the H line, where the 
photoelectric curves give no indication of intensities as high as 10 
and 5 assigned by the Rowland table. The scatter is no doubt 
indigenous in the Rowland intensities themselves in this difficult 
region dominated by the H and K lines where intensities of other 
lines must be measured from the contours of the wings of H and K. 

The intensity of the K line is 11,820 mA while H is 10,050 mA, 
a ratio of 1.18 against the theoretical ratio 1/2 = 1.41 adopted by 
Mulders to obtain the value of H from that of K. 

CARNEGIE INSTITUTION OF WASHINGTON 

MounT WILSON OBSERVATORY 
October 1939 








GALACTIC STAR CLUSTERS' 
ROBERT J. TRUMPLER 


ABSTRACT 

The principal methods for determining intragalactic distances greater than 500 par- 
secs are discussed and their application to galactic star clusters is outlined. A compari- 
son of distance moduli derived from photometric spectroscopic data with distances 
based on angular diameters yields an average coefficient of photographic absorption 
of o™7 per kiloparsec. Preliminary results for the radial velocities of 39 clusters lead 
to a value A = +0.015 km/sec - parsec for the coefficient of differential galactic ro- 
tation. The good agreement of this result with other determinations confirms the cor- 
rectness of the adopted distance scale of galactic star clusters. The radial velocities 
also indicate a negative & term increasing in amount with distance. 

The appearance of galactic star clusters seen at great distances projected on the 
background of faint stars is studied. At distances of more than 5000 parsecs even the 
most favorable objects would have been missed by the surveys of the past, while the 
poorer and less luminous objects may have been overlooked at much smaller distances. 
This explains why the space distribution of the known star clusters shows little relation 
to the structure of our galaxy. Attention is drawn to the promise of a search for more 
distant clusters on photographs taken in red light with instruments of large aperture 


and large field. 
I. INTRODUCTION 
Among the many problems with which the study of galactic star 
clusters is concerned it seems appropriate here to take up those which 
bear directly on the structure and dimensions of our galactic star 
system. On the one hand, we have the distance determination which 
involves the study of interstellar light-absorption and galactic rota- 
tion, and, on the other hand, the space distribution of known clus- 
ters and the question of incompleteness in our present knowledge of 
these objects. 
For the determination of galactic distances greater than 500 
parsecs the following methods are available: 
I. Photometric methods 
a) Luminosity from spectral class 
b) Luminosity from periods of variable stars 
(Requires a knowledge of the relation between absolute magnitude 
and spectral class, or the period-luminosity relation, and of the absorp- 
tion of light in interstellar space.) 
t This paper was presented at the Symposium on Galactic and Extragalactic Struc- 
ture, held in connection with the dedication of the McDonald Observatory on May 5-8, 


1939. 
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II. Angular diameters of extended objects 
a) Galactic star clusters 
b) Planetary nebulae 
(Requires a knowledge of the linear diameters of various classes of such 
objects.) 
III. Differential galactic rotation effect on radial velocities 
(Requires a knowledge of the Oort constant 4.) 

Strictly speaking, all these methods are of a statistical nature; the 
first two, however, may within reasonable error limits be applied to 
individual objects, while the third can be used only to find the mean 
distance for groups of galactic objects. Moreover, these methods are 
not absolute; they require for calibration a knowledge of certain con- 
stants or statistical relations which have to be established empiri- 
cally. Some of these, like the luminosity-spectral class relation or 
the zero point of the period-luminosity relation, result from a study 
of the nearer stars for which trigonometric and proper-motion paral- 
laxes are available. Others, such as the average coefficient of inter- 
stellar light absorption, the linear diameters of clusters and nebulae, 
or the galactic rotation constant cannot be reliably determined from 
the nearer objects and have to be studied by intercomparing the re- 
sults of different methods of distance determination. The distances 
of the more remote galactic objects depend considerably on the val- 
ues adopted for the absorption coefficient and for the galactic rota- 
tion constant A. Star clusters offer the possibility of applying all 
three methods—a particularly favorable opportunity for the calibra- 


tion of the distance scale. 


ll. THE PHOTOMETRIC-SPECTROSCOPIC DISTANCE DETER- 
MINATION OF STAR CLUSTERS 

The procedure for the photometric-spectroscopic distance deter- 
mination of a galactic star cluster is as follows: 

1. The apparent photographic magnitudes m and the spectral 
classes s of the cluster stars are found by direct observation. 

2. After the giant stars have been separated from those of the 
main sequence and after the background stars have been eliminated 
as far as possible, the mean apparent magnitude m, of the cluster 
stars is formed for each spectral class. This is compared with the 
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mean absolute magnitude M/, of the same spectral division in the 
Hertzsprung-Russell diagram of the nearer stars. 
3. The distance modulus y of the cluster is given by 


y=m,— M,, 
and its mean value for all cluster stars is adopted. 
4. The relation between the distance modulus y and the distance r 
(in parsecs) is 
y= s5logr—5+A(r), (1) 


where A(r) is the amount of interstellar absorption affecting the 
light of a star at distance r. 

The determination of the distance modulus involves the assump- 
tion that the luminosity-spectral class relation in a cluster is the same 
as for stars in general. Kuiper’ suggests that the chemical compos'- 
tion (hydrogen content) of different members in the same cluster is 
very similar, while it may vary considerably from one cluster to an- 
other. A difference in hydrogen content, however, affects both the 
spectral class and the luminosity in such a way that a star of the 
main sequence is principally displaced parallel to the sequence. Al- 
though our assumption may not be strictly fulfilled for individual 
clusters, the error introduced will, in general, be small. 

The absorption law A (7), which must be known in order to find the 
distance r from the modulus y, does in reality vary with galactic 
longitude and latitude. Our information on the irregularities of the 
interstellar medium is as yet too incomplete to allow for these. Two 
simplifying assumptions about the absorption may be introduced as 
an approximation: 

1. For the true absorption correction we substitute its statistical 
average for low galactic latitudes, and this average absorption is as- 
sumed to be uniform, 

A(r)=a-r, (2a) 


where a is a constant representing the average absorption of photo- 
graphic light per kiloparsec near the galactic plane. 

2. The photographic absorption A (7) may be taken as proportion- 
al to the color excess E, which is the difference between the photo- 


2 Ap. J., 86, 176, 1937. 
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graphic and the visual absorption and can be ascertained by ob- 
servation. 
A(r)=x-E, (2b) 


where the constant x is the ratio of photographic absorption to color 
absorption. 

The second assumption is probably a much better approximation 
than the first, since it takes some account of the irregularities in the 
absorbing medium. The data on color indices in star clusters are, 
however, at present not sufficiently complete or sufficiently homo- 
geneous for its application. In either case the constants a or x must 
be evaluated empirically with the aid of one of the other methods of 
distance determination which is not affected by interstellar absorp- 


tion. 


Ill. DETERMINATION OF DISTANCES AND OF THE ABSORPTION 
COEFFICIENT FROM ANGULAR DIAMETERS 

The diameter method is based on the principle that star clusters 
which have the same constitution C have the same linear diameter 
D, (in parsecs), whatever their distance from the observer. Clusters 
of the same constitution can be grouped together by a classification 
which takes into account the number of cluster stars as well as the 
degree of central concentration. If the angular diameter d (in min- 
utes of arc) of a cluster has been estimated by the observer, the dis- 
tance r is obtained by 


r= — 3438 


or 
log r = log D, — log d + 3.536. (3) 


For each class of clusters the value of D. has to be found by means 
of clusters with known distances. The photometric distances of the 
nearer clusters for which the absorption is negligible might be used 
for the purpose; but the number of such objects is rather too small. 
If the assumption in equation (2a) of uniform absorption is adopted, 
the constant a and the values of D. can be determined by a combina- 
tion of the photometric method and the diameter method. 
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Each cluster furnishes two equations of the forms (1) and (3): 


4 


slogrta-r =yt+5 


(4) 
—log r + log D, = log d — 3.536, | 

in which the terms on the right-hand side are given by observation. 
For n clusters divided into k classes there are thus 2” equations with 
n unknown distances r; k unknown D,’s; one unknown a. Since k is 
much smaller than , there are fewer unknowns (7 + k& + 1) than 
there are equations of condition, and a statistical solution is possible. 
Owing to the transcendent nature of the first equation, however, the 
solution has to be accomplished by successive approximations. 

The first step is to find the relative values of the mean linear 
diameters 

D. 


Di, = D, (5) 


of the different classes, where the constant D, represents an unknown 
scale correction to the adopted relative diameters D/.s When the 
value of log r from equation (3) and of D. from equation (5) are in- 
troduced into equation (1), the latter becomes 


d 
A(r) +c= yt 5 log — 12.68, (6) 


where 


c= 5 log i; oe 


It is thus possible to compute for each individual cluster the amount 
of photographic absorption, except for a constant zero correction on 
which the scale of the distances by the diameter method depends. 

For 97 clusters? for which reliable and homogeneous data are 
available, the values of A(r) + c have been plotted in Figure 1 as a 
function of the distance modulus y, which is also a function of r and 
is obtained by the photometric method. The vertical scattering of 

3 In the present solution these were taken from Lick Obs. Bull., 14, 160, 1930. 

4 Of the 100 clusters listed in 7bid., p. 156, Table 3, many of the southern clusters 


with incomplete data were omitted, while clusters with more recent observations were 


added; about one-third of the clusters are new or have new data. 

















GALACTIC STAR CLUSTERS IgI 


the points is, in the first place, due to observational errors in the 
estimated angular diameters and in the classification which enters 
into Di. As these are approximately percentage errors, the results 
for A + c are practically of the same accuracy for near and distant 
clusters. The errors in y, as far as they are due to incompleteness of 
data or statistical uncertainty of the photometric method, are con- 
siderably smaller. Part of the scattering of points is probably due 
also to irregularities in the absorption. Normal points found for 
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Fic. 1.—Photographic absorption of star clusters as function of distance modulus. 
The abscissa represents the distance modulus obtained by the photometric method, the 
ordinate the absorption (plus an unknown zero correction c) derived from equation (6). 
Each cluster is marked by a full dot, group means by open circles. The full curve shows 


the relation between absorption and distance modulus in case of a uniform absorption 
with a coefficient a = o™7o per kiloparsec, the dotted curve for a coefficient a = 1™5 
per kiloparsec. The unknown zero correction c is for both curves so determined as to 


fit the nearer clusters. 


groups of 12-20 clusters are marked by open circles; the full-line 
curve shows the relation between A + ¢ and y for uniform absorp- 
tion and a coefficient 

a = o™70 per kiloparsec , 
which gives the best fit with the observations. The dotted curve in- 
dicates the case of uniform absorption with a coefficient a = 1™5 per 


kiloparsec; for distances greater than 700 parsecs such a large coeffi- 
I 7 
cient seems quite incompatible with our data. 
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In a similar way it would be possible with the assumption of equa- 
tion (2b) to determine the unknown ratio x and the values of D, by 
a combination of the photometric method and the diameter method. 


IV. DISCUSSION OF SYSTEMATIC ERRORS 
To obtain a correct result for the absorption coefficient by this 
method it is essential that the estimates of angular diameters and 
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Fic. 2.—Areal star density in Praesepe and NGC 6645. The number of stars per 
S I 


unit area (ordinate) is plotted as a function of the angular distance (abscissa) from the 
center of the cluster. The dotted horizontal line indicates the mean density of the 
background stars; L, the estimated limit of the cluster. 


classification be homogeneous and free from systematic errors de- 
pending on the distance of the objects. 

In estimating the angular diameter, we have attempted to find the 
diameter of the smallest circle outside of which no excess of cluster 
stars over the background stars is noticeable by mere examination of 
a photograph or chart. It is of interest to compare the results of such 
estimates with those of star counts. Figure 2 gives such a comparison 
for Praesepe and NGC 6645, two clusters of similar constitution 
(class 1r) but at very different distances. From star counts in con- 
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centric rings drawn around the center of each cluster the star density 
per unit area was derived and plotted against the angular distance 
from the center. For Praesepe the counts were made on the Bonner 
Durchmusterung chart, for NGC 6645 on a photograph reaching the 
photographic magnitude limit 15. The ‘‘estimated” limit of the clus- 
ter is indicated by the vertical line at Z and is considerably inside of 
the true limit shown by the star counts. The scale of the abscissas is 
so chosen that the estimated limiting radius has the same length for 
both clusters, and the circle with this radius is taken as unit area. 
The scale of ordinates is made inversely proportional to the total 
number of cluster stars within the limiting circle, and the areal den- 
sity of the background stars is given by the dotted horizontal line. 
The two characteristic curves obtained are remarkably similar for 
the two clusters. 

This and other tests show that if the observer enlarges different 
clusters of the same class to approximately the same size and in- 
creases the magnitude limit of the photographs accordingly, the 
estimated angular diameters are well comparable and are little af- 
fected by the distance of the cluster. Since the conditions stated 
have, on the whole, been realized, a direct effect of distance on the 
estimates of diameter and central concentration such as suggested by 
Ten Bruggencate’ is out of the question. 

It will be noticed in Figure 2 that the observer’s estimate over- 
looks that part of the cluster where the excess of cluster stars be- 
comes less than 50 per cent of the background stars. The richness of 
the background on which a cluster appears projected may, therefore, 
have an effect on the diameter estimates. Local irregularities in the 
star field of the Milky Way affect the statistical data like accidental 
errors of observation. But the question we have to investigate is 
whether the number of background stars within the limiting circle 
increases with increasing distance of the observer when the magni- 
tude limit is raised according to the distance. 

The example of Praesepe will illustrate this point. On a photo- 
graph with the apparent photographic magnitude limit 9.6, Praesepe 
will appear as a group of 35 stars occupying a circle go’ in diameter, 
and such a field will contain 9 additional background stars. If we 


Naturwissenschaft., 18, 725, 1930. 
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imagine this cluster removed ten times farther away (distance 1500 
parsecs) along the galactic plane, it will reduce to an area of 9’ in 
diameter, and to show the same 35 cluster stars our photograph must 
reach the magnitude limit 15.7 (assuming an absorption of o™8 per 
1000 parsecs). In low galactic latitudes the average number of stars 
brighter than 15™7 in a circle of 9’ diameter is 30, and the 35 cluster 
stars will now appear projected on a field containing 30 background 
stars. 

If the same upper limit of absolute magnitude is reached in all 
clusters, we find in low galactic latitudes the more distant clusters 
projected on a richer background. The angular diameters of distant 
clusters are thus probably estimated somewhat too small as com- 
pared with the nearer clusters. On the other hand, a richer back- 
ground would also lead an observer to underestimate the number of 
cluster stars and would affect the classification. As the adopted 
linear diameters D. are smaller for the poor clusters than for the rich 
clusters, an underestimate of the angular diameter d for distant clus- 
ters would in many cases be accompanied by the assignment of too 
small a linear diameter D.. The ratio D, : d which determines the dis- 
tance by the diameter method will thus, on the whole, be much less 
affected. If the two effects do not quite compensate each other, it is 
more likely that the diameter distances of the distant clusters are 
slightly too large and the absorption coefficient a derived therefrom 
slightly too small. We have therefore retained the value a = o™8 per 
kiloparsec, with which the distances of Lick Observatory Bulletin, No. 
420, were computed. 

Van Rhijn® suggested that the clusters with spectroscopic data 
suitable for such an investigation are necessarily selected so as to 
favor objects seen through the more transparent regions of our stellar 
system. This is undoubtedly true for clusters more distant than 1500 
parsecs, and the average photographic absorption coefficient in low 
galactic latitudes may well be as large as 1™o per kiloparsec. Even 
if such a selection is admitted, the absorption coefficient determined 
from the clusters is that which applies to these objects on the average 
and must be used for the correction of the photometric distance 
measures. 


6 Groningen Pub., No. 47, 1936. 
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V. GALACTIC ROTATION SOLUTION 

The third method, based on the differential galactic rotation effect 
in radial velocity, can now be used to compare our distance scale of 
star clusters with that of other galactic objects (B stars, Cepheid 
variables, planetary nebulae). Although the radial velocity meas- 
ures in 75 galactic star clusters are far from complete, preliminary 
data are at present available for 39 clusters with distances ranging 
from 430 to 4600 parsecs. 

The observed radial velocities of the clusters, corrected for solar 
motion (19.5 km/sec toward A = 18%o, D = +30°), which we shall 
designate as V, were represented by the formula 


V = k, + kz [r cos? 6] + A[r cos? b sin 2(1 — 325)], 


in which ¢ is the adopted distance (mean of photometric and diam- 
eter methods), / and 0 are the galactic longitude and latitude of each 
cluster, while the longitude of the galactic center was assumed to be 
325°. The last term represents the effect of differential galactic rota- 
tion when higher harmonics are neglected. The k term was assumed 
to consist of a constant k, and a term proportional to r cos*b. The 
determination of three unknown constants by a least-squares solu- 
tion furnished the following values: 


A = +0.0150 + 0.0010 km/sec + parsec , 
k, = —3.5 km/sec , 
kz = —0.0043 km/sec - parsec . 


The residuals of the observations from the formula give a disper- 
sion (standard deviation) of +12 km/sec for the peculiar velocities 
of galactic star clusters. That the value of k; mostly represents a 
systematic error in our preliminary radial velocities is confirmed by 
observations of bright standard stars. The term with k,, however, 
appears to be real and suggests a local contraction along the galactic 
plane of the observed star clusters. A similar change of the k term 
with increasing distance is indicated by the radial velocities of O and 
B stars,’ and it is remarkable that the galactic rotation solutions of 


7 See Paris Pismis, Pub. Istanbul U. Obs., No. 10, 1938. 
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other distant galactic objects, such as planetary nebulae® and Ceph- 
eid variables,’ also lead to negative k terms. 

Table 1 lists the results obtained by various observers for the ga- 
lactic rotation constant A (third col.) and for the average coeflicient 
of photographic absorption a per kiloparsec (fourth col.). The second 
column indicates the class of objects studied, and the last column, 
the methods used for their distance determinations. Our result for A 
is in close agreement with that of Plaskett and Pearce’® and with that 


TABLE 1 


Author | Objects { a Distance Scale 


(km/sec + par- | (per 


sec) kpc) 
Plaskett, Pearce*| O-B3 stars 0.0155 +0.0009 Proper-motion parallaxes 
OMSEET «sos < + 34:00} RRS o1g +0.003 | 1™8 | Proper-motion parallaxes 
TOV E..« Cepheid vari- 0209+0.0008) 0.85) Period-luminosity rela- 
ables tion 
Berman$..... Planetary nebu- O140 0.72} Proper-motions; photo- 
lae metric; angular diame- 
ters 
Trumpler.... Galactic star 0.0150+0.001 | 0.70; Photometric; angular di- 
clusters ameters 
* Pub. Dom. Ap. Obs., §, No. 3, 1033 t Mt. W. Conir., No. 607, 1939. 
t For A: B.A.N., 4, 82, 1927; for a: ibid., 8, 233, 1938 §$ Lick Obs. Bull., 18, 64, 1937 


of Berman but differs considerably from that of the Cepheid vari- 
ables. 

A former galactic rotation solution by Miss Hayford," based on 
25 galactic clusters, yielded a considerably smaller value for A 
(0.009 km/sec - parsec), and this has induced some investigators” to 
suggest a scale error in our cluster distances. The difference between 
the result of the new solution and the former one is partly due to the 
additional and improved data and partly to the different treatment. 
Miss Hayford’s solution included a constant k term as well as the ele- 
ments of the solar motion and led to a large value of 30 km/sec for 
the solar velocity. Our new solution adopts the usual value of the 

8 Lick Obs. Bull., 18, 64, 1937. 10 Pub. Dom. Ap. Obs., 5, No. 3, 1933- 

9 Mt. W. Contr., No. 607, 1939. 11 Lick Obs. Bull., 16, 53, 1932. 

12 Raimond, Groningen Pub., No. 46, p. 32, 1934; Oort, B.A.N., 8, 245, 1938. 
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solar motion but introduces a k, term proportional with the distance 
in addition to a constant f;. 


VI. THE OBSERVED SPACE DISTRIBUTION OF THE CLUSTERS 
The distances of all known galactic star clusters (334) which are 
sufficiently well defined were derived by either or both of the meth- 
ods outlined in sections II-III, and these data permit a study of the 








Fic. 3.—Distribution of known galactic star clusters in the projection on the galactic 
plane. Within the circle of dark shading the clusters are nearly uniformly distributed; 
in the ring of light shading they thin out rapidly. The large dotted circle drawn around 
the sun indicates the limit of 10,000 parsecs to which a search of clusters might be ex- 


tended; the full circle shows the assumed limit of the galactic system. 


space distribution. Figure 3 gives an outline of the region (shaded 
area) occupied by the known clusters as compared with the outline 
(large circle) of the galactic system according to the generally ac- 
cepted hypothesis. The projection of the individual clusters on the 
galactic plane is illustrated in Lick Observatory Bulletin, No. 420, 
p. 179. Essentially, the clusters group themselves around a center 
situated near the sun, and their arrangement shows no relation to 
the general structure of the galaxy. Up to a distance of 2000 parsecs 
from the sun the clusters are in their projection on the galactic plane 
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fairly uniformly distributed. As the distance increases beyond 2000 
parsecs they thin out rapidly, and at a distance of 5000 parsecs we 
practically reach the limit of the known clusters. 

The question immediately imposes itself as to whether these fea- 
tures in the apparent distribution mainly reflect the incompleteness 
of our present cluster list. To admit such an interpretation it must 
be shown, in the first place, that even the richest and most luminous 
clusters would not be observable at distances greater than 5000 par- 
secs with the instruments used for the search, such as Herschel’s re- 
flector or the Franklin-Adams camera. In the second place, an ex- 
planation has to be given of how a large number of clusters could 
have escaped discovery even at distances of 2000-3000 parsecs. 

The well-known local clouds of obscuration (dark nebulae) in the 
Milky Way are undoubtedly obstructing the view of some of the 
more distant clusters, but their effect can hardly explain the rapid 
thinning-out, starting at 2000 parsecs, or the lack of clusters beyond 
5000 parsecs. The existence of a general absorption of light in a 
layer stretching along the galactic plane will, on the other hand, 
greatly affect all distant clusters. Although this layer appears to be 
of considerable irregularity, we have assumed an average rate of o™8 
per 1000 parsecs for the absorption of photographic light. 

To study the extreme distance limit to which galactic clusters are 
observable on Milky Way photographs with a given equipment, we 
have selected three of the richest and most luminous clusters: Prae- 
sepe, the Pleiades, and h Persei. The data of Table 2 illustrate how 
these clusters would appear on a photograph containing stars bright- 
er than magnitude 18, if seen from distances of 3000, 5000, and 7500 
parsecs. The third column of the table shows how much fainter the 
cluster stars would be than at their actual distances. The fourth 
column gives the apparent angular diameter, column 5 the number of 
cluster stars brighter than eighteenth magnitude, and column 6 the 
number of foreground or background stars within the field of the 
cluster. Column 7 indicates the percentage of all stars within the 
area of the cluster which are cluster members; a cluster would prob- 
ably not be recognized as such if this percentage is less than 25 or 30. 
The magnitude of the brightest cluster star is found in the last col- 
umn. 

Praesepe would evidently no more appear as a cluster at a dis- 
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tance of 5000 parsecs; the Pleiades would still be recognized at 5000 
parsecs but not beyond; h Persei might possibly be noticed at a dis- 
tance of 6000 parsecs. None of the systematic surveys of the Milky 
Way, such as the Franklin-Adams chart, reaches the eighteenth- 
magnitude limit, and these most favorable examples leave little 
doubt that clusters more distant than 5000 parsecs would have been 
discovered in only a few exceptional cases. 

We have now to examine the effect of distance on some of the 
poorer clusters containing no stars of high luminosity. The star clus- 

TABLE 2 


APPEARANCE OF CLUSTERS AT GREAT DISTANCES 


(Absorption o™8 per 1000 Parsecs) 








| | . . : | 
| | | | 
mpg < 18 
| ESTIMAT | : 
| Per Cent} Bricut- 
A DISTANCE | Am | ED ANGI ie : 
CLUSTER | ss | | CLUSTER | EST CLUS- 
| (PARSECS) (Pc) | car D1 ; Back- : 
| Cluster STARS TER STAR 
| AMETER ground 
| | Stars . | 
| | | Stars 
| , . | ’ ; j iss si ap or 
Praesepe 3000 | 8.8 | 4°5 26 44 | 37 15.0 
5000 | I1.5 | 2.9 2 16 Be winwate ws 57.9 
; 3000 | 8.8 | 6 49 78 ( fy 
Pleiades 45000 | I1.5 | 3.6 15 28 5 14.4 
7500 | 14.4 2.4 | I 12 37.3 
h Persei 5000 5.9 5.2 07 59 62 12.7 
7500 | 8.8 | KEL. 5 37 16 15.6 
| 











ter in Coma Berenices is the best-known example of this type. 
Table 3, in an arrangement similar to Table 2, describes the appear- 
ance of the cluster as seen at various distances, projected on an aver- 
age star field in low galactic latitude. Three photographs with limit- 
ing magnitudes of 18, 16, and 14 are assumed to be available. Al- 
though the five brightest cluster stars would be visible on the first 
photograph at a distance of 4000 parsecs, they would be lost among 
the great number of background stars. On the second photograph 
the percentage of cluster stars is somewhat more favorable for dis- 
tances of 2000-3000 parsecs, and on the third photograph the cluster 
would be fairly noticeable up to a distance of 1000 parsecs. On a long 
exposure a poor cluster at great distance is “drowned” by the back- 
ground and foreground stars; on a shorter exposure the cluster stars 
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are only visible if the distance is not too great. From the data of 
Table 3 we conclude that the Coma Berenices cluster might be recog- 
nized as a clustering to a distance of 1500 or 2000 parsecs. For poor 
clusters without stars of high luminosity the distance limit to which 
they are noticeable is mainly conditioned by the rapid increase of 
background stars with increasing magnitude limit, and most of such 
clusters would be overlooked at distances greater than 2000 parsecs. 

The observed space distribution of galactic clusters seems, there- 
fore, satisfactorily explained by incompleteness and selection of the 
known objects. 

The absorption of light in the vicinity of the galactic plane is the 
most serious obstacle in any search for more distant galactic clusters. 


TABLE 3 


APPEARANCE OF COMA BERENICES CLUSTER AT DIFFERENT DISTANCES 


™ £15 M 1¢ m I 
ANGu- BRIGHT PE - = = pe a _ J = 
DISTANCE Am LAR EST | | 
(ParsEcs) (PcG) | DiaAmM-/CLustER Clus- | Back P Clus- | Back Per | Clus- | Back- | p,. 
ETER | STAR ter | ground y st ter | ground | , ta | ter | ground | ¢ a 
Stars Stars | f Stars Star | Stars Stars sis 
1000 6.4 | 18’ 11.6 31 70 3 147 14 I 26 2 
2000 8.7 rs) 13.9 I 13 I } a 3 23 0 
3000 Io ) 15.¢ 12 8 I i ae 
4000 Ir.d } 17 I | 


As this absorption is greater for photographic light than for longer 
wave lengths, it is of interest to investigate whether observation in 
yellow or red light would be more advantageous for the discovery of 
distant galactic objects. That this is the case is clearly indicated by 
the data of Table 4, which describes the appearance of Praesepe, the 
Pleiades, and h Persei at distances of 5000~-10,000 parsecs, as seen 
on photographs taken in red light (AX 6000-6500), reaching the sev- 
enteenth red magnitude. On such photographs, clusters like the 
Pleiades may be recognized at a distance of 7500 parsecs, and such as 
h Persei even somewhat farther. The red magnitude limit 17 is thus 
more effective than the photographic magnitude limit 18, because 
the highly luminous cluster stars at such great distances are much 
reddened by selective absorption, much more than the background 
stars, which are in reality mostly foreground stars. 

It is a remarkable fact that the introduction of photographic ob- 
serving methods during the last fifty years has added very few new 
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galactic clusters to the large list of those discovered by the Herschels. 
Their use of visual observations evidently gave them an advantage 
over ordinary photographs utilizing light of short wave lengths. 

An extension of our knowledge of clusters to a distance of 10,000 
parsecs seems, at present, quite feasible if the Milky Way belt is pho- 
tographed in red light with an instrument of sufficient resolution, 
light-gathering power, and a large field. A large Schmidt reflector 
would seem particularly adapted to such purpose and may possibly 
record stars as faint as the eighteenth or nineteenth red magnitude. 

TABLE 4 
APPEARANCE OF CLUSTERS AT GREAT DISTANCES IN RED LIGHT 


(Absorption 0"45 per 1000 Parsecs) 








| 
| | mred < 17 
| . 
| Estimat-|_ | en Coes BRIGHT- 
: DISTANCE | ED ANGU- | as “| Est CLus- 
CLUSTER | | CLUSTER 
PaRSECS)] LAR D1 Back- . TER 
| Cluster STARS . 
AMETER ground STAR 
| Stars } 
| Stars | 
| 
| 
Praesepe 5000 7 ey a 6G: | «4 15.6 
Pleiades 7500 2.4 | 7 = fF a 14.9 
i isinaad 7500 |} 3.5 | 26 28 48 13.8 
: 10000 | Oo | 4 15 21 15.5 











The diagram of Figure 3 gives a good idea of the results that might 
be gained by a systematic search for distant galactic clusters. A 
study of their space distribution to a distance of 10,000 parsecs from 
the sun (dotted circle) would bring out some of the significant struc- 
tural features of our galaxy. In galactic longitudes 80°-200° the 
limit of the system could be traced, and in this direction we would 
not expect an appreciable number of more distant clusters. A point 
of great interest would be to ascertain whether the region of the ga- 
lactic center is devoid of clusters. Such a possibility is suggested by 
an examination of the nearer spiral nebulae of intermediate and later 
types, in which nuclei of condensation or clusterings are mainly found 
in the spiral arms but not in or near the central nucleus. If our pres- 
ent views on the dimensions of our galactic system are correct, we 
must expect a great number of undiscovered clusters in galactic 
longitudes o°—7o° and 220-290. 

UNIVERSITY OF CALIFORNIA 
November 3, 1939 








THE CONTINUOUS ABSORPTION OF LIGHT BY 
NEGATIVE HYDROGEN IONS 


H. S. W. MASSEY AND D. R. BATES 


ABSTRACT 

In view of the possible importance of negative hydrogen ions in absorbing long 
wave-length radiation in stellar atmospheres, it is desirable that a reliable theory of the 
continuous absorption by such ions be available. Such a theory is provided in this 
paper, and the resulting calculated values of the absorption coefficient in the wave 
length region of astrophysical importance are shown to be accurate. This is made pos- 
sible, on the one hand, by the existence of the accurate wave function for the ground 
state of H- due to Hylleraas and, on the other, to a detailed analysis of the factors de- 
termining the form of the continuous wave function representing the state of the system 
in which one electron is ejected. The effects of distortion of the wave representing this 
electron by the static and polarized atomic field, of electron exchange, and of the non- 
separability of the wave function are all considered in detail and are found, for the par 
ticular case of H~-, to be incapable of affecting the calculated values. 

It has been pointed out by Wildt' that, in discussing the absorp- 
tion of radiation in the outer layers of stellar atmospheres, account 
must be taken of the possible existence of negative ions in such 
regions. In particular, in an atmosphere including both meta! and 
hydrogen atoms the ionization of the former can supply electrons 
necessary for the formation of negative hydrogen ions. The energy 
of attachment of an electron to a hydrogen atom is 0.70 e. volts? 
(there is only one stationary state of the attached electron). As a 
result the presence of H~ ions provides a source of continuous ab- 
sorption of long-wave-length radiation (wave No. 5674 cm~ and 
upward) and Wildt? and Strémgren* have shown that this suggests 
a way out of certain difficulties concerning the magnitude and the 
wave-length dependence of stellar absorption coefficients. An accu- 
rate theory of the continuous absorption of light by H~ ions is there- 
fore of considerable interest, and it is fortunately possible to develop 
such a theory because of the relative simplicity of the problem, 
which involves the interaction of two electrons only. Calculations 
in this direction have been already made by Jen> and by Massey 

tA p. J., 89, 295, 1939. 

2 Hylleraas, Zs. f. Phys., 60, 624, 1930. ‘ Unpublished. 

s Ap. J., 90, O11, 1939. 5 Phys. Rev., 43, 540, 1933 
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and Smith,° but neither result can be regarded as indicating more 
than the order of magnitude of the effective cross-sections involved 
until a more detailed analysis is employed. Jen, although using an 
accurate wave function for the negative ion, did not use the correct 
form for the continuous wave function of the detached electron, 
while Massey and Smith, although using the correct form of con- 
tinuous wave function, used only a rough approximation for that of 
the ion. In this paper the problem is treated accurately as one in- 
volving two electrons, and a number of effects not previously taken 
into account in the calculation of absorption coefficients are dis- 
cussed in detail and are found to lead to no appreciable modification 
of the results. 
GENERAL FORMULATION OF THE PROBLEM 

Consider a stream of light quanta of definite frequency traversing 
a block of matter containing V atoms/cm. The fraction of quanta 
absorbed in a distance 6x is given by VQ6x, where Q is the effective 
area or absorption cross-section for absorption of a quantum by an 
atom. The intensity of the beam in traversing a distance x is there- 
fore diminished by a factor e~*?*, and the absorption cross-section Q 
is identical with the atomic absorption coefficient of the material. 
Our problem is to determine ( for H7~ ions as a function of the fre- 
quency v of the incident quanta. This involves determining the 
probability of a photoelectric transition from the normal state of the 
ion to one of its continuous states. If £, is the electron affinity of the 
ion (the energy excess of the neutral atom over that of the negative 
ion), and £ the energy 3mv? of the ejected electron, we must have 


E=h-—E,. 


If &(r,, r,) is the wave function of the ground state of the negative 
ion and W;,(r,, r,) that for the state of energy E + Ey (where Ez is 
the energy of the neutral hydrogen atom in its ground state), the 
absorption cross-section is given by the usual quantal formula 


3274m’e? 


(j= 
7 3hsc 


vv | ff@*(n, r2)(1 + r2)Ve(n, 12)dridr, |? , (1) 


6 Proc. R. Soc., Ser. A, 155, 472, 1936. 
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provided ® and Wz are suitably normalized, i.e., ® to unity and VW; 
to represent an electron bound to the proton, together with an out- 
going current of unit density. This formula can certainly be relied 
upon to give accurate results provided the wave functions ®, VW; are 
accurately known. For most absorption problems the accurate de- 
termination of wave functions is a matter of some difficulty, but in 
the particular case of H~ we can obtain these in a satisfactory man- 
ner by methods which we now proceed to discuss. 


DETERMINATION OF WAVE FUNCTIONS: THE GROUND- 
STATE WAVE FUNCTION OF 1 
A thorough study of the structure of H~ has been made by 
Hylleraas,’ using the variation method which has been so successful 
in determining the energy in the two-electron problem of the helium 
atom. The wave function® which he gives has the form 


to 
— 


P(r, T2) = Ne~e1ta) } I + V(r: = $2)" + Bry. ) ( 


where a, 8, y have the values 0.770, 0.308, 0.119, respectively. The 
term JN is the normalizing factor given by 


' 


oy 7 2 a Ray 2 ) —(1/2 
a. 3 + 66° 77BY OY jr ett 


8a a? 8a3 a‘ 


a3 r(1 + 


and 


The corresponding wave function for the helium atom leads to a 
total binding energy of the two electrons equal to 78.45 e. volts, dif- 
fering only very slightly from the observed 78.48 e. volts.’ It is very 
unlikely, therefore, that the form in equation (2) is seriously in error 
in those regions of r,, 7, which contribute appreciably to the energy 
of the negative ion, i.e., in the region where ®(r,, 7,) is large. One 
cannot be so sure of its accuracy in other regions (particularly, at 
large distances from the nucleus), but it so happens that in the fre- 
quency range of astrophysical importance, the main contribution to 

i Loc. cit. 

8 Throughout this paper the wave functions are given in atomic units. 


9 Baber and Hassé, Proc. Cambridge Phil. Soc., 33, 253, 1937. 
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the integral in equation (1) comes from such values of 7,, 7, that we 
can be confident that the function is satisfactory for our purpose. 


THE CONTINUOUS WAVE FUNCTION 


In the absence of any interaction between the hydrogen atom and 
the ejected electron, we would have 


WVe(n, r,) = 2 2)! xs0(71)e***s + X10(72)e**1} ’ (3) 


where x: is the ground-state wave function of the hydrogen atom, 
nr /e-, and k is the momentum of the ejected electron. The sym- 
metrical wave function alone concerns us since the ground state of H~, 
being a singlet, is also symmetrical in the space co-ordinates of the 
two electrons. This simple approximation can be shown to give re- 
sults of quite high accuracy for the problem with which we are con- 
cerned, so we will consider first the calculation of the absorption co- 
efficient with its use and then show that the corrections introduced 
by higher approximations to VY, are unimportant. 


CALCULATION OF ABSORPTION COEFFICIENT 
We have, using equations (2) and (3), 


3 214m?7e? 


0 = vo|I|?, 


3hsc 
where 


I =(2r)—-@/N ffe-erta) {1 + y(r, — 12)? 
+ Brij (ni + 12) (et*4a-1 4. e‘*:a~") dr, dr, > 


Using the expansion 


a tie. “ae 
e tke = ( 1"(2n + DS a1 (kr) P,(cos 6) ’ 
— 
n 


2kr 


and carrying out the angular integrations involved, we can readily 


see that 


I = (I, + BIp)16m?NR-@”) 
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where 
© co 
T4 = f f é ~er,—(e-ta)raf I + v(r; — r2)?\r3/*r3 J) (kri)dridr, ’ 
U 0 
© co 
Ip aS f J e~elrr tray 3? { o4( ry, rae, "/2) J, (rs) 
0 rt) 
+ prs, retry /)S ynlkrs) \dridrs , 
with 
grz7t+ rz | 
Pos, r2) = : aoe 
3M 
fs > NG 
"2 5rils 
pi(r:, r,) = ee 5 ile 
5!1 ) 
Since 


2 \*s sin kr 
J 3(kr) = (=.) (cos kr — | ; 


the integrations over r,, 7, involved in 74 and Jz can be performed 

without difficulty to give 

I = (204873)/*Nosk~[6yp4(4kp, cos 5: — sin 4£:) + 2(8 — 6ya)p} 
(3kp; cos 4; — sin 3&:) + (1 + 12707)p?(2kpx cos 3&: — sin 2€,) 
— 2B{op?(2kp, cos 3&2 — sin 2§.) + 607p,(kp, cos 2& — sin &) 


+ 12 o5p,k cos £ — 12004 — £2)}],, 


where 
p; = (Rk? + a?)—@), py = {k? + (1 + 2a)?}-/2) , 
k k 
g, = arctan -, £, = arctan ———— , 
a I+ 2a 
o =(1+ a)". 


In Table 1 the numerical values of the continuous absorption co- 
efficient for H-, calculated for a series of wave numbers of the inci- 
dent quanta, are given. In the range of importance in astrophysics 
it will be seen that Q is in the neighborhood of its maximum. Before 
proceeding to discuss the errors arising from the use of the simple 


~~ 
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form in equation (3) for the continuous wave function, it is of inter- 
est to compare our results with those obtained previously by Jen" 
and by Massey and Smith." Jen obtains results somewhat similar to 
ours for the higher frequencies, but his absorption coefficient tends 
to a finite value at the low-frequency limit, whereas ours tends to 
zero. This difference arises from our use of a plane wave to repre- 
sent the behavior of the ejected electron, whereas Jen uses a wave 
modified by a Coulomb field of potential — 0.4e?/r. For small en- 


TABLE 1 


CALCULATED ABSORPTION COEFFICIENTS FOR H- 


| | 


Wave Number | Absorption | Wave Number | Absorption 
} 


of Radiation in | Coefficient in || of Radiation | Coefficient in 


103 Cm~! 10-'?7 Cm? || in 103 Cm7! 10~'? Cm? 
| . 
5.674 (head) 0.00 ee or See : 
ea ea } °o.86 | 45 1.66 
SEs... 1.87 | 50 1.41 
20... 2.49 60 i I.02 
ee 2.62 70 nets 0.76 
30.. 2.49 80 0.56 
Ce 2.26 | QO... 0.42 
40 1.94 Foo... . er 0.33 


ergies of the ejected electron the major contribution to the integral 
in equation (1) comes from large distances. At such distances his 
assumption of a Coulomb interaction between the neutral hydrogen 
atom and the second electron becomes grossly in error, and he ob- 
tains too large a value of Q. On the other hand, our results differ 
from those of Massey and Smith in the other direction. They find 
()—o0 at the low-frequency limit as do we, but their maximum 
occurs at much higher frequencies. This is due to their use of a 
rough approximate form for the ground-state wave function, con- 
siderably more compact than that given by equation (2). 


EVALUATION OF CORRECTING TERMS 


As the accurate determination of the continuous wave function 
WV(1;, 7.) is a problem essentially the same as that which arises in 


10 Loc. cit. tt Loc. cit. 
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discussing the collisions of electrons with hydrogen atoms, so it is 
convenient to adopt the method of collision theory." We write 


Ve(r,, t2) = {Q(r, m2) + Q(r, m1) } 2-8/2) , (4) 
where Q(r,, r,) has the asymptotic form for large r, 
etky (rz) . 


The expression Q(7,, r,) must satisfy the wave equation 


(Vit V) +E+—+———[0%n,1) = 0. (5) 
1; 2 r 


h? 


812m 


where E is now the total energy of the system. 
As it must be a proper function of 7,, 7., it may be expanded in a 
series 


SE ae T2) = (Za + )xnilt2)Frilts), (6) 


where the xn(r.) are the sequence of wave functions for the different 
states of the hydrogen atom and so satisfy the equation 


~~ 
~— 


he 2 : 
( Vit Ext = )xu =o. ( 
2 


8r2m ~ 


The integral sign arises from the contribution of the continuous 
wave functions of the hydrogen atom. The suffixes m/ represent 
the total and azimuthal quantum numbers, respectively, of the 


discrete hydrogen states, so 
X nil T2) = Uni r2) V (2, o2) A 
where Y; is a spherical harmonic of order /. Expanding 
Fri(t:) = U(2s + 1) Fils) ¥ (61, $1), 


t2 Mott and Massey, The Theory of Atomic Collisions, chap. x, Clarendon Press, 1933. 
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and substituting in equations (6) and (1), we find, since 
f{VsV. sin 0d6dd = o, $ 9s’, 
a. ae —. 
a2s+1’ 
(friV.(01, 1) O(, 2) Vel, 12)drdr, 
= (xasynen [Of Mri a(rytra)((r + yr — 72) Dn { Uno( 2) F1,(11) 


+ tni(r1)Fe,(r2)} + BUnZif oils, 72) Uni(r2) 
(uF (ri) + oF (rs)) + ttni(rs) (uprss(rs, 72) Fee *(r2) 
+ vpri(rr, r2)Fi7*(r2))}|dridr. + I. , 


where p, v are numerical coefficients 


and 
ro = T(25 + 1)p,(7:, 72)P.(COS O,:) « 


[. represents the contribution from the terms of Vz(7%,, 7.) involving 
the continuous wave functions of the hydrogen atom. It is apparent 
that only certain terms in the expansions of the functions F,,; give 
nonvanishing contributions to the absorption coefficient, and this 
simplification is important. To determine the functions /’,, we sub- 
stitute equation (6) in equation (5) and use equation (7) to give 


e? e? 


h? 


(Snr t Pxnilt) Vi+tE— Ey) Fialt:) = (= = ar, r.). (8) 


ry 


8r?m 
Multiplying both sides of this equation by x,,(r.) and integrating 
over r, gives, on account of the orthogonal properties of xn2(1.), 


: VV+E—E, ) Fit) = ef(; ape Joan, r.)x*(12)dt2. (9) 


Sa?7m ri2 r; 


In view of the form of equation (6) for 2, this gives an infinite set of 

simultaneous integro-differential equations for the functions Fy. 

Solutions are required which have the asymptotic form for large 1, 
Fyo(t,) ~ et + rete ®fio(81, bx) , 


Fil) ~ rite knfn(Or, d:) , (n,l # 1,0), 
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where 


; 81?m 
fo 





iB). 


The functions F,,,(7,) correspond to electron waves scattered after 
exciting the atom to the state m/l. In the energy range of importance 
in our problem the ejected electrons have insufiicient energy to pro- 
duce excitation. This does not mean that all the F,:(r,) vanish, but 
they have the asymptotic form 


Ce Yn" fil 1, ¢1) ) 


where 72 = — k?. 
Terms of the form 


Xe(12)F (rs) , 


where x.(r,) is a continuous hydrogen wave function, represent 
capture of the outgoing electron and ejection of the atomic electron. 
This corresponds to electron exchange, and it is possible to take 
such terms into account directly without explicit introduction of the 
cumbersome continuous hydrogen wave functions. For convenience 
we shall refer henceforward to the term of Q involving F,,(r,) as the 
bound-free term, those of form x.(r.)F.(r,) as the free-bound terms, 
and the remainder as bound-bound. Then we can say that the 
bound-free term includes the unperturbed plane wave and the modi- 
fication of this wave by the atomic field (including polarization), the 
free-bound terms represent exchange phenomena, and the bound- 
bound arise from the nonseparability of Q. 


THE BOUND-BOUND TERMS 
Our approximation in equation (3) amounts to taking 
Fi(ts) = 0 ; n,t#1,0; F,.(1,) = e** 


To obtain the next approximation we substitute these values on the 


right-hand side of equation (g) giving 


l_- 





a 82?me? I 
(V? — ya)Fu(ti) = ref (E— FE )etxlrdxairddrs § (10) 
WN 
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so, using Green’s function, we should have 


F,(r) = mf je (— y.\7 — r1|) eits( 2 - *)xsolryxdi(rabdrdrs. (11) 


h? lr — 7,| ‘ tee 





Now, it is well known from the theory of collisions that the F,; are 
largest for n, ] = 2,1. To obtain an estimate of the error involved 
by neglect of the bound-bound terms we therefore proceed to evalu- 
ate F,, from equation (11), using 
(sin 6, cos ¢, 
Xa1(72) = (324)—@/2)7,e—© 24 cos 0, 


sin 6, sin @2. 


As only the zero-order and second-order harmonics in the expansion 
of F,, are required, the calculation is not difficult. Use is made of 
the expansion 


; r+ Ts 
({% — %| eT! = (9, 7.) IZ (on + Ps ( , ) mm, 
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where 
tn = Kassjalyrs) lL n41/2(Yr2) 3 1 >72. 


K nst/2, 1 ns1/2 are Bessel functions which can be obtained in closed form 
so the integration may be carried out analytically except for integrals 


Tr ~ 

Cos 
_. are=r—"e , 

sin 


which must be evaluated numerically. The result of the calculation 
shows that the effect of the bound-bound terms is not greater than 


of the form 


1 per cent of the total. 


THE BOUND-FREE AND FREE-BOUND TERMS 


Although the method of successive approximation is satisfactory 
in evaluating F,, where nl # 1, 0, it is well known in the theory of 
collisions that it cannot be relied upon to estimate either the second 
approximation to F,, (the effect of distortion of the plane wave by 
the atomic field) or the effect of exchange. This is primarily due to 
the fact that, with n,/ # 1,0 the magnitude of F,,;, is determined by 
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the “‘transition”’ potentials which are nondiagonal matrix elements 
of the interaction. The correction to F,, involves the usually much 
larger diagonal matrix element. To avoid this difficulty, however, 
one may proceed by ignoring all terms in WV, which do not involve 
waves in either 7, or r,, of the same wave length as the ejected elec- 


tron. That is to say, we write 
Ve(n, r,) = 2" BS esi r>) Fs(¥:) + X10(7 (7) Fy0(12) } 


and substitute in equation (5) directly. As has been shown by Morse 
and Allis,"3 this leads to the equation 


(V2 + R v) Fi o(t1) = Ux10(11) (12) 
where 
Sxr?me? (/1 I 
h? a 
and 
81? me? _ i. , ; 
aiias ial 24 x10(r2)Fio(%2) — + (Rk? — Ex) fx:0(12) Fro(12) dr, 
“ "12 


The term v is the potential of the undisturbed atomic field and w is 
the additional term due to exchange. 
In our case F;,(7;) is alone of interest and for it we may write the 


equation in the form 


Po 2 
(“-2- oth) I Fig(ts) = UrX10("1) , CEs} 


a“ 
where 
ke a ee er sli th eal 
“4; = 3h? nf X0(":) Fi (r)r dr + 7; X10(7) FI (r)rdr 
. oJ! . 
Now, if the simple approximation in equation (3) is valid, we must 
show that the substitution 
sin kr, 
kr, 


v,F*(7;) = cos kr, 
is really a very good approximation to the solution of equation (13) 
i.e., we must show that, writing V for v + wx.0(7)/nF to! 


2 = 
-> I 
2 


13 Phys. Rev., 44, 269, 1933. 
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In Figure 1 the relative magnitudes of the terms involved are illus- 
trated, and it is quite clear that our simple approximation is a very 
satisfactory one. The effect of exchange is to improve the accuracy 
of the approximation as it reduces the effective atomic field. As a 
final test of the validity of equation (3), we solved equation (13) 
exactly, with exchange neglected, thus obtaining an upper limit to 
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Fic. 1.—Illustrating the relative magnitude of the terms 2/r? and V(r) which de- 


termine the form of the wave function of the ejected electron. 


the correction necessary. It was found that no error greater than 
4 per cent arises from neglect of exchange and distortion of the plane 
wave by the undisturbed atomic field (7). 

The effect of polarization of the atomic field on the function 
Fi.(r,) is much smaller still, as it arises only in the higher approxi- 
mations. Thus, substituting equation (6) for Q on the right-hand 
side of equation (g), we find in the equation for F,,(7,) an additional 
effective potential 


> aS nl\ T; ) Fart r ) ; 
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; ‘(1 I 
J ia a | (2 acc? -) X10(%2) xe) (T2)dr. . 
V12 rT; 


The largest contribution to this will again come from n, / = 2, 1, 
and it was estimated to be less than ro per cent of the exchange and 


where 


distortion corrections. 

It is of interest to notice here that if the zero-order harmonic of 
F,, had been required instead of the first order, the problem of accu- 
rately determining it would have been much more difficult. For 
then we would have / = 0, the dominant term /(/ + 1)/r’ of equa- 
tion (13) would vanish, and a large perturbation from the atomic 
field and exchange would result (possibly also polarization effects 
would be important). 

It is therefore established that our approximation in equation (3) 
for Vz is never likely to introduce an error of more than a few per 
cent at most in our calculated absorption coefficient. As we have 
given reasons to expect that the discrete wave function is also to be 
relied upon, at least in so far as it determines () in the region of im- 
portance in astrophysics, the values given for the absorption co- 
efficient can be regarded as accurate to a much higher degree than in 
any previous calculation of continuous absorption coefficients by 
atoms or ions with more than one electron. In no other cases have 
the effects of exchange, polarization, and nonseparability on the con- 
tinuous wave function been discussed, and the fact that they may 
be neglected for 7~ by no means insures that they can be neglected 


for other elements. 
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METASTABILITY OF HYDROGEN AND 
HELIUM LEVELS 


G. BREIT AND E, TELLER 


ABSTRACT 


The mean life of the 2s state of hydrogen in the absence of collisions is primarily 
determined by the simultaneous emission of two photons. Owing to this cause the 
mean life is ~} sec. It is expected that the mean life of 1s 28'S of He 1 is of the same 
order of magnitude and is much shorter than that of 1s 2s 3S of He 1. The non-adiabatic 
collisions with electrons are found to be more important than the static effects at small 
electron densities, leading to a transition probability ~, } sec ~' for an electron temper- 
ature ~10,000° C and an electron density 30 cm~3. The effect of static electric fields on 
the 2s state of hydrogen is worked out in more detail than previously by taking into 
account the hyperfine structure. The mean life in static electric fields which are strong 
enough to make double emission of little importance and yet weak enough to give a 
mean life > 1077 sec is found to be the sum of two lives: (1) the mean life calculated 
neglecting hyperfine structure and (2) the mean life calculated neglecting radiation 
coupling. In both (1) and (2) the mixing of 2p into 2s by the Stark effect is taken 
into account. 

Radiations due to quadrupole and magnetic dipole effects are considered qualita- 
tively for 1s 2s 3S of helium and quantitatively for 2s of hydrogen. For the latter there 
exists a small radiative magnetic dipole effect on Dirac’s theory of the electron which is 
absent in the ordinary spin model. 


INTRODUCTION 

Struve, Wurm, and Henyey' consider the mean life of the metasta- 
ble level of atomic hydrogen (2s) using Bethe’s calculation? of the in- 
fluence of the electric field. At the very low pressures of interstellar 
space the mean life exceeds the values obtainable in the laboratory 
by such a large factor that it becomes necessary to consider some 
additional effects. There are in the first place dynamic effects of elec- 
trons. The collision cross-section for the transfer from the metastable 
level to the 2p,,. level gives a transition probability of the order of 
5 X 10-5 sec“ at an absolute electron temperature of 10,000°. Here 
n is the number of electrons per cm’. The possible importance of 
this process was mentioned to us first by Professor E. Wigner. The 
static effect of the electric field is smaller at the supposed values of 
T and n. It is considered somewhat in detail on account of the in- 

' Proc. Nat. Acad., 25, 67, 1939; see also O. Struve and K. Wurm, Ap. J., 88, 84, 
1938, for discussion of electron collisions in He, and O. Struve, Proc. Nat. Acad., 25, 36, 
1939, for electron temperature. 

2 Handbuch d. Phys. (Geiger-Scheel), 24, Part 1, 452, 1933. 
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teresting physical problem presented by the overlap of the hyper- 
fine-structure pattern by the radiation width of the 2p,,, level. 

Second, the simultaneous emission of two photons? makes it im- 
possible for the 2s level to have a life longer than ~,') sec even in 
the absence of effects of collisions and of electric fields. Double emis- 
sion is presumably also important for the 1s 2s'S state of He 1 in 
which it may be expected to be of the same order of magnitude as for 
hydrogen. For 1s 2s 3S of He 1, however, selection rules make double 
emission less probable. We understand from Professor Struve that 
there is some astrophysical evidence for a shorter mean life of 1s 28'S 
as compared with that of 1s 2s 3S. 

I. EFFECTS OF ELECTRONS 
A, THE EFFECT OF THE STATIC ELECTRIC FIELD 

It will be found later that the dynamic effects are more important 
than the static ones. In the interest of system in presentation, static 
effects will be considered first and the effect of hyperfine structure 
will be worked out. From a hasty examination of the problem it 
would appear that the effect of electric fields may be quite different 
with and without hyperfine structure (hfs) because the accidental 
degeneracy of the 2s,/2, 2pi/. states is removed by the nuclear mo- 
ment. This is not so, however, because the width of the 2p,,. level 
due to radiation broadening covers appreciably the hfs pattern. 

Even though the average electric field assumed to exist statisti- 
cally in interstellar space causes a practically negligible distortion 
of the hyperfine-structure pattern, its effect on the average mean life 
turns out to be relatively important. For small fields the mean life 
for a given value of the fine-structure magnetic quantum number will 
be seen to be the sum of two parts. The first of these is the mean life 
due to the admixture of the 2p,,, state caused by the electric field and 
obtainable by neglecting the effect of coupling between stationary 
levels due to their natural width. The second contribution is due to 
the radiation breadth of the 2p,/. level which diminishes the ad- 
mixture of 2p,,. to 2s and increases the mean life. This increase ap- 
pears in the result as an additive term, provided the electric field is 
weak. 

3 Maria Goeppert, Naturwiss., 17, 932, 1929; Maria Goeppert Mayer, Ann. d. Phys., 


9, 273, 1931. 
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Assuming the proton magnetic moment to be 2.8 nuclear magne- 
tons, the hyperfine structure of the 2s,,/2., 2p:/2 levels is shown in 
Figure 1. The hyperfine-structure quantum number is denoted by f 
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Fic. 1.—Hyperfine-structure pattern of the 2s, 2p:/2 group of levels. The dashed line 
represents position of the group in the absence of hyperfine structure. Frequency 


unit »/¢ = 0.00049 cm™. 


in the figure and below. It is equal to the angular momentum of 
electron+proton in units #. The s:/2, pry. levels are split by the in- 
teraction with the nucleus into components having a separation of 
0.006 cm~ for the first and 0.002 cm~ for the second. 

The average electric field used is 


C=,, 


:( 3 \¥s ws. 
2\4rn 3.22 | 


where x = number of electrons per cm*. The Stark effect splitting 


(1) 


II 


caused by this field in the absence of hyperfine structure is 


,, (an\? e? 
AN3 lL} 2a 


2Aw _{aa\? 
7-5 X 108 (2) ow. (. 


i) 
— 


= 
II 
= 
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Here e?/ 2aq expressed in wave numbers is 1.09 X 10° cm~ and ag 
is the Bohr radius. For = 30 cm~ one has / = 0.10 cm and 
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2Aw/hc = 2 X 10°? cm. This is very much smaller than the sepa- 
rations between hyperfine-structure levels. 

The mean life of the 2p level of hydrogen is 1.6 & 107° sec corre- 
sponding to a half-value breadth of 10° sec~' in frequency units or 
0.003 cm~' in wave numbers. This breadth is sufficient to cover the 
hyperfine components of the 2s level by the natural width of the 
2P1/2 level. 

The mean life under the foregoing conditions can be treated by 
the method of the damping matrix.4 This method considers the emis- 
sion and absorption from groups of atomic levels separated by energy 
intervals of the order of magnitude of the radiation width. The pho- 
tons emitted when the atom leaves the level b’ may be reabsorbed 
so that a neighboring level b’’ becomes populated. An exponential 
decay of 6’ occurs only in conjunction with an exponential decay of 
b” having the same time constant. Definite mean lives do not be- 
long, therefore, to 6’ or b’’ but rather to linear combinations of the 
atomic states. The calculations concern themselves with finding the 
linear combinations and the mean lives. 

The probability amplitudes of the atomic states b’ are denoted by 
c(b’). For solutions corresponding to definite mean lives 


c(b’) = $°(b)A(b’, b)a(b) . (3) 


A definite mean life corresponds to any given 6. The state corre- 
sponding to it is obtained by setting a(6) = 1 and a(6’) = o for 


b’ # b. The probability amplitudes for such a state are c(b’) = 
A(b’, 6). Their values can be determined by solving 


A(b’, b)[w(b’) — w(b) + ty(b)] = i> (b')y(0’, b’’)A(b", b) , (3.1) 


where w(b’) = 27v(b’) and » stands for the transition frequency to 
the ground level corresponding to the level b’. The y(b’, 6’’) are ma- 
trix elements of the damping matrix and can be computed by 


7(b’, b”) = yoD_(a’)[x(b’, a’) x(a’, 6”) + y(b’, a’) v(a’, b’’) 


Ce) 
) + 2(b’, a’)z(a’, b”’)] | 
anc 
_ 32m4e7y3 
- 3he3 


1G. Breit, Rev. Mod. Physics, 5, 91, 1933 (see pp. 118-125). 
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the summation being performed over the lower levels a’ located out- 
side the group of interacting levels b’. Direct transitions between 
b’ due to single photon emissions and absorptions are here neglected. 
The complex numbers w(b) — iy(b) are determined by the require- 
ment of compatibility of the set of linear equations for A(b’, d) 
similarly to the conditions for real characteristic values. The ratios 
of the A(b’, 6) for fixed 6 and variable b’ determine that linear com- 
bination of atomic states which has a time dependence 


exp [—iw(b)t — y(b)t] 
for which the mean life is 


(6) = —~. (3-4) 


The general theory will now be applied to the special case of the 
28, 2p group of hydrogen levels split by hyperfine structure and in a 
weak electric field. The magnetic quantum number of the hyperfine- 
structure level (in the absence of an electric field) will be called yp, 
and the stationary states corresponding to given f and yw will be 
called y.;, and ,,, for sand p states, respectively. In an electric field 
the stationary states of the atomic system are 


Ysp = C(S, S)Wesfu + c(s, P) f'n 


(4) 


Pou = C(P, S)\Wstu + C(P, P)Wos'n » 


where c(s, s), c(s, p), etc., are coefficients of a unitary transformation. 
Since the electric field is weak, c(s, s) and c(p, p) are nearly unity 
while c(s, p) and c(p, s) are nearly zero. The linear combinations 
in equations (4) contain two wave functions only, because the 
electric field directed along the z axis is a matrix diagonal in uy. 
There are four sets of equations (4) corresponding to (f, uw) = 
(1, 1), (1, 0), (1, —1), (0, o). The matrix elements y(d’, 6’’) must 
be referred to the states ¢ since these are the stationary states of 
the atom. The selection rules of this matrix are simpler when re- 
ferred to the y. In this representation the nonvanishing elements of 
y belong to a row and column with the same y,,,, and all these di- 
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agonal elements are equal to each other. The stationary states repre- 


sented by equation (4) are 


Sic = A,S, + B,P,, S_.« = A.S..+B.P..; 
Pig = —B,S, + A,P,, Pies -~B.S..+ AP... ; 

Soe = A’S, + B’PS, So’e = A”SS + BPS ;sx sad 
Pig = —B’S, + A’P2, Pg = —B"S2 + A”'Pe. 


Here 
S: = Vs1 , So = Wato ; Sar = Ws, So = Wsvo 


and similarly for P;, etc. The elements of y are zero except in the 
squares determined by (S,¢, Pig), (S-.¢, P-1e), (Soe, Poe), and (Sy, 
P3¢). The wave functions y, ¢ will be chosen so that all A, B in the 
foregoing equations are real. The nonvanishing elements of y are 


y(Se, Se) = 2yB? ; v(Se, Pe) = y(Pe, Se) = 2vyAB ; 
(4.2) 


where indices and primes are omitted. The constant y is related to 
the mean life of the 2p level in the absence of the electric field by 


Toa = ; (4.3) 
1 


In a strong electric field A? = B? = 5, and the mean life of all levels 


is }y. For intermediate fields equation (3.1) gives 


w(b) — iy(b) = —iy + 4[w(Se) + w(Pe)] 
+ {[$w(S¢e) — 4w(Pe) + iy}? — 2tyB[w(S¢e) — w(Pe)]}! 


The w(S¢), w(Pe¢), B*? can be calculated by standard methods, and 
hence 7(6) is determined by equation (3.4). For weak fields B varies 
linearly with ©. Including effects of order © only, it suffices to ex- 
pand the foregoing formula including only the term in B*. The minus 
sign of the radical gives y(6) = 2y+ term in B’. For vanishing € this 
solution has the mean life r,,, and the correction due to the term in 
&? is negligible for the small © under consideration. The A (0’, 6) cor- 
responding to this root is small for the coefficient of ¢, and nearly 
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1 for gp. This solution corresponds to the atom’s being primarily in 
the p state during the emission. The plus sign of the radical gives 
y(b) = o for & = o. This root gives A(b’, 6) equal to unity for ¢, 
and zero for g». If © # o, the ¢, and hence also the y, are present 
only to a small extent. Calculating y(b) by expanding equation (5), 
one obtains for small € 


oe _ _27B'w(Se) — o(Pe)) , 

¥(s) = 7(b) = [lw(Se) — w(PoP + 47°’ (5.1) 
so that by equation (3.4) 

T2s = TapB {1 + 377 [w(Se) — w(Pe)]} . (5.2) 


If 7?,[w(Se) — w(Pe)]? > 1, the second term in the braces is negli- 
gible, and then 7,, = 7.)B~?. This is just the result which would be 
obtained neglecting the overlap of the s level by the natural width 
of the p level and attributing the mean life to the presence of By, 
in the wave function. The fractional error in 7., due to the neglect 
of the second term in braces becomes negligible if [w(S¢) — w(Pe)|? 
is increased sufficiently. On the other hand, the absolute error due 
to this term is independent of [w(S¢) — w(Pe)]’, because B’lw(Se) — 
w(P¢)]? does not depend on [w(S¢) — w(Pe)]? for small ©. For small 
[w(Se) — w(Pe)|? this term gives the main effect, since then 7,,B-? 
is negligible. The mean life for small © may be thus thought of as 


rt = 7 (electric field mixing) + 7 (radiation coupling), (5.3) 


where 7 (electric field mixing) = 7,,B-? and may be computed by 
neglecting radiation coupling but taking into account the hyperfine 
structure; 7 (radiation coupling) is, on the other hand, independent 
of hyperfine structure and may be computed assuming the hfs split- 
ting to be zero. 

Calculation shows that the four B’s of equation (4.1) are given by 


B? = 3C’e’aj,h[w(Se) — w(Pe)], (5.4) 


where dy, is the Bohr radius. In this formula the frequency difference 
belonging to the pair of interacting levels with which B occurs in 
equation (4.1) is supposed to be used. The difference between 
w(S¢) — w(Pe) and w(S) — w(P) is negligible on account of the 
weakness of the electric field. The values of w(S) — w(P) are 27 + 2p, 
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for (s, f = 1, = +1;p,f = 1, = +1) and for (s,f =1, 4 = —1; 
p,f =1,u#= —1);itis 2m - 6yfor(s,f = 1, un =0;p,f =0, uw = 0) 
and —2m + 10y for (s, f = 0, uw = 0;p,f =1, u=0). The first two 
levels give short mean lives, and the last two longer ones. Substitut- 
ing equation (5.4) into equation (5.2), we find 





h?[w(S) — w(P)}? hey 
— 12yG7e7a}, ” 36772 7 (5-5) 
Here 
YY = 1.57 X 1o*sec“', 
and 








 ~ ry ; oe 4304'3 sec? 
ll = s21-ti— = 430n4'3 se 
h? l 2agh ‘ 


where equation (1) is used, and 


e? 





= 3 X 10 X 1.10 X 10° sec '. 
2anh . 


Substituting these numbers, we find the two short mean lives to be 
T5 = (1.2 + 3.6) XK 105n-4/3 sec = 4.8 XK 105n-4/3 sec —_ (5.6) 


2s 


and the two long lives to be 


T35 = (11 + 3.6) X 1050-43 sec = 15 & 105n~4/3 sec 


and (=. 


T's = (31 + 3.6) XK 105n-4/3 sec = 35 X 105n-4/3 sec. 


wal 
~ 


The first number in parentheses is in each case the contribution due 
to 7.,B-?. The mean life is seen to be longer than without taking 
into account the hyperfine structure on the average in the ratio 
15:3.6 = 4. 
B, DYNAMIC EFFECTS 

The foregoing discussion is incomplete inasmuch as it replaces the 
actual field by an average value and neglects the time dependence 
of the field. If the electrons were at rest and the hydrogen atoms 
were drifting very slowly through the electron atmosphere, the static 
formulas (5) and (5.5) could be applied directly. Atoms located in 
regions of small © would then have a longer life than those in regions 
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of large ©. If the number of 2s atoms produced per unit volume per 
second be /, the equilibrium equation in such a static model is 


p = N(a2s)[r28(€)]7 
where (2s) is the local density of 2s atoms. The average density of 
2s atoms is then 


Spra(C)dxdydz__ ig 
Sdxdydz — = pralG) ; 








N(2s) = 


the integrals being extended over a region containing many elec- 
trons and atoms. In such a static picture the important quantity 
is the space average of 7... The statistics of the variation of € in 
space have been worked out by Holtsmark and by Verweij.5 Accord- 
ing to Holtsmark the chance of € being in €, © + d€ behaves as 
&-S/*d€ for large € and as G7d€ for small ©. Equation (5.5) gives 
therefore convergent results for r(2s), and the difference between 
the more accurate equation (5) and the approximation (5.5) cannot 
be very important. The most probable electric intensity is, accord- 
ing to Holtsmark, ~1.5(4.2)?/3n?/3e~4n?/%e and corresponds crudely 
to Bethe’s consideration which gives equation (1) as used above. 
At a temperature of 10,000° the velocity of electrons is ~10° 
cm/sec, and for nm = 30 cm~ the distance ] ~ 0.1 cm. The time 
during which this distance is covered by the electron is ~10~° sec. 
Attributing roughly the most probable field to the influence of 
one electron at the distance /, the order of magnitude of the 
time during which this field lasts is, therefore, ~10~? sec which 
is of the order of magnitude of r., and is somewhat smaller than 
the time |»,,,|~' corresponding to the transition frequencies between 
the hyperfine-structure components of 2s and 2p. The latter times 
range from 3 X 10~* sec for the interval 2v, to 6 X 107° sec for the 
interval 1ov,. The effect of admixing the 2p level by the quadratic 
Stark effect cannot be fully developed, therefore, at this T and n, 
although it should become more important at smaller T and 1. 
At larger densities and temperatures the most probable field lasts 
only a fraction of 7,, without changing its value appreciably. The 
5 J. Holtsmark, Ann. d. Phys., 58, 577, 1919, and Verweij, Pub. Astron. Inst. Amster- 


dam, No. 5, 1936. 
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static effects should be noticeably decreased in this case. In fact, 
for an electric field alternating sinusoidally with the time, the transi- 
tion probability, neglecting hyperfine structure, is proportional to 


where v is the frequency of the field. For high electron densities 
and temperatures (47v7,,)? > 1 and the equations (5) cannot be 
applied. There is perhaps an offsetting influence in the statistical 
character of the time dependence of © which occasionally permits 
the field to last at ~e//? for relatively long times. These effects are, 
however, difficult to estimate. It is doubtful that they can affect 
appreciably the decrease in the transition probability due to the 
alternating character of © which has just been discussed. 

In addition to effects due to electrons at a distance ~/ which vary 
with 2*/3 there are also definitely non-adiabatic effects due to col- 
lisions with electrons. These effects give a transition probability 
proportional to ” and will be discussed next. Their possible impor- 
tance has been mentioned to us in this connection by E. Wigner. 
For He the effect of electron collisions has been considered by Struve 
and Wurm.’ In hydrogen this effect is relatively more important 
because of the proximity of 2s to 2p. Even a comparatively distant 
collision transfers the atom from the 2s into the 2p state with an 
appreciable probability. From this state it radiates with the transi- 
tion probability (7,,)~' either after or during the collision. 

The process is very similar to that considered in the theory of the 
loss of energy of alpha particles in going through matter, which has 
been treated by Bohr, Mott, Williams, Bloch, Bethe, and others.° 
The present problem is characterized by the extreme smallness of 
the excitation energy and by the fact that the incident particle has 
a velocity comparable to that of the atomic electron. The impact 
parameter method can be expected according to Mott and Williams 
to lead to much the same result as the Born method and will be 

®N. Bohr, Phil. Mag., 25, 10, 1913; N. F. Mott and H. S. W. Massey, The Theory of 
Atomic Collisions, chaps. xi, xiii, Oxford University Press, 1933; N. F. Mott, Proc. 


Camb. Phil. Soc., 27, 553, 1931; F. Bloch, Zs. f. Phys., 81, 363, 1933, and Ann. d. Phys., 
, 285, 1933; H. Bethe, Ann. d. Phys., 5, 325, 1930; E. J. Williams, Proc. R. Soc., 139, 


163, 1933- 
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used here since the calculation can be arranged to follow Bohr’s 
classical treatment rather closely. The influence of the electronic 
and nuclear spin will be neglected at first but will be taken into ac- 
count later on. The probability that after a collision the atom has 
been transferred into the p state is then 


le PR sie iw td ‘ sie iwt a 
ie | tep| a ,e at a mz Gye at = C ; 


where G,, G, are the components of & produced at the atom by the 
colliding particle, taken perpendicular and parallel to the velocity. 
The frequency w/27 is the transition frequency between the s and p 
levels while x,, is the matrix element of the x co-ordinate of the 
atomic electron between the s level and the p sublevel directed 
along x. The electron velocity will be called v, and the distance from 
the atom to the electron path will be called p. One finds then 
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C,e“"dt = — —_—_——- dx = 2f | — 
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where the functions f, g are those introduced by Bohr. The collision 


and 


cross-section for excitation is 
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Pp '*) 


where 


P(z) = f(z) + g7(z), 


again in Bohr’s notation. The lower limit of integration will be 
taken as 
h 


mv 
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This choice is somewhat arbitrary, but it is natural since for it 
C ~|2.,/ay|? so that the approximations used in deriving C cease 
to apply for smaller p. The results are not very sensitive to the 
choice for po, since p, enters logarithmically. The effect of collisions 
with p < pp, will be neglected. Even though these collisions should 
be considered for the stopping-power of an atom on account of the 
large energy transfer, they are not of special importance in the 
present problem, since only the very small energy transfers are of 
interest. Proceeding as in Bohr’s paper and using his evaluation of 
his constant k = 1.12, one obtains 


874! Xsp|? I.I mv? 
¢= —— J | —— R (S) 
hv? hw 











The quantity in the logarithm cannot be considered as certain on 
account of the arbitrariness in the choice of po. The factor multi- 
plying the logarithm is just twice that obtained for the excitation 
of an atom in the normal state by a high-velocity electron. This 
formula is 


re4| x.,|? 2mv? 
= 4ret| ep” ln (2) , (F ) 


where AE = excitation energy. For a fast electron the wave length 
is < dy, and the integration should be limited at ~a,y since other- 
wise the electron penetrates into the atom. In the foregoing calcula- 
tion modified for fast electrons the quantity under the logarithm 
is, therefore, essentially v/(waq) = vme?/(wh?); for wh ~ e/ay 
this is 


mvady vk (meee 1/2 (me 
—_— == ——_- = ——__--- ~ - ~ 
h e? e? AE 


gy mv? 
in | —) ~ 3iln|— }. 
Way AE 


In this way the use of ay for the lower cutoff together with the rela- 
tion wh ~ e?/ay changes the factor 8 in formula (S$) to the factor 





so that 


4 in formula (F). In the present problem wh < e?/ay, and 
h/myv ~ 10° °T~? cm. For T < 10,000°, h/mv > 10-* cm, so that 


—————E— 
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the condition of the incident electron staying out of the atom is 
satisfied partially. At temperatures above 10,000° formula (3S) 
becomes poor because the velocity of the incident electron is of the 
order of magnitude of the velocity of the atomic electron, and p, 
penetrates into the atom. 

A systematic description of the Born method is given by Mott and 
Massey in The Theory of Atomic Collisions (Oxford University Press). 
On pages 176 and 177 the relation to the present results can be seen. 
The factor 5 is introduced essentially in their inequality (38) which 
limits the momentum change to [2m| E,| /h’]'/? = Ko. If, however, 
the calculation leading to their equation (40) is made for a case 
in which one may set e’“* = 1 + iKx for the computation of 
feX~~.Wxdr, then one obtains from equation 13, page 161, essen- 
tially twice the value of Q., as given by equation (40). From this 
point of view the essential condition for preferring equation (S) to 
equation (F) is that mv < mc/137 or v < 2 X 10° cm/sec. This is 
also the result of the consideration with the impact parameter 
method. 

Neither equation (S$) nor equation (/) is exactly applicable. For 
a more exact discussion one would have to evaluate few pdr. 
For present applications formula (S) is probably good enough since 
the primary question is one of order of magnitude. In view of the 
fact that the impact parameter method and the Born method indi- 
cate that the true result is between formulas (S) and (F), the former 
will be used, and the results will be considered in the nature of an 
overestimate for the transition probability. 

On account of the splitting of the 2p level into 2p,,2, 2p3/2 and the 
further splitting of these levels into hyperfine-structure components, 
the calculation leading to formula (S) must be modified. The pro- 
cedure is straightforward and leads to 
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The last form is convenient for numerical substitution, while the 
first is intended to show the relation of formula (S’) to formula (S). 
The | x5, |? would be replaced by 347; if there were no hyperfine struc- 
ture. The coefficients }, {, } are the relative strengths of transitions 
from f = 1 tof = 1, from f =o to f = 1, and from f = 1 to f =o. 
The energy denominators 2/w,, 1o/w,, Ow, are the corresponding 
energy differences between s and p components. The frequency 
w,/ 2m corresponds to »,/¢ = 0.00049 cm of Figure 1. 

For transitions from 2s,,. to 2p;,. by neglecting hyperfine structure 
one obtains 


nh? mv? 
(rcon)2p,;, = 48x —— In (- 
coll/2P; /2 4 vm? hw(sp, = 


4 
= “—— a -nln(5.7T) . 


(S”’) 


. . ry. re) 
For n = 30 cm-3, T = 10,000 one has 
I I 
—I a ad — 
(Toon) 2p, /2 a ne a eZ bd (Toon) 2p, 2 = 7 » * 
/ 1400 sec 1000 sec 

and 

Toon = (toon)2 + (to) = 

coll coll *P1/2 coll *P3 2 600 sec - 


The value of the static effects corresponding to the foregoing n and 
T is tr ~ 1/5000 sec. For higher n the formulas (5) and (S$) indi- 
cate an increased relative importance of the static effects. As has 
already been discussed in the beginning of this section, however, the 
static effects cease to have a definite meaning at higher electron den- 
sities, because at the smaller average distance of the electron from 
the atom, the electron distance changes by a large fraction of its 
value in a time small compared with the mean life of the 2p level. 

It is clear that the foregoing estimates could be improved in many 
ways regarding the chance of transfer in a single collision, the effect 
of radiation before the transfer is completely effected, as well as 
with respect to the statistical nature of the field and the velocity 
distribution of electrons. Since at the density and temperature of 
electrons supposed to exist in interstellar space the effects of elec- 
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trons appear to be of secondary importance, such improved estimates 
do not appear to be essential at the present time. 
II. SIMULTANEOUS EMISSION OF TWO PHOTONS 
The chance per second for the simultaneous emission of two pho- 
tons with one photon in the frequency range dv’ is 
(TU!) n'n”(TUT nn 
Van tv” 


(ru) nn (TU) nn | 
Van + v’ 





fo hh 
10247 esp 3y"'3 x 
Af SS 
h?c® — 
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) aor | 


Here v” is the frequency of the second photon emitted simultaneous- 


a 








ly with v’. Conservation of energy holds, so that 


The upper level (2s) is designated by n, the lower (1s) by n’. The 
intermediate levels are denoted by nm”. They are in this case the 
p levels. Matrix elements are denoted by suffixes such as (1u’) wn”. 
The polarization vectors u’, u’’ have absolute value 1 and are paral- 
lel to the electric intensity of the photons Jw’, hy’’. The displace- 
ment vector of the electron is denoted by r. The averaging in the 
foregoing formula is supposed to be made over the directions of 
propagation and over the polarization directions independently for 
hv’ and hy’. The formula is supposed to be used for the probability 
of simultaneous emission of a pair of photons, and there is no dis- 
tinction with respect to the order of emission. The total probability 


~=A,= bf Ade’. (6.1) 


Equation (6) applies only if the electric dipole transitions to in- 


of emission is thus 


termediate states predominate over other ways of obtaining these 
transitions, such as quadrupole and magnetic dipole effects. It ap- 
plies in the case of the 2s state of hydrogen on account of the pre- 
dominance of dipole transitions to the p levels. In applying equa- 
tion (6), one may neglect the effects of hyperfine structure and of fine 
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structure since the operator r is diagonal in the nuclear and electron- 
ic spin. The s levels are then single, and every p level is triply de- 
generate. It is convenient to use for each energy three linearly in- 
dependent p states having the transformation properties of the Car- 
tesian co-ordinates x, y, and zs. These eigenfunctions will be denoted 
by &, », and ¢. The sum over n”’ in equation (6) will be taken first 
over &, y and ¢ for each energy E,”. One has 


~ y 
SN (ru’) nin (TU) n''n _ (xuz)n'¢(XUy int (yaty)n’n (Vey on 
nt + (2u;)n’s(stz’)en 


= Sn'pen(uu, + uyuy’ + ujuy’) 


because 


Xn'tXin” = Va'nVan"” = Sn't2tn" 5 Vali = Vet an FD 
by symmetry. 
The sum over n” is, therefore, 


“1 
“nme-mn + 


Vnn + v" 


wa’) > i — bees 
rome Yan TY 
Here the principal quantum number of the p state is denoted by m 
and the matrix elements z,,,, are between the state 7’ and the state 
of type ¢ with principal quantum number m. The factor (u’u’’) 
shows that there is a statistical relationship between the directions 
of polarization of the two photons and that the emission probability 
is proportional to the square of the cosine of the angle between the 
polarization vectors. Such a relationship is natural, since angular 
momentum must be conserved. The average over the polarization 
directions for (u’u’’)? is 3. Substituting into equation (6) one has 
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3hc° pam Vnn + V Vmn + V 
m=2 
The summation over the principal quantum number must be extend- 
ed also into the continuous spectrum. The integral over the con- 
tinuous spectrum is not indicated in the foregoing formula but is 


understood to be present. 
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In the evaluation of the sum one is helped by the summary of 
values of the squares of radial integrals given by Bethe’ for the 
discrete spectrum and by the closed formulas of Gordon*® and Stobbe? 
for the continuum. By using the Bohr radius as the unit of length, 
the wave functions will be used as follows: 


1a 


™~ 


; ae E \ 1/2 
Ys = (47) ’ Wap _ (5) Rew ’ 


4r 
f rPRdr=1, 


: r 
Ry = 247° re"? , Ry = (* )orelir (1 — 3) 
6 
The notation 
Rip = if PRisRmpdr (6.4) 
0 


will be used. One has 








(m > 1) Rup = — (6.5) 


(m _ ce +5 


28(m — 1)?™ | 2 





(m > 2) Rip = 





2'7m7(m? — 1)(m — 2)?™-6|1/2 
Wal (m a 3)" ’ (6.6) 


and 
Roa» = —3V3; (6.7) 


where positive values of the square roots are understood. The radial 


integrals in Stobbe’s notation are 


2/x) arc tan 2 x 2 








, ; e€ d 
C™ 1 = Cy = 4? — ) 6.8 
0,0 I 4 (1 —€ 25 da 2 (1 ote x?)s 2 ( ) 
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a ' ae (2/x) are tan2z 1 (1 + =) 2 
C 6. = ( 28 a cals »>—2m/xr\1/2 a -2)\3 (6.9) 
(1 —¢ ) (1 + 4x7)3 
7 Op. cit., Table 15, p. 442. § Ann. d. Phys., 2, 1031, 1929. 


9 Ann. d. Phys., 7, 661, 1930 (see Tables 2 and 3). 
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where x = kag, dy = h?/me?, k = (2mE)*/?/h. The superscript m, is 
the radial quantum number for the continuum. To obtain the con- 
tribution to }onRnpRmpy due to the continuum with ay as unit of 
length one takes 

Cal le. 


The quantity x used here is Stobbe’s x for the ts state. 
In terms of x the energy in the continuum is 


: : e? 
E=x—. 
20H 
Expressing the energy in Rydberg units e?/2ay, the energy is 2’. 
The smooth joining of the discrete and the continuous is verified, 
giving at the series limit 27e~4 as the contribution per unit energy 
range with ¢?/2ayz as the unit of energy to }om(Rmp)? and 2™e-* to 
28 \ 

>. (Rin) 
The relative signs of the radial integrals are important for the 
evaluation of equation (6.2). These can be checked by working out 


in terms of them 


(9*).5, 2 = f r1R,.R.dr = —31%V 2 = —2.098, 


which should also be expressible as 


(7?) 1s, ae pe se J Crs ht : 


Using Bethe’s table on his page 442 and estimating the contribution 


~ : ‘ 
from m = 8 to m = © aso.19, one obtains > RupRmp = 2.38 and 
geome t ! . 
m= 3 
R;>R3> = —6.70. The integral over the continuum is evaluated us- 
ing Stobbe’s numerical tables? of C,,, C2, as +1.33. This gives 
(r?):s)28 = 2.38 — 6.70 + 1.33 = —2.99. The good check with the 


directly computed value indicates that no important slips have been 
made in the values of the matrix elements. 
The radial integrals can be used also for the evaluation of the 
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products Zn'mZmn = 3RmpRmp occurring in equation (6.2). The most 
significant contribution comes from m = 2. This contribution is 
negative, while the other contributions are positive. The sum, in- 
cluding the continuous spectrum, is negative. It will be convenient, 


therefore, to consider 


co C 
~ P : : 
: Te a. ee ee rs 
re) ( ad + iia (— + py’ Vnn y" 


m=2 


The quantity C is positive, and the contribution to C due to m = 2 
is positive, while the other contributions are negative. Expressing 
radial integrals in units dg one has 
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The quantities in parentheses have a minimum for y = 3 and in- 
crease steadily toward y = o and y = 1. Keeping the term for 
m = 2 as it stands and letting y = 3 and o in the other terms, one 
obtains, therefore, as upper and lower limits for C, 


GscCs<c,, 
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and 
m | 
* 
b, a . R*. pall cS | poe < Jae 
— PN — ms 1 — Mm 
m=3 
+ (cul a roe 
8 238 -— a ‘a ace eaptic” rf x . 
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Substituting upper and lower limits of C into equation (6.2) and | 


then into equation (6.1), one obtains as upper and lower limits 


I 34 (< °/ ¢ a? aidi b? 
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6 15 140 


the upper limit corresponding to a,, 6, and the lower to a., b,. The 
number a, = a, can be calculated exactly and is (4/3)6.70 = 8.93. 
The evaluation of 0,, 6. has been made only approximately by using 
the radial integrals up to m = 5, neglecting 1/m? for the remaining 
contributions of the discrete spectrum, and using « = 0.4 in the 
denominators of the integrals over the continuum. This gave }, = 


12.2, 6, = 22.2, and 


I ) ' 
4.4 sec' <—~ < 8.7 sec. j 


The values of the denominators involving frequencies used for 0, 
correspond to the maximum of the statistical weight factor vv} at 
/ 


v’ = 3¢?/16ayh. It is, therefore, likely that the upper limit for 1/7 
is closer to 1/7 than is the lower limit. It is probable, therefore, that 


= I , 
Osh SOC 5. =< Bec 


III. DOUBLE EMISSION FOR HELIUM AND OTHER EFFECTS 


In He the (1s 2s) 'S, and (1s 2s) 3S, levels are also metastable. 
Direct transitions with emission of a single photon from (1s 2s) 'S, to 
the ground level (1s)? 'S are strictly forbidden for dipole or multipole 
radiation of either the electric or the magnetic type. Nevertheless, 
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double emission from (1s 2s) 'S to the ground level should take 
place with roughly the same probability as from the 2s state in hy- 
drogen. For the °S state, however, the spontaneous emission of two 
photons is very much smaller, since the intermediate states either 
combine poorly with *S or with 'S. The factor due to this is of 
the order of magnitude of the admixture of *P, in 'P,, which is 
2000 cm~' ~ 10-3 for matrix elements or ~10~° for the in- 


r em 


tensities. Double emission for *S should be, therefore, much less 
probable than for 'S. 

While single photon emissions are strictly forbidden from (1s, 2s) 
‘S, so that collisions and double emissions are the only effects, 
there are minor radiation effects in *S, which compete with double 
emission. Without any attempt at precision the general order of 
expected magnitude of these effects will be considered. 

The ordinary spin-orbit interaction is of the form 


A,o, + Ao, , 


where ¢,, ©, are the Pauli spin vectors for the two electrons, and A,, 
A, are orbital vectors. According to standard-selection rules, this 
interaction mixes 'S, only with *P,, and 4S, is mixed only with 'P, 
and 3P,. The terms under consideration are even and, therefore, the 
P terms when expanded in terms of product wave functions are also 
even. They cannot be, therefore, of the (1s, np) or (2s, np) types, 
and it suffices to consider admixtures of pp’ and dd’ configurations. 
Configurations p? contribute only *P, and *P, but no 'P,. Taking 
1 cm as the order of magnitude of the nondiagonal perturbation 
matrix element and using ~ 300,000 cm~' for the distance between 
(1s)? and (2p)? and ~150,000 cm~' for the distance between 1s 2s 
and (2p)’, the coefficient of the *P, wave function can be expected 
to be 10~ or less for (1s)? as well as for 1s 2s. The admixture of 'P, 
in (1s 2s) 3S, is not likely to be bigger and will be considered as 
being taken care of by a coefficient 1o~> in the wave function. 

In addition to the ordinary spin-orbit interaction there is present 
the spin-orbit-spin interaction of the type 7, = 3[3(re,)(re.) — 
r°(o,0,)|. It produces no mixing of 'S, since 7, applied to a singlet 
wave function gives zero. In 4S, it produces an admixture of *D, in 
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addition to 3P,. The character of the terms is thus 'S, + a *P, and 
38, + B'P, + y 3P, + 6 3D, where a, 8, y, and 6 are ~10~5, and the 
spectroscopic symbols stand for normalized wave functions. 

In the foregoing estimates the (2p)? configuration was used in esti- 
mating the order of magnitude of the coefficients. The important 
coefficients will be seen to be a for (1s)? and y. 

The presence of terms in the pp’ configurations is suggested not only by the 
selection rules but also by the approximate form of the wave function. The 
orbital function for 'S, has the form Us(r,, r2; 7:72) wherer,, 7, are, respectively, 
the vector distances from the nucleus to the electrons 1 and 2. The interaction 
energy A,o, + A.o, operating on the electronic wave function gives a linear 
combination of triplet spin functions and components of (A, — A,)U’s. The 
vector A, consists of two parts 


he? {3 ; ’ 
At = ee (7: X pi] + r-3[(r1 — 2) X (2p2 — p)]} ; 
The first part contains the angular momentum L, = [r, X p,|._ Because of it 


the contribution of A, — A, occurs only on account of the dependence of U’s 
on r; + r, and not onr, orr,. The orbital functions arising in this way are of the 
type [r, X r.|dUs/0(r.7.). The second part of A, gives rise to orbital functions 
of the type 
Fi) Fi) 20 ; 

rats X Wal) oar, + rar, t O07) | US 
If Us did not depend on r,r,, only the second part of A; would matter, and the 
admixture would be entirely of the pp’ type. This is in agreement with the fact 
that if the dependence of Us on rr, is neglected, the function U’s is entirely of 
the ss’ type so that d states of single electrons cannot be brought in by the spin 
orbit interaction. If the dependence of U’s on 7,7, is taken into account, the pp’ 
configurations are not the only ones entering into the admixture. 


Quadrupole radiation between *S, + 6 'P, +.... and 'S, + 
a 3P, is excluded by the exact selection rule which prohibits direct 
jumps between any two levels with angular momenta /’,7”’, for which 
a triangle with sides 7’, 7’’, 2 cannot be constructed. 

The magnetic dipole operator consists of a part proportional to 
the total angular momentum which does not give rise to transitions. 
The remaining part of the magnetic dipole contains o, + ¢, and is 
proportional to the total electron spin. This part gives zero when 
applied to singlets and leaves the space function unchanged when 
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applied to other terms. The multiplicity of a term is also unchanged 
by 7, + ¢,. Hence the effect of the magnetic dipole is to give a transi- 
tion only between y ?P, which is admixed to *S, and a *P, which is 
admixed to 'S,. Using the estimate 10o~ for a and y, one obtains the 
very small value 








105 X 105A |? E 
ees ~ 10% 
aHymc 


for the ratio of the transition probability to that expected for an f 
value (oscillator strength) 1 with electric dipole. This effect is so 
small that one does not need a more accurate estimate. 

In the absence of external fields and nuclear spin the 2s,/. term of 
hydrogen is pure. The electric quadrupole and magnetic dipole ef- 
fects vanish in this approximation, the first because one deals with an 
s — s jump and the second because the ts and 2s states are strictly 
orthogonal. On account of the proton spin there is a slight perturba- 
tion of s,/. and d;,.. The matrix element responsible for this does not 
vanish because the coupling of the electron to the proton involves 
Dirac’s a, which in the Pauli approximation is neither a pure orbital 
nor a pure spin vector. One may expect therefore s,/.(f = 1) to be 
weakly mixed with d,,.(f = 1), but s,,.(f = 0) cannot be affected 
by this perturbation since d,,, has no f = o. Even if the life of 
2S:/.(f = 1) should be shortened by electric quadrupole radiation 
from the 3d;,,.(f = 1), the f = o part could still show absorption 
lines. This perturbation appears to have only theoretical interest 
and should besides be small since the closest d term occurs for n = 3. 

The proton presumably does not act on the electron quite as a 
point charge and may have a finite volume. As long as its action 
may be represented by a fixed central field, the 1s and 2s functions 
are orthogonal, and the magnetic dipole effect vanishes. If, however, 
the proton field depends on the velocity of the electron, this is no 
longer true, and other types of radiation could conceivably occur if 
the proton and electron were more intimately coupled. This already 
takes place weakly for Dirac’s electron equation. 

The approximate Pauli representation of the electron spin gives 
no magnetic dipole radiation for 2s,,. of hydrogen. A slightly differ- 
ent result is obtained by using Dirac’s equation. The matrix ele- 
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ments of the components of 3[r X a] take the place of the matrix 
elements of r in the formula for the transition probability due to 
the electric dipole. The radial integral {PP (fisrs + frsgis)dr occurs 
as a common factor in the matrix elements. Here g is the “‘large”’ 
radial function’? and is approximately equal to the nonrelativistic 
one. Approximately f is a constant multiple of dg/dr. In this ap- 
proximation the foregoing radial integral vanishes as is seen by par- 
tial integration using also the approximate orthogonality of g,, and 
£25. The radial integral is not exactly zero, however, but only small. 
Its order of magnitude is (e?/fic)?(h/mc). The emission probability is 


Met : : 
31 (eV ~ 5 x i oe, 
which is negligible. The foregoing relativistic magnetic dipole effect 
takes account of all terms in the usual expansion of the factor e’*” 
in the expression for the vector potential up to 1 + 7kr. It is sensi- 
tive to the nuclear charge on account of the extra factor (Ze?/hc)* 
due to the difference between the Dirac and the Schroedinger radial 


functions. 


The writers are grateful to Drs. Struve, Henyey, and Unsold for 
stimulating discussions. 
UNIVERSITY OF WISCONSIN 
AND 
THE GEORGE WASHINGTON UNIVERSITY 


10 P. A. M. Dirac, The Principles of Quantum Mechanics, p. 252, Oxford: Clarendon 


Press, 1930. 


THE GENERALIZED THOMAS-FERMI METHOD 
AS APPLIED TO STARS 


R. E. MARSHAK AND H. A. BETHE 


ABSTRACT 

A treatment is given of stars of intermediate density such as dense red dwarfs and 
subdwarfs, using the statistical method of Thomas-Fermi to calculate the free-electron 
density, the gas pressure, and the opacities. The Thomas-Fermi equation has been 
generalized to take account of the high temperature and a mixture of elements. By 
taking finite boundary conditions as for metals, it is then possible to integrate the star 
equations. 

In this note we wish to indicate a method of integrating the equa- 
tions of stellar equilibrium for stars of moderately high central den- 
sity (of the order of 1000 gm/cm). The method could have been used 
to determine the temperature-density distribution in the interior of 
a white dwarf star in a more elegant fashion than has been done;' 
but the labor involved in doing so would have been needlessly great 
since strong degeneracy and, therefore, complete (pressure) ioniza- 
tion prevail throughout most of the star.? However, it seems that it 
would be very worth while to investigate the internal constitution of 
the dense red dwarfs and the recently discovered subdwarfs} on the 
basis of a “‘generalized”’ Thomas-Fermi model; the information thus 
derived should be more accurate than that given by any of the usual 
stellar models and would therefore throw light on the nature of the 
energy-production process of this interesting group of stars. 

The starting-point, then, is the Thomas-Fermi equation for elec- 
trons, generalized to include the effect of high temperatures and 
subject to finite boundary conditions as for electrons in metals. Solu- 
tions of this equation lead directly to expressions for the free-electron 
density, the gas pressure due to electrons, and the screening radius 
which enters the formula for the conductive opacity. Moreover, the 

« R. E. Marshak and H. A. Bethe (to be published shortly). 

2 Moreover, it would be necessary to take account of relativistic degeneracy; this 
introduces additional difficulties for the Thomas-Fermi method. In this note we 
restrict ourselves to the nonrelativistic Thomas-Fermi atom. 

3G. P. Kuiper, Ap. J., 89, 548, 19309. 
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mass density can be expressed in terms of quantities which the gen- 
eralized Thomas-Fermi equation yields, and the star equations can 
then be integrated. These considerations hold for one element but | 
can easily be extended to a mixture of elements such as the “Russell” 
mixture. 

We write for the electron density n | 


2 


ot a 2 
= (27)3 2h3 3V x (log a ): : + 


us 


(1) 





8 log? a’ 


Here 4 = Planck’s constant divided by 27, m is the mass of the 
electron, k is Boltzmann’s constant, T is the temperature, and a’ is 
related to the Fermi energy, ¢, and the potential energy in the 
neighborhood of the atom, eV, by the equation 





et ee 6 gg 2 
log a’ = LT -a(i+v)-¢. (2) 


Then equation (1) becomes 


- alin sad call 








i= 


V3/2 [ 4 sat] | 


Bel? 


3m7h3 


This electron density when inserted into the Poisson equation 
AU = 47mne leads to an equation of the Thomas-Fermi form 


(4) 





an a 
dx? g? | 


7 


4In this approximation only the leading temperature-dependent term is retained; 
if all orders were included, the expression in the brackets would be replaced by 


jr: :¥ (6 ~ 26g aA =, 


mS, £0 +4 6:2 


where £, is the Riemann-zeta function of order v (cf. Sommerfeld, Zs. f Phys., 47, 9, 
1928). Also the exchange term is neglected and justifiably so (cf. H. Jensen, Zs. f. Phys., 
IIT, 373, 1938). 
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In equation (4) we have set® 


——— 
l = - € p r= ax, (4a) 
r 
_ - or 1/3 7—(1/3) _ he ai Da 
il (=) AZ ; &ee* Bohr radius , (40) 
2 4/37710/3 k \2 mead 
si aa (x) Ry = — = Rydberg energy. (4c) 


This equation is to be integrated subject to the boundary conditions 


dn : 
(A) (0) = 1 (B) —=o0 at *+=X. 
dx 
The first condition implies that in the neighborhood of the nucleus 
the potential goes as Ze/r, the second, that the atom occupies a 
spherical volume with radius R = aX. The second condition is 
equivalent to the requirement that Z electrons be contained in the 
"R 

sphere of radius R, i.e., arf nr-dr = Z. This can be seen either 


by direct integration using n = AU /47e or from the fact that the 
electric field must vanish at the surface of a neutral atom; this means 
that dV /dx = o at the boundary and, since for a given temperature 
V is a unique function of , also dn/dx = o. Equivalent to dn/dx = 
o is the condition® d®/dx = ®/x. For a given X the boundary 
condition dn/dx = o at x = X will require a certain value c(X) for 
the derivative d®/dx at x = o. 

Once equation (4) is integrated for given temperature and mass 
density, the free-electron density can be found immediately from 


= — 3 2z* - Py : I 4 pI°x* é (5) 
or3a3X3/2 < Py 
and the electron pressure from 
c/a Frofe.s 5/2 ~QT2V2 
elec 279 3Z' 3p? at 587 2X 2 
lec mepees eee 19 fe eee 6 
PG Cy oe 3q4 xX 5/2 Py ( ) 


5 FE. Fermi, Zs. f. Phys., 48, 74, 1926. 
6 J. C. Slater and H. Krutter, Phys. Rev., 47, 559, 1935. 
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Equation (6) follows from the definition of pressure as the mean rate 

of transfer of momentum across unit surface.? A solution of equa- 

tion (4) would give the screening radius R = aX directly, and this 

expression could be used to evaluate the conductive opacity (cf. 

eqs. [33] and [41] in n. 1). Finally, the mass density p is given by 
_ 3Amz . 


= —— + A is atomic weight . ¢, 
p 47R3 : 6 ( 


sl 
— 


The above formulae hold for one element; when a mixture of ele- 
ments is present the only modification introduced is one of self- 
consistency: the gas pressure due to the electrons and, therefore, the 
free-electron density must be the same at the boundary of each 
atom.* The mass density is then given by 


eo 3.89 (7a) 
can aia Si = a — re yen = Te Taal . i¢ ? 
7 Loe. we 
47 Az @ - ann 


In equation (7a) cz is the concentration of element Z. 

We have shown how certain astrophysically significant quantities 
can be expressed in terms of solutions of equation (4). Since this equa- 
tion can be integrated only by numerical means, we point out two 
ways of doing this. The first is straightforward and involves a nu- 
merical integration by a simple polynomial approximation.’ The 
second is to treat the temperature-dependent term in equation (4) 
as a perturbation and to make use of tables’? for the unperturbed 
Thomas-Fermi function which have already been calculated. In the 
latter case, if we let 6 = & + ®,, then ®, satisfies the differential 
equation 

Cd, 30%? — BT2xs? 


eC. SER, ita 8 
de 2x, BN? ” 


7p= +f N(q) + q + vqdg, where N(g)dq denotes the number of electrons per 


unit volume with momenta between g and gq + dq, and 2q is the velocity associated with 
momentum g. The integral can be evaluated by means of Sommerfeld’s formulae (cf. 
Chandrasekhar, An Introduction to the Study of Stellar Structure, 1939). 

8 This essentially fixes the X for each atom and replaces the condition obtained from 
arguments about the free atom (cf. n. 6). 

9M. F. Manning and J. Millman, Phys. Rev., 53, 673, 1938. 

10 Professor Slater possesses many copies of these tables; we are indebted to him for 


sending us a copy. 
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with the boundary conditions 


, ra) = . Or. _ 
(A‘) &(0)=0 (B) ~¥ -(2),. 


In this derivation account has been taken of the fact that 





with 


d®, ®, 
P,(o) = tI, (=). = (2), . 


In terms of ®, and ®, we have for the free-electron density 


n = 322? (Po\ , 162? (Me 4 322°8 ( x \"" 
omaz \ x /x 3m3a3 \ x3/? Jy omsaz \@/x ’ 


and for the pressure 





f 


25/3 F10/3 v2 5/2 47/6 Z10/: 52 3/2 5/2 

p' lec _ ? en ¢ 1G 4 sae ae '¢ 5)? ®, \%/ 
oro ; J 2 8/3a-11/3 4 5/2 , 
x 1s Ga\ & X 


5.38 3mis af 





X i 
9225/3 Z10/3 eB ®, 1/2 
gins as \a)x 


Some rough calculations we have made indicate that mass densities 
of the order of a thousand fall within the range of the Slater- 
Krutter tables. The Thomas-Fermi procedure is thus quite feasible 
for the dense red dwarfs and the subdwarfs and should reveal a 
great deal about the temperature-density distribution of these stars, 
provided the hydrogen contents are not too high. 

UNIVERSITY OF ROCHESTER 


AND 
CORNELL UNIVERSITY 








PROPER MOTIONS IN THE GALACTIC CLUSTER M 67 
E. G. EBBIGHAUSEN 


ABSTRACT 

Proper motions are derived on Yerkes plates for an area 30’ in diameter centered on 
the cluster M 67. The results are contained in Table 2. Plots of the motions are given 
in Figures 1, A and 1, B, and the position of the cluster members in the cluster in 
Figures 2, A and 2, B. Figure 3 gives the spectrum-magnitude diagram of the cluster, 
and Figure 4 the luminosity function of the cluster members. 

This paper gives the results of the determinations of proper mo- 
tions in the galactic cluster M 67 (a = 8'45™8; 6 = +12°11'; 1 = 
184°; 6 = +33°). For this cluster Trumpler gives 2~3a as its type, 
meaning that the cluster has many giants and that the upper limit 
of the main sequence is at spectral type A. For M 67 he gives 740 
parsecs as its distance and 16’ as its apparent diameter. Because of 
its large distance from the galactic plane (+406 parsecs) and the 
large number of late-type giants, the cluster is of considerable in- 
terest. 

Previous work on the cluster has been done by Olsson at Stock- 
holm in 1898" by Fagerholm,? Shapley,’ and van Rhijn.* The results 
of these papers will be discussed later when the occasion demands. 


I. THE MATERIAL 
The material consisted of four pairs of plates, each pair having 
only one image of each star. Table 1 is a list of these plates. I wish 


TABLE 1 
' Interval Limiting pv 

Old Plate New Plate : Symbol 
F : ‘ in Years Magnitudes . , 

F 179 F 599 26.1 13.0 P ] 
F 174 F 600 A 26.1 Ez..0 Pi 
F 305 F 597 A 17.0 12.3 P Il 
Sch 241 F 601 A 35.0 [1.7 Pilly 


to express my thanks to Dr. J. Titus and Dr. C. Hetzler for having 
taken these plates in my absence. 

t Pub. Stockholm Obs., 6, No. 4, 1898. 3 Mt. W. Contr., No. 117, 1916. 

2 Diss. Upsala, 1906. + Groningen Pub., No. 33, 1922. 
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II. MEASUREMENT AND REDUCTION 

The procedure adopted in this paper is the same as that used in 
the author’s paper’ on the galactic cluster NGC 752, where it is fully 
described. The only difference in procedure is that the measurements 
were made only with the north edge upward and not again with the 
same edge to the right. This resulted in two measures of Ax and two 
of Ay, which were then averaged to form for each star a mean value 
Ax and Ay. A least-squares solution using a linear formula was then 
made for each set of differences (Av and Ay), using all the stars on 
P I. With the aid of the resulting constants, the motions uw, and p, 
were derived. The stars of large motion were removed, and the dif- 
ferences (Ax and Ay) of the remaining stars were used in a final solu- 
tion of the reduction constants and in the derivation of the motions 
of the stars on P I. These comparison stars were used whenever pos- 
sible on the remaining three plates. 

The removal of the magnitude equation was accomplished in the 
following manner: For all stars common to all four plates the mo- 
tions in the x- and y-directions were averaged and then plotted. From 
the resulting concentration near the origin, a number of stars closest 
to the center of gravity were taken as being cluster stars. Then for 
each plate the individual values of uw, and wy, were plotted against the 
magnitude, and through the points on each plot a straight line was 
drawn. For five of the resulting eight plots the line had an appreci- 
able slope, and the average value of the magnitude equation so de- 
rived was about 070025/year mag. The magnitude equation for the 
stars on P Land P II fainter than 11.7 was obtained by extrapolating 
the straight line defined by the brighter stars. 

After the motions from the four plates had been corrected, they 
were then in a form to permit them to be combined after the weights 
of the individual plates had been determined. In order to do this the 
differences PI- PH, PI-— PIII, PIl-— PIII, PHI - PIV 
were formed for all stars common to all four plates and for the mo- 
tions in both x and y. Each of the eight groups of differences was 
divided into two subgroups according to whether the star was bright- 
er than or fainter than m,, = 11.2. Then, for each of the sixteen sub- 
groups the quantity Yv?/(# — 1) was formed. For none of the sub- 


SA p. J., 89, 431, 1939 
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TABLE 2 

















Star No. | Mov Spec. Mer | My Wt. 
2, | ee Mem eru rete + 52 — 19 I 
2: | ? ey ae ee 116 fe) 2 
4 $209 122 | = 21 + 14 2 
S. 220° Es | =o a ate {> I 
6 29 «Is = 30 = 5 2 
8. 7 ad — 28 + 48 2 
12 [300 4 + 27 =. 6 I 
10. 12.0 |. = 20 oe 2 
18. 2S | “a as + 44 2 
2 2:5 | = 22 + 3 2 
22. 12.8 + 20 + 30 2 
23. f244.. | | — 86 +145 2 
24. 13.0 | —= 17 aie. I 
a6. . 12.8 } =i. | aot ae 2 
29. . 13.2 je oe ae et I 
| 
30 12.0 | F2 — 20 | + 3 | 3 
oo 13.0 | a a7 =) I 
30. i. + 68 eae ae I 
ee : 
37-- 12.7 | + 39 33. | 2 
39. ‘3.3 | + 26 + 29 | I 
| 
Ax, isin | +o) +a I 
44. 13.1 — 22 | — 28 | 2 
45. 13.0 + 45 [ oO | 2 
46. 13.0 + 32 + 18 | 2 
47 11.5 +126 —198 | 3 
48 2.7 “+ Ar | — 6 | 2 
49. 12.6 + 96 | + 36 | 2 
51. t2..7 Go | — 20 | — 17 | 2 
" | 
53 13.2 7 20° | a a I 
54. [2.9 ° ce 2 
Coe 11.6 A4 + 46 ye 4 
56 [3:2 + 99 aie I 
65. 12.8 F4 + 16 — 20 2 
68 13.2 +129 } + 75 I 
69. 13.2 - 67 a a. I 
70 “te Ry F2 +103 | + 90 3 
72 12.4 K - 24 apes I 
78 13.3 3 + 20 I 
78 13.6 +104 + 6 I 
79 r2.0 = 2 —- 40 4 
81 10.2 Bo ° + 18 4 
83 13.0 ites + 31 — 29 I 
84 10.6 gKr + 10 — 2 4 
88. 12.4 Go — 4 — 12 2 
goa. II.2 F3 + 2 + 25 4 


Membership 


wnNKR SPW 


N NM Ww 


n~ahr 


“Ne 


m th 
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TABLE 2—Continued 





























Star No. Moy Spec Mz | My Wt. Membership 
908. | RO Ree aan + I | — 10 I 2 
OI... | 12.7 F2 + 42 | — 30 2 3 
g2a. | UP a Sa Barer gee meee + 14 + 27 2 2 
04.. | 329 F8 + 24 + 18 2 2 
95. 12.8 F2 ° + iI 2 2 
06... 12.6 F8 + 29 + 3 2 2 
98 12.7 Fs + 3 — 3 2 2 
IOI. 1250 beinss + 17 + 25 I 2 
102. Ly UE Cee —= 30 — a 2 2 
104.. II.4 gGs + 4 — 19 3 I 
105 | 10.4 gKo — 5 — 6 4 I 
108 9.9 gK3 — | + 8 4 I 
II! 12.6 errr — 4 — I a 2 
112 ae i: ae ene ey + 20 — 6 I 2 
115 | 412.7 F6 — I — 4 2 2 
116 12.9 + 10 — 22 I 2 
TI7. 12.6 + 22 + 3 2 2 
119 12.8 F8 ° + 28 2 2 
124. 12.2 Fy — 2 + 28 3 I 
127 2.7 — 10 — 25 2 2 
128. 13.0 M- — 30 — 34 2 3 
130. 12.9 G: + 20 —- 9g 2 2 
131 ke Fr + 19 ° 4 I 
£23. 13.3 G: — 16 — 27 2 2 
134 | -32.3 Fs + 22 ett 3 I 
135 | brs gKo + 18 + 36 3 2 
136 11.4 Fo + 18 — 4 4 I 
141 10.4 gKo — 13 + 9 4 I 
143 11.6 gG8s + 7 — 10 4 I 
144 12.9 + 82 — 2 2 3 

} 
145 12.8 Go: — 52] — 10 2 2 
147 13.0 Go: + 22 | — 32 2 2 
148 13.0 Go | + 22 + 20 2 2 
149 12.7 Fs + 16 + 4 2 2 
151 10.6 Go | + 18 + 4 I 
| 
153 11.5 3 — 9 — 2 4 I 
155 10.5 gGo | — 49 > 9 4 3 
156 FI.2 | Al ) = eee | — 12 4 I 
157 ‘2.9 F8 + 14 — 4 2 2 
160 12.8 — 7 +158 2 4 
161 12.8 | F3 + 10 | + 20 2 2 
162 12.8 | F8 + 4 + 32 2 2 
163 ‘2.7 Go + 6 — 3I I 2 
164 10.8 | gGo + 16 — 7 4 I 
165 | 12.7 | Fs + 53° | — 12 2 3 





| 
| 
| 
| 
| 
| 
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TABLE 2—Continued 

















| | | | 
Star No. Mov Spec. | My My Wt. | Membership 
| | 
166 a4 12.8 Gs + 40 — 24 2 2 
170 rai 9.8 gK2 —- 7 + 2 4 I 
171 : 12.9 + 58 + I0 I | = 
172 12.6 + 6 — 36 2 2 
173 122 F8 + 28 — 3 3 | I 
| | 
173 13.0 | + 26 — 4 2 2 
174 | ‘see | F8 + 39 — 6 2 2 
175 | 13.4 G: +31 | -— 4 I 2 
176 12.7 F8 + 25 | — 22 2 2 
180 12.6 F5 — 10 — 22 2 | 2 
181 12.8 + 30 — 10 2 | 2 
182 | oe ae Loe Fy ° 2 2 
184 2 | Fo << = i + 21 | 2 I 
185 | oe A4 | + 9 | + 2 | 4 | I 
185 | 7 | F6 | + 20 | ~ a 4 I | 2 
| 
187 13.0 F8 | — 17 + 6 | 2 2 
189 | 432.7 aati 2 2 | 2 
190 11.2 A8 | emi =~ 4] 4 I 
1g2 12.8 + 16 | o | 2 2 
193 12.3 | gGs5 | TT. | 6 | 2 2 
| | 
195 Sy aa Go | + 42 | T 23. 2 | 2 
199 13.0 | 40 | — 15 2 | 3 
200 io | G8 —172 | + 5: 4 4 
202 | 249 4 Fs | +20 | — 36 2 2 
203 12.9 ; +18 | + 10 2 2 
| } 
206 b eaay | —198 | — 20 } 
207 12.3 sy — 12 3 3 
210 £2.23 F8 Se 3 I 
215 ae, F5 + 20 | + 27 y 2 
216 / 33:6 F5 | + 80 | — 24 | ? 3 
217 ae ae Ko | +22 | + 18 1 i 
217 f2.2 — 4 + I! 3 I 
218 iF + 6 | — 6 4 I 
219 12.6 | + 60 | ++ 23 | 2 3 
221 12.6 F8 + 68 | — 98 | { 
| | | | | 
223 | 10.7 | gG8 } +12 | 7 { I 
224 | ET] [oes oe a + 4 3 I 
225 12.9 | + 10 | = 66. | 2 > 
226 a7 4 Go + 1 | — 22 2 
227 12.7 Go + 27 + 10 2 2 
225 12.9 - 9 | — 4 | I 3 
220 12.8 F2 | +156 | —139 | 2 4 
231 | 2.0. |. @Ko = 8 | o | 3 I 
234 : ee. | | — 27 | —509 | 3 4 
236 12.4 Go | + 43 — i | I 2 
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TABLE 2—Continued 




















| | 
Star No. | Moy | Spec | My | My Wt Membership 
i Sree) Ros. 
- =" a - 
237 > aa : | + 2 | + 27 2 2 
238. 11.2 | A3 | +16{| + 6 4 I 
241 ¥3.7 ; + 28 | — 6 2 2 
242 8.0 gKo — 60 | + 20 3 4 
243 12.7 F8 | ui + 2 2 2 
244 tke (| Go | -— 2] + 3 4 I 
248 12.3 | ‘2 | — 12 | ° 2 2 
255 [2.7 Go | + 28 | — 6 2 2 
256 12.6 + 18 | ° 2 2 
257 ri6- | F8 —128 | —366 4 4 
| | 
} 
61 m.7 . | A3 + 12 — 78 4 4 
262 12.7 | | + 59 | + 26 I 3 
265 12.5 F8 — | Page 3 4 
266 10.7 | gGs + 13 | — 24 4 I 
207 ee | + 16 + 62 2 3 
269 } 13.0 + 26 + 16 I 2 
271 12.8 | + 33 — 22 | 2 2 
277 | mo | “— 74 +110 | 2 4 
280 } 11.0 A2 + 33 + 34 | 2 2 
281 } 12.9 + 30 ae Oe | I 2 
285 2.5 +210 | — 86 2 4 
286 | 10.7 gG5 | <—ogn | Roe 4 I 3 
| | i | 
289 | 12.7 . eae | a et 3 I 
| | | | 





groups did this quantity differ significantly from the mean of all, and 
hence all plates were given unit weight regardless of co-ordinate or 
magnitude. On this basis a weight of four corresponds to a mean 
error of +0*%0012/year in either co-ordinate. 

With the aid of these weights the motions from the various plates 
were combined and are given in Table 2. 

The columns give, respectively : 

t. The number according to Fagerholm’s catalogue.* His paper 
also gives a chart of the cluster. 


to 


The photovisual magnitude according to Shapley.’ His paper 
also gives photographic magnitudes. 

3. The spectral type, kindly communicated in advance of publica- 
tion by Dr. R. J. Trumpler. 

4 and 5. The proper motions in x and y expressed in units of 
0’ ooo1 / year. 
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6. The weights of the motions in the preceding columns. Weight 
4 corresponds to a mean error of +0%0012/year. 

7. A measure of the probability that a star is a member of the 
cluster. Unity indicates that the star has a high probability of being 
a cluster member; 2 indicates that the probability is lower but that 
the star has a reasonable chance of being a member; 3 signifies that 
the star is probably not a member; and 4 denotes that it is definitely 
not a cluster member. The assignment of weights is discussed in sec- 
tion IIT. 
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FIG. 1 
III. DISCUSSION OF THE RESULTS 
The data in Table 2 are complete down to m,, = 13.0, within a 


radius of 15’ from star No. 155, which was considered by Shapley to 
be the center of the cluster. 

Figure 1, A, is a plot of the motions of the stars with a weight of 
3 or 4, and Figure 1, B, is a plot of those whose weights are either 
EOF 3. 

The assignment of probability of cluster membership was made in 
the following manner: The average weight of the stars in Figure 
1, A, corresponds to a mean error of +0%0014/year. About the cen 
ter of gravity of the concentration three circles were drawn with radii 
equal to 2, 2+/2, and 4 times the mean error just given. Stars within 
the smallest circle were given class 1 membership, those between the 














oo 
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first and second circles class 2, those between the second and third 
circles class 3, and those outside the largest circle class 4. The aver- 
age weight of the stars in Figure 1, B, corresponds to a mean error of 
+o%oo21/year. In this case the circles were given radii of two and 
four times the mean error just given, and no class 1 membership was 
assigned. Those stars within the first circle were given class 2 mem- 
bership and those outside the second circle class 4. 

Figure 2, A, gives the position in the cluster of all stars with 
weights 3 or 4, and Figure 2, B, gives the positions of all whose 
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FIG. 2 


weights are 1 or 2. In both cases the filled dots correspond to stars 
whose membership class is 1 or 2. 

The appearance of the cluster on a photographic plate is that of a 
nucleus of bright stars about 16’ in diameter, superimposed on a 
background of fainter stars. On P III and PIV only the nucleus 
of brighter stars was present, while on P I and P II a considerable 
number of fainter stars outside the nucleus were measured also. For 
this reason, stars with weight 3 or 4 are concentrated toward the 
center of the measured area, and stars of weight 1 or 2 are spread 
over the whole measured area. On the basis of star counts Shapley 
has found that the cluster of bright stars is ‘merely a well marked 
nucleus of brighter and redder stars in a much larger system,” which 
he estimates may extend to a radius of as much as half a degree from 
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the center. This is the ‘‘shoulder”’ effect which has been found for 
several clusters by Trumpler. 

Figure 3 is the magnitude-spectrum diagram for the cluster. The 
ordinates are bolometric magnitudes, and the abscissas are given as 
both log Te and spectral type. The open circles refer to stars whose 
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membership class is 3 or 4. Attention is called to star No. 193, whose 


position in this diagram is log Te = 3.66 and my.) = 12. The dis- 
tance modulus being 9.3 mag., the star should have an mj) = +2.7 


which, if the star belongs to the cluster (membership class 2), makes 
it rather faint for its spectral type, gGs5. 
The brightest star in the cluster, No. 242, appears not to be a 
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member. However, since it is 2 mag. brighter than the next faintest 
star in the cluster, the corrections for magnitude equation had to be 
extrapolated, and therefore the corrections are very uncertain. The 
corrected values of motions as derived from the three plates on which 
the star was measured are quite discordant and their mean cannot 
be very accurate. Since the magnitude equation is uncertain, it 
might be better to use the uncorrected values. These happen to be 
very much more in accord with one another, but their mean still gives 
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the star a class membership of 3. The considerable deviation of this 
star from the center of the cluster (7’) might also be considered as 
evidence that it does not belong to the cluster. The meridian proper 
motion of this star as derived from only two observations is wz = 
—o"’037/year and wy, = —o.006/year. 

The spectrum-magnitude diagram of M 67 resembles, to a con- 
siderable extent, that obtained for globular clusters, in that a promi- 
nent late-type giant sequence is present. 

Since the late-type giants have a smaller concentration toward the 
galactic plane than the B- and O-type stars, those clusters which con- 
tain late-type giants might be expected to behave in the same fashion 
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with respect to clusters that have no such giants but whose brightest 
stars are of types O, B, and A. To examine this expectation, the ga- 
lactic clusters in Trumpler’s® list were divided into three groups, as 
follows: first, those whose Trumpler type is 1, i.e., no giants; second, 
those whose type is 1-2, i.e., uncertainty exists as to whether the 
giants belong or not; and, third, 2, 2-3, or 3, i.e., giants are certainly 
members of the cluster. For each of these groups using Trumpler’s 
values for the z-distances the arithmetic mean was formed, and these 
values together with their standard deviation are for the three 
groups, respectively, 48 + 15, 60 + 16, and 102 + 32 parsecs. In 
order to avoid selection effects, only clusters between 200 and 1000 
parsecs were used in the computation. The run of decreasing con- 
centration toward the galactic plane is in the direction to be expected 
for early-type and late-type stars not belonging to clusters. 

The luminosity function is given in Figure 4. The rapid increase 
of stars beginning at m,, = 12.0 may, toa small extent, be due to the 
fact that the average motion of the cluster stars differs but little from 
the mean of the noncluster stars, and some noncluster stars may be 
included as members. This is partly offset by the fact that some stars 
which were given class 3 membership may actually belong to class 2. 


IV. THE PROPER MOTION OF THE CLUSTER 

Two previous attempts have been made to determine the proper 
motion of the cluster; one by Fagerholm’ and the other by van 
Rhijn.4 Fagerholm compared the positions of 118 stars common to 
his catalogue and that of Olsson’s,* made eight years before, but 
found no distinct evidence of motion. 

From two carte du ciel plates of the cluster with an interval of 27 
years, van Rhijn determined the absolute proper motion of the clus- 
ter to be 0%0073/year in position angle 287°. The mean error of an 
individual star’s motion in either co-ordinate is +o’009/year. Al- 
though this direction differs by 13° from the direction of the conver- 
gent for the Hyades (274°), van Rhijn thought that the two clusters 
might have the same motion in space. None of the stars measured in 
the region of M 67 is found in the Boss General Catalogue, and only 


© Lick Obs. Bull., 14, 154, 1930. 
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one star, No. 242, isin Schorr’s Eigenbewegungs Lexikon as mentioned 
before. 

An attempt was made to determine the proper motion from the 
measures of faint stars in the cluster. The motion of the cluster rela- 
tive to 29 faint noncluster stars whose mean photographic magnitude 
is 13.6 was computed to be uz = —o%0017 + o%0014/year and py, = 
—o"”o008 + o%0012/year (m.e.). For stars of mp = 13.6 at galactic 
latitude 33°3, Oort gives a mean parallax of o”0021 + 0o”%0002 (m.e.), 
where, because of the few stars used in our case, the mean error is a 
lower limit. The antapex is 58° from the cluster, in position angle 
214°. On the basis of a solar motion of 20 km/sec, the parallactic mo- 
tion of the group of faint stars used is 00075 + 0%0008/year in posi- 
tion angle 214°. Hence, the proper motion of the cluster so deter- 
mined is 07009 + 0%002/year in position angle 220° + 10° (m.e.). 
On the basis of this determination, the difference in direction be- 
tween the motion of M 67 and the direction of the Hyades conver- 
gent is 54°, and hence van Rhijn’s suggestion seems improbable. 

It is of interest to reinvestigate the space motion of the Praesepe 
cluster in order to determine whether or not it is parallel to that of 
the Hyades. Of the stars listed by Klein Wassink’ as being members 
of the Praesepe cluster, 33 were found to have proper motions in the 
General Catalogue. The weighted mean proper motion is pz = 
—o"034 + o%001/year and wy = —o’018 + o%oor/year or a total 
proper motion of 07038 + o’0o01/year in 242°1 + 1°5. In Moore’s 
catalogue of radial velocities 6 stars were found. Two of these are 
binaries; for one there is no orbit, and for the other a determination 
has been made, and hence its space motion is known. The weighted 
mean of the five determinations of radial velocity is +33.2 + 1.3 
km/sec. The parallax of the cluster as derived by Trumpler is 
0”0067. On the basis of these data the calculated position of the con- 
vergent is a = 95°2 and 6 = —o°6, with a mean error of about one 
degree in either co-ordinate. The computed space velocity is 42.9 
+ 1.5km/sec. The distance between this convergent and that of the 
Hyades as determined by Smart is 9°5. An attempt was then made 
to reduce this distance by changing the adopted parallax of Praesepe. 
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By changing the parallax to 070072, the distance between the con- 
vergents was reduced to 9°3, but no more. Hence it would seem im- 
probable that the magnitude and direction of the space motion of the 


two clusters is the same. 


I wish to express my thanks to Dr. Struve for having placed the 
facilities of the Yerkes Observatory at my disposal during the sum- 
mer of 1939 and gratefully to acknowledge the advice which Dr. 
Kuiper has given to me. 


YERKES OBSERVATORY 
September 1939 
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THE HYPOTHESIS OF THE EXISTENCE OF 
CONTRATERRENE MATTER 


V. ROJANSKY 


ABSTRACT 

Symmetry arguments suggest the possibility of existence of atoms (here called 
“contraterrene”) consisting of negatively charged nuclei surrounded by positrons. In 
this note it is assumed that a fundamental contraterrene particle and the corresponding 
ordinary (here called “‘terrene’’) particle annihilate each other on collision. Annihila- 
tion effects should make it possible to identify contraterrene matter (which presumably 
cannot be distinguished from terrene matter spectroscopically) provided that it can 
be observed in the process of colliding with terrene matter. Within the solar system 
perhaps the only bodies that are not certainly terrene are some meteors and the comets. 

The quantum theory of the electron suggests a symmetry between 
positive and negative electric charges, and it has been conjectured 
that in some parts of space there might exist atoms whose nuclei are 
charged negatively and whose extranuclear members are positrons.* 
The purpose of this note is first to put the speculative notion of such 
atoms into a reasonably clear-cut form and then to consider briefly 
some of the astrophysical implications that their existence might 
have. For brevity, we shall call these atoms and their constitu- 
ents and aggregations “contraterrene,”’ and call ordinary bodies 
“terrene.”’ 

Guided by symmetry arguments, we shall assume that the contra- 
terrene neutron, unlike the ordinary neutron, has parallel spin and 
magnetic moments and that the contraterrene proton and the ordi- 
nary proton are symmetric to the same extent that the positron and 
the electron are. The spectra of an ordinary atom and of the cor- 
responding contraterrene atom would then be identical, and the 
possibility of direct spectroscopic identification of contraterrene 
matter would be ruled out. At this point, however, symmetry con- 
siderations suggest a further hypothesis, namely, that a fundamental 
terrene particle and the corresponding contraterrene particle destroy 
each other on collision, just as an electron and a positron are known 

tP. A. M. Dirac, Die moderne Atomtheorie (by W. Heisenberg, E. Schrédinger, and 
P. A. M. Dirac), p. 45, Leipzig: Hirzel, 1934; G. Gamow, Structure of Atomic Nuclei 
and Nuclear Transformations, p. 14, Oxford, 1937; V. Rojansky, Phys. Rev., 48, 108, 1935. 
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todo. This hypothesis introduces the possibility of detecting contra- 
terrene matter if it exists, provided that it can be observed in the 
process of colliding with terrene matter. The radiant energy, given 
by Einstein’s equation E = mc’, released in an electron-positron 
annihilation is shared by two photons, each of about half a million 
electron volts. In a proton-proton (one terrene, the other contra- 
terrene) annihilation the energy per photon would be about a billion 
electron volts. A neutron-neutron annihilation would either yield 
directly a pair of billion-volt photons or reduce to a proton-proton 
and a positron-electron annihilation. Since electrons attract positrons 
and ordinary atomic nuclei would attract contraterrene nuclei, the 
mutual annihilation of terrene and contraterrene bodies in contact 
would be rapid. 

We now turn to possibilities of collisions between terrene and 
contraterrene bodies and begin with the earth. If contraterrene 
matter exists in some parts of space, then some of the meteors that 
have their origin outside the solar system might be contraterrene. 
A contraterrene meteor that does not escape after entering the at- 
mosphere would become completely converted into radiation; in par- 
ticular, if it should reach the ground, it would be annihilated ex- 
plosively; hence we must conclude that if contraterrene meteors 
exist at all, those that are large enough to penetrate deeply into the 
atmosphere are extremely rare.” Incidentally, the capture of contra- 
terrene meteoric matter by the earth at the rate of 1 gm per square 
mile per year would result in the release of about 0.003 gm cal per 
square centimeter per minute. 

The effects of a contraterrene meteor in the upper atmosphere 
would differ from those of a terrene meteor for two main reasons: 
the high penetrating power of the primary products of annihilation 
and the large quantity of energy released per unit mass. In particu- 
lar, a contraterrene meteor would display ionizing power which is 
very great compared to that of a terrene meteor of the same mass 
and speed, and some ions would appear far from the path of such a 

2 If we disregard the effects of heating and make the rough assumption that a contra- 
terrene meteor annihilates all the air molecules which it would strike if it did not affect 
their motion, we find, for example, that unless the initial top-to-bottom dimension of a 


contraterrene iron meteor, falling vertically without tumbling, exceeds 130 (i.e., 76X 
13.6/7.9) cm, the metzor would be entirely radiated away before reaching sea-level. 
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meteor sooner than they could arrive there by diffusion. Therefore 
contraterrene meteors, if they exist, can probably be detected by 
experiments combining ionization measurements over large portions 
of the atmosphere, like radio signal observations, with measure- 
ments, made in the upper atmosphere, of the energies of the annihi- 
lation products. Since the energy of the primary photons would be 
dissipated in the lower atmosphere through mild showers, it is not 
likely that the release of these photons could be ascertained on the 
ground by means of a single detector of a cosmic-ray type. Without 
a detailed study it is of course not possible to decide whether a con- 
traterrene meteor could be identified by visual observation. 

Next we shall consider a hypothetical contraterrene body revolv- 
ing about the sun. Were it moving in a void, it would behave in a 
very ordinary way; but since it would in fact be bombarded by 
terrene meteoric material, it would undergo annihilation, transmu- 
tation, and mechanical disintegration. Annihilation follows from our 
fundamental premises. Transmutation is a concomitant of annihila- 
tion; it is certain to occur, for example, when a terrene atomic 
nucleus of atomic number A is annihilated in collisions with contra- 
terrene nuclei of atomic number B which is not an integral submulti- 
ple of A, and it is likely to occur quite generally, since nuclear 
annihilation would consist of individual proton-proton and neutron- 
neutron annihilations. Mechanical disintegration, in addition to 
structural changes that might be caused by the sun through tides 
and heating, is to be expected in view of the large amounts of 
energy released locally by the annihilation of the larger meteoric 
fragments. Because of transmutation the chemical composition of 
the body would be gradually shifting toward the lighter elements; in 
particular, volatile compounds might thus be continually produced 
and released. 

If an object having these characteristics should be observed from 
the earth, it would probably be classified as a comet; and the ques- 
tion, therefore, arises whether the comets, or at least some of them, 
may not actually be contraterrene bodies. The hypothesis of the 
contraterrene structure of comets implies the idea that, although 
the comets are now genuine members of the solar system, they had 
their origin outside the system, an idea that was first suggested for 
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quite different reasons and is now well known. The present hy- 
pothesis, however, does not require that, when the comets were first 
acquired by the solar system, they had, apart from being contra- 
terrene, any unusual chemical composition or mechanical structure. 
On the basis of this hypothesis, meteors of cometary origin become 
of special interest; the existence of a meteorite that is definitely 
known to be of cometary origin would of course prove that the 
parent comet was terrene. 

The following computation, suggested by Professor John A. 
Wheeler, is of interest. If a contraterrene body, bombarded uni- 
formly by terrene matter but not appreciably heated by the sun, 
annihilates all this matter and absorbs half of the annihilation 
energy, its equilibrium temperature, computed in round numbers 
from the Stefan-Boltzmann law without regard for transmutation 
and mechanical effects, turns out to be 


7° > 4 
T = 1800V uz, 


where yu is the rate, expressed in grams per square mile per minute, 
at which terrene matter impinges upon the body. To illustrate, we 
shall make two estimates of the equilibrium temperature of a contra- 
terrene body, not heated by the sun, bombarded by terrene matter 
at the same rate per unit surface as that at which the earth is bom- 
barded by meteors. Nininger’ gives the rate of capture of meteors by 
the earth as 50,000 tons per day; on this basis u is about 2 X 107 
and T is about 1200° absolute. Wylie4 gives the rate as 10 gm per 
square mile per year; on this basis yw is about 2 X 10% and T is about 
12c° absolute. 

The idea of contraterrene matter suggests implications on several 
levels of the astronomical scale. We shall, however, stop with the 
few phenomena that seem most likely to provide an early proof 
either of the existence of this matter or of its nonexistence within the 
solar system. 

UNION COLLEGE 
SCHENECTADY, N.Y. 
3 Our Stone-pelted Planet, p. 91, Boston and New York: Houghton Mifflin Co., 1933. 


4 Phys. Rev., 47, 192, 1935. 
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Ha [Or] 3727 
A, Prime-Focus PHoTtoGRAPH OF M 101 TAKEN WITH THE 82-INCH 
McDONALD REFLECTOR 
Positions of the emission nebulae are indicated by arrows. Ten-minute exposure on 
Cramer Hi-Speed plate. Contact print. Scale, 25”6 per mm. 


B, ENLARGEMENT FROM THREE FILMS TAKEN WITH THE CASSEGRAIN SLIT SPECTRO 
GRAPH OF THE 82-INCH REFLECTOR SHOWING (a) OByect NGC 5471 IN M 101 
(EXPOSURE 90 MINUTES), (0) ORION NEBULA (ExposurRE = Minutes). (c) PLANI 
TARY NEBULA NGC 6543 (ExposuRE 2 MINUTES) 

Agfa Super Pan Press film used for all three spectrograms. Enlargement 17 

Dispersion on original spectrograms is 345 A per millimeter at //} 


EMISSION NEBULAE IN MESSIER 101* 
CARL K. SEYFERT 


ABSTRACT 

Ten emission objects in the spiral nebula M 1o1 have been found on slitless spectro- 
grams taken with the 82-inch telescope of the McDonald Observatory. Attention is 
directed to one of these emission objects in particular, a small group of stars which 
appears to have an emission spectrum of an unusually high degree of excitation. 

Patches of emission nebulosity have been found in several extra- 
galactic nebulae, principally by Edwin Hubble,’ by N. U. Mayall and 
L. H. Aller,? and by H. W. Babcock.’ A to-minute exposure with 
the slitless spectrograph of the 82-inch reflector of the McDonald 
Observatory revealed ten emission objects in the large late-type 
spiral M ror (NGC 5457; 14"1™4, +54°35' [1950]). The positions of 
the emission patches are indicated by arrows on Plate II, A. In the 
table are listed: (1) the current number; (2) the NGC designation 
(the most prominent condensations of M 1o1 have NGC numbers); 
(3 and 4) the rectangular co-ordinates of the emission objects, referred 
to the nucleus of the spiral; and (5) the emission lines observed in 
the order of their apparent strength on the film (Agfa Super Pan 








Press). 
Current No. NGC No. | AX AY Emission Lines 
I | (—301” —262”" | Ha 
5447 | 1-378 | —277 | Ha 
3 = 309 | —284 | Ha 
{ 5455 —103 | —387 | Ha, Nx, 3727 (O I) 
: +126 | —134 | Ha 
( +164 | + 17 Ha. 
5401 +250 | —II3 | Ha, N,, 3727, HB, Hy, N., Hé, He 
| 
8 sci +346 + so | Ha 
9 . +359 | + 68 | Ha 
er 
I 5471 | +670 | +170 | Ha, N,, N2, HB, Hy, 3727 
| 





* Contributions from the McDonald Observatory, University of Tex1s, No. 17. 
1 Mt. W. Contr., No. 304, 1925; No. 310, 1926; Ap. J., 62, 409, 1925; 63, 236, 1926. 


2 Pub. A.S.P., §1, 112, 1939. 
Reported at the A.A.A.S. meeting in Alpine, Texas, in May, 1939. 
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Ha is, without exception, the strongest feature of each of the 
emission spectra; in fact, it is the only line visible for all but the three 
brightest condensations. Among these latter objects, the apparent 
relative strengths of the lines on the film for NGC 5455 and 5461 are 
characteristic of galactic emission nebulae as typified by that in 
Orion. Ha is very much stronger than N,; the O m doublet \ 3727 is 
quite strong; and for NGC 5461 in which the lines are visible, 78 
is stronger than V,. These apparent relative intensities are, of course, 
only rough indications of the true relative intensities, since the ap- 
parent values are strongly influenced by the type of instrument and 
emulsion employed. 

On the other hand, in NGC 5471, NV, appears almost as strong as 
Ha; N, is definitely present; but only traces of HB, Hy, and d 3727 
could be detected. Hence, it is probable that the degree of excitation 
of this object (NGC 5471) is considerably greater than that of the 
other condensations in emission which have been found in M tot. 
The high excitation of NGC 5471 is further illustrated by Plate II, B, 
which shows a slit spectrogram of this object, together with spectra 
of the Orion nebula and the planetary nebula NGC 6543 for purposes 
of comparison. All these spectra were taken on the same type of 
emulsion (Agfa Super Pan Press), with the same instrument (Casse- 
grain spectrograph of the 82-inch reflector; f/1 camera). It will be 
seen that the spectrum of NGC 5471 is not unlike that of the plane- 
tary nebula. Relative to the Balmer lines, the forbidden O 111 lines 
(NV, and N,) and the Ne 11 line at \ 3869 are considerably stronger 
in NGC 5471 than the corresponding lines in the Orion nebula. Us- 
ing, as an excitation criterion, the ratio N,/3727 (i.e., O ur/O 11), 
we see from Plate II, B that this ratio is approximately one for NGC 
5471, whereas it is considerably less than one for the Orion nebula 
and greater than one for the planetary. Hence, the degree of excita- 
tion of NGC 5471 is intermediate between that of the diffuse nebula 
in Orion and the planetary nebula NGC 6543. 

NGC 5471 is well separated from the main body of the spiral; and 
the observed radial velocity, as determined from the slit spectro- 
gram (+500 km/sec), is somewhat greater than that of the spiral as 
determined by Humason‘ at Mount Wilson (+300 km/sec). In 


4 Handb. d. Ap., 7, 555, 1936. 
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spite of this, it is probable that NGC 5471 is a part of M 1ot, since 
the radial velocity of the former is far too large for it to be a member 
of our own galaxy and too small for it to be a separate group of ex- 
tragalactic nebulae showing emission lines. 

On the prime-focus photograph NGC 5471 looks like a close 
group of five or six objects of semistellar appearance with no appreci- 
able nebulosity involved. The other patches of emission, on the 
other hand, are definitely nebulous in appearance and are all within 
the main body of the spiral. NGC 5471 has a diameter of 15’, or 
about 40-50 parsecs, as against a diameter of 10-25 parsecs for the 
nebulous emission knots. At least two of the objects comprising the 
cluster NGC 5471 appear in emission on the slitless spectrogram, to- 
gether with a faint continuous background, which can be seen on the 
original film. 


McDONALD OBSERVATORY 
September 15, 1939 
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A METEOR SPECTRUM OF HIGH EXCITATION 
On the morning of October 20, 1939, in the course of an exposure 
made with the 10-inch Cooke prismatic camera (f = 1:4.5),a very 
bright meteor appeared in the field. The meteor was not observed 
visually, but its velocity must have been high, as a very short flash 
of light was noticed inside the building. A bright train was observed 
about 10° long; at first it was straight, but subsequently it curved ir- 
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regularly; it drifted backward about 10° and faded away after 5 min- 
utes of total duration. 

On the plate the spectrum first appears (Pl. III) at a = 3556", 
5 = +42°9, where it is barely visible. The intensity increases con- 
tinuously until the edge of the plate at a = 3"24™, 6 = +46°0 
(1939.0). As the direction of the meteor path traced backward 
passes close to the radiant of the Orionids, the meteor may be safely 
identified with this shower. It appeared at 4"7™ local sidereal time, 
or 7530™, G.M.T. 

The spectra on the plate are short, the dispersion being 350 A/mm. 
The positions of sixteen lines were measured on the Gaertner ma- 
chine, and their intensities were estimated at five different points. 

t At this point the apparent magnitude was about — 7.5, as determined from trails 


of Jupiter. 
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METEOR SPECTRUM, OCTOBER 20, 1939 (ORIGINAL SIZE) 
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The measured wave lengths are given in the table and the intensities 
are represented schematically in the figure. 





Line A Meteor 
I 3820 
2 3850 
3 3880 
+ 39390 
5 397° 
6 4040 
7 | 4070 
8 4140 
9 .| 4240 
10 4280 
Il 4390 
12 4490 
| 
13 4050 
| 
14 | 4950 
15 5060 
ays 
16 | 5200 








Element 


F e 
Fe 


Fe 
Mg 


Fe 


os" 
Fe*: 


Fe 
Mg 
ir 
Fe*: 














X Laboratory 


3820, 3826 
3815 


3860 
3838, 3832, 3829 


3579 
3934 
3968 


4046 
4031, 4933, 4034 


4064, 4072 


4131, 4125 


4227 
4230, 4250 

4233 

4254, 4275, 4289 
4251, 4272, 4308, 4326 
4353, 4405 

4481 

4508 

4029 


4646, 4651, 4652 


4958, 4920, 4891 


4924 

5057, 5041 
5018 

5192, 5233 
5183, 5173, 5107 
5208, 5206, 5205 
| 5169 





| 
| 
| 





Remarks 


Considerably out of 
focus 


Slightly out of focus 


Faint, fuzzy, in the fog 
of lines 4 and 5 at in- 
tense end 


Very sharp all through 


Fairly sharp, but 
blended 


Fairly sharp at the be- 
ginning; broad and 
fuzzy at the end 

Sharp at the beginning; 
possibly double at 
the end 


Faint and fuzzy 


Sharp 


Fuzzy 


Fuzzy 


It should be noted that some lines, such as 6, 7, and 10, appear 


from the very beginning, yet they never become very intense; others, 
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such as 4, 5, 8, and 12, while faint or invisible at the start, are very 
intense at the end. Evidently the lines of the first group are of low- 
excitation potential, while those of the second group have distinctly 
higher excitation potentials. In fact, the lines at the intense end are 
very similar in position and relative intensity to those in the spec- 
trum of an A2 star, with the exception of the hydrogen lines, which 
do not, of course, appear in the meteor spectrum. This similarity has 
been used as a guide in the identifications which are given in the ta- 
ble. Thus, the lines of the ionized elements all show a rapid increase 
in intensity compared with the behavior of the neutral elements. 
Compare, for instance, line 12 of Mg* with 10 of Fe, or, again line 8 
of Si+ with 9 of Ca. The unusually high excitation of this spectrum 
is doubtless due to the high velocity of the meteor, as is evidenced 
by the almost instantaneous flash and by the appearance and long 


duration of the train.” 
A. N. VYSSOTSKY 
LEANDER McCorMIcK OBSERVATORY 
UNIVERSITY OF VIRGINIA 
October 25, 1939 





NOTE ON THE SURFACE TEMPERATURE 
OF VENUS 

Recent measurements of the total radiation emitted by gaseous 
carbon dioxide afford an interesting application to the determination 
of the surface temperature of Venus, by way of establishing an upper 
limit to the greenhouse efiect produced by the carbon dioxide con- 
tained in this planet’s atmosphere. Spectroscopic observations have 
so far failed to reveal the presence of water vapor or any other gases 
which are opaque in the infrared. By independent methods T. 
Dunham! and A. Adel? have estimated the amount of CO, in the 
atmosphere of Venus to be equivalent to 400 and 200 m, respectively, 
under standard conditions of temperature and pressure (T = 273° K, 
p =1atm). The uncertainty of these estimates is irrelevant to the 
argument to be given presently, as so large a mass of CO, has an 
optical depth practically infinite in the ranges of wave lengths 

2 Cf. with the conclusions of P. Millman, Harvard Ann., 82, 173-174, 1935. 

* Carnegie Institution of Washington Y ear-Book, No. 31, 154, 1932. 


2 Ab. J., 88, 345; 1037. 
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covered by the emission bands and is completely transparent be- 
tween these regions. 

The total emission of carbon dioxide has been studied by T. 
Dreisch? and E. Eckert.4 The work of the second author is quite 
comprehensive, his observations covering the range of temperatures 
between the boiling-point of water and 1200° C. The thickness of 
the radiating layer was varied between 2 cm and 3 m, and the carbon 
dioxide was diluted by varying amounts of nitrogen, while the total 
pressure was kept constant at 1 atm. Thus it was proved that the 
total emission of CO, obeys Beer’s law, i.e., the total emission is de- 
termined solely by the amount of CO, per square centimeter, which, 
incidentally, is not the case for water-vapor radiation according to 
Eckert’s findings. The validity of Beer’s law for CO, is a highly 
fortunate result, as it permits us to disregard the other hypothetical 
constituents of the atmosphere of Venus, such as nitrogen and the 
noble gases. Eckert has published his measurements of the total 
radiation Eco, as fractions ¢ of the total radiation of a black surface 
of equal temperature, so ¢ is defined by 

_ Eco, 
€ = ore . 


While ¢ increases with the thickness of the radiating layer, expressed 
in meter-atmospheres of CO,, it is practically independent of T for 
layers thicker than about 0.10 m-atm. In other words, a thick layer 
of CO, behaves like a gray radiator with an absorption coefficient e, 
which is independent of the temperature. Eckert’s observations 
extend up to 2 m-atm, and their extrapolation to 200-400 m-atm 
leads to the value € = 0.35 + 0.05. These limits represent, in the 
writer’s opinion, the uncertainty of both the extrapolation and the 
actual amount of CO, present on Venus. 

On a slowly rotating black planet the subsolar point reaches the 
maximum temperature T = 392°/V R, at a distance of R astro- 
nomical units from the sun. If the planet has a spherical albedo, A, 
for the incoming solar radiation, while behaving like a perfect radia- 
tor in the infrared, the maximum temperature of the subsolar point 

3Zs. f. Phys., 79, 711, 1932. 

‘VDI, Forschungsheft No. 387, Beilage zu Forschung auf dem Gebiete des Ingenieur- 


wesens, Berlin: VDI Verlag, 1937. 
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is reduced to T = 392V (1 — A)/R?, or T = 366°K for Venus 
(A = 0.60),5 by neglecting the effect of the planet’s atmosphere. 
Actually, the atmospheric CO, will intercept a certain part of the 
radiation emitted by the surface, which, according to Eckert’s re- 
searches, is equal to the fraction e characterizing the absorption co- 
efficient of a thick layer of CO,. The complemental fraction 1 — e 
can escape directly into space between the CO, bands. The raising of 
the surface temperature, resulting from the re-emission of the radia- 
tion intercepted by the atmosphere, would be most difficult to evalu- 
ate rigorously. However, an upper limit to the surface temperature 
follows in a simple manner from the energy balance of the total 


radiation, which reads 
(1 — €)T¢ + eT4 = 366'. 


On the left-hand side the first term is the radiation escaping between 
the CO, bands from the surface (temperature 75) and the second 
term is the emission of the atmosphere directed outward, which can 
be ascribed to a layer having the effective temperature 7. This 
outgoing radiation is balanced by the insolation, producing on the 
atmosphereless planet a temperature of 366° K, as shown before. 
For obvious physical reasons, 74 must be smaller than 75. Hence, 
an upper limit for 7's can be established by putting the term 74 
equal to zero. Physically this would be equivalent to the assump- 
tion that the atmospheric CO, acts like a perfect reflector for its 
characteristic radiation and that consequently the surface behaves 
like a gray radiator with the absorption coefficient (1 — €) for 
infrared radiation. With 74 = 0 and (1 — €) = 0.65 + 0.05, 7's is 
found to be 408° + 8° K. The result of this admittedly quite crude 
procedure is nevertheless of interest, as it indicates that the at- 
mosphere of Venus cannot raise the temperature of the surface by 
more than 50° at most. The actual surface temperature of the sub- 
solar point will be somewhere between 366° and 408° K, probably 
closer to the second value. In any case, it appears to be higher than 
the terrestrial boiling-point of water. 
Rupert WILDT 
PRINCETON UNIVERSITY OBSERVATORY 


November 1939 


5H. N. Russell, Ap. J., 43, 173, 1916. 
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NEW WHITE DWARFS, SUBDWARFS 
AND BINARY STARS* 


This note is a continuation of one published ten months ago.’ 
Only some of the more striking objects are given here. A complete 
list of the stars observed will follow later. 

Table 1 contains six new white dwarfs. Although no trigonometric 
parallaxes are as yet available for these stars, their white-dwarf na- 
ture is established from the spectra in four cases. For the two F stars 
the possibility of the stars being subdwarfs cannot be entirely ex- 
cluded without more data. 

Next follow seventy-nine probable subdwarfs. Most of the objects 
for which spectra were published earlier are included here, for the 
sake of uniformity, with the writer’s spectral types. The present 
list, together with the preceding one," forms a fairly complete list 
of the known subdwarfs with annual proper motion in excess of 
o” 30. Not included in this list are the many objects for which the 
subdwarf nature is based only on the spectral characteristics. Our 
list is selected on the basis of two criteria: transverse motion™ and 
reliable trigonometric parallax. 

It is possible that a few objects now classed as subdwarfs do not 
belong to that class because the proper motion found in the blink 
microscope is considerably too large, or perhaps even spurious. Such 
cases may be eliminated by a redetermination of the proper motion 
or by a closer study of the spectrum, which in that case should not 
show the subdwarf features.* In this manner the spurious proper 
motion of +54°2461, 10™2 Ao, was detected. Cases of this kind will, 
however, be quite rare. 

The assignment of some faint G stars to the subdwarfs rather than 
to the white dwarfs is tentative and is based on the presence of the 
G band in the spectrum, which seems to be absent in white dwarfs. 

The present paper contains the first M-type subdwarfs definitely 
established. All three stars appear to be about 2 mag. below the 

* Contributions from the McDonald Observatory, University of Texas, No. 18. 


t Ap. J., 89, 548, 1939. 
? For a discussion of these features, as well as of other properties and their interpre- 


tation, cf. Ann.les d’astrophysique (now in press). 
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/ 
TABLE 1 
Fw 4 | | eo | 

Star | a | 6 | Mvyis “ Sp x (Sp) | m (tr.) 
Wolf 1. | 0:08 5 | — o:14 | 15.0 0"53 | Ao | 
Wolf 1516.... | 1 12.7 | +15 40 13.6 0.69 | Con | 
Ross 548..... | 131.2] —a1 51 | 13-5 | © a | Ao | 
Ross 22 | 3 57-9 | +48 56 | 14.1 | 0.77 F 
Ross 640 2.) 26 24.9: | --36 58 | 13-2-| 0-80:| F 
Wolf 672A....../ 17 13.6 | + 204] 14.4 | 0.56 | Wh.* | 

| | + = 

+71°31.. | © 37-2 | +71 38] 10.2 | 0.35 | Fo | of003 | o%004+ 
Wolf 1504.. | o 39.1 | +29 22 | 13.0 | ©.49 | Gs | .003 | 
Wolf 56 | 1 o1.9 | +62 59 | 11.5] 1 4 .2+ |] 010 
Wolf 79 | 1 32.6 | + 9 49 | 12.0 | 0.40 | Go | .003 
Wolf 1065 | x 37:4 | --36 38 | 12.3 | 0:48 | Go: | 2003 
+72°94 | zr 38.9 | +72 58 | 10.1 | 0.30 | Fo jl 200.1 
Wolf 107 rss.7 | + 3 33 | 11-7 | 0.42 | Ge 004 
Wolf 110 157.4| + 5 14| 11.9 | 2.41 K4 } -O12 | 038 
Wolf 125 208.2 | +15 31 | 13.4 | 1.08 | Ki | .004 |— .0064 
Ross 18 209.2 | +31 56] 11.8] 0.48] G2 | .004 
Ross 557 | 2.38.4 | --e7 1s | 1327 | 0.54 | ‘Gs | .002 
+ 33°529 | 245.9 | +34 00] 9.7|1.39 | K3+ | .029 066 4 
Wolf 134 | 307.6] +18 29 | 14.0] 1.74 | K3e | 004 050 
Ross 348 | 3 10.1 | +51 56 | 13.4 | 0.41] G | .002 
Ross 373 | 3 82:0:] 4-23 56 | -I4:0 | 0.70 | “Ges 002 
Ross 374 if 3°23.1 | --23 26] 10:8 | 0.43 | Fs 003 
+66°268 3° 258: hh 66 Zeus) 029° |-5.505) 5: ) || OFS 023+ 
Ross 578 | 3 33.4 | —1I 45 | 12.1 | 3:16 | M2 025 o60 + 
Ross 581 3 36.2 | — 8 09 | 12.8] 0.35 | G8 | .004 
Ross 31 4 26.5 | +50 23] 13.6} 0.44] K2 | .004 
Wolf 1540 4:490:0 | “+ 6:54.) 1350 10.35; FS. | -.c02 
Ross 385 4 56.6 | +15 2 3.6.1 0.35.1 G2 | ;002 
Kapteyn’s star 5 07.7 | —44 50 8.8 | 8.78 | Mo 115 262+ 
Ross 65 5 16.6 | +33 06 | 12.0] 0.74 | Ki+ | .007 
+19°1185A 5 57.3 | +19 23 9.210.903 | Gr | .o13 028 + 
+19°1185B S54. 3.) =¥O. 22.1 2208510.07 | Bee | 0067 028+ 
+31°1684. 7 47.2 | +30 55 $23.1, 2:06: | ‘Ge 020 036 
Ross 618 8 04.5 | +10 28 | 13.6] 0.40 | G5 002 | 
+54°1216 | 811.7 | +54 25 9.5|0.64| F2 | .006 | 008 + 
—15°2546.. 8 36.2 | —15 50 9.6|0.60| F2 | .005 025 
— 12°2669 | 8 42.0 | —13 00 | 10.0 | 0.38 | Az2 002 = .O15 
Ross 683 8 45.1 | + 8 00 | 11.4 | 0.62 | G5 005 
— 3°2525 |  Sgace } Hg. 37 9.8] 0.76] F5 006 026 
Ross 885 9 08.6 | +20 27 | 13.4 | 0.46 | Ki 004 
+1°2341 | 9 35-6| + 1 29 | 10.6| 0.54] Fo | .003 
+44°1910 | 9 43.2 | +44 46 | 10.9 | 0.21 | As | .002 
+14°2151 | 9 43.5 | +14 14| 8.4] 0.83 | Fo | .008 | 026+ 
Wolf 336 0.55.5 | 33.25 | 12:7 | 0.38 | G8 | .oc4 
Wolf 338 9 56.2 | +35 09 | 15.6 | 0.32 | K3: | .002 
Ross 891 [0°21..7 | — © $i | 43.2] 0.34.1 Ge-.}. 002 
+29°2091 | 10 42.0 | +28 56 | 10.4 | 1.0 | F8 006 o164 
Wolf 365 | 11 05.9 | + 6 59 | 10.8 | 0.82 | G2 006 0.0224 
(L1259-59) |} rr 12.4 | +17 48 | 13.7 | 0.88 | Gs 0.002 t 





* Provisional; a spectrum in the visual region shows the star to be white; from the companion 
we find the absolute magnitude to be about 11. 


t Identification uncertain. 
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TABLE 1—Continued 









































| 
Star a ry | Mvyis | mn | Sp x (Sp) w (tr.) 

¢ UMa(M)......] 11:12.9 | +32:06 4.6 | 0773 | Fo "09 0"138+6 
Ross 111.. ..| 11 25.8 | +59 45 | 13.0 | 0.68 | K3+ re a a eee 
AC77°4245.... Il 26 +77 12 | 11.3 | 0.60 | Go COR ATs eee 
Ross 451.. | t1 34.7 | +67 53 | 22.3.1 3.20) Ba O10 O41+7 
Ross 452 woeeef EE 84.2 | +68 23 | 22.2) 0.56 | Ge Calas Dec on ees 
Ross 636 12 o1.1 | +10 46 | 12.4 | 0.38 | Go os Re el nee 
Wolf 1440 ..| 12 12.2 | -2r 37 | 2.1 | 0.74 | Gs Ce am Geer ee ean ae 
Wolf 458 12 55.2 | + 607] 14.5 |(0.6) | Go .OOI t 
+10°2519A.. 13 06.4 | +10 09 8.6 | 0.60 | G2 .018 .051+6 
+10°2519B ..| 13 06.4 | +10 09 | 12.3 | 0.60 | K5 O14 .051+6 
+34°2476. ..| 13 54.8 | +34 23 | 10.3] 0.54] Fo 003 .023+7 
+6°2932..... 14 38.4 | + 615 ]| 10.4] 0.91 | G3 .009 |— .046+9 
— 21°4009 14 54.2 | —21 36 $:5 | 0:99 F2 009 .030+9 
—15°4042 1 15 04.7 | —1I5 54 9.1 | 3.68 | G8 .020 .040+4 

—15°4041 ..| 15 04.7 | —15 59 9.4 | 3.68 | Ko 020 .040+4 
—10°4149 ..1 38 37.7 | —10 36 7.2|1.18| Fa .o16 .033+6 
+42°2667 15 59.9 | +42 32 | 10.0] 0.51 F4 Cie bos ares set 
L1130-91 16 07.3 | + 5 47 | 12.1 | 0.69 | Kz CE: AS a ane a 
Ross 530 16 15.6 | +22 53] 11.4] 0.47 | G2 a a See ere 
+ 13°3683 18 28.7 | +13 03 | 10.6] 0.31 F3 003 |— .002+10 
+ 35°3650sf 19 27.6 | +35 57 | 10.6 | 0.56] Fa2 003 007+7 
L—Ebb. 21 20 05.8 | +56 52 | 13.4 | 0.50] Ka OBA Perkaatinn eee 
Wolf 1106 21 05.6 | +59 22 | 13.0] 2.14 | Mo O17 037+8 
Wolf 923 ar 28.5 | — 7 17 | 15.2 Leese | Ra: no a7 Me Rs 2 e 
Her. 39 21 50.8 | +32 10] 11.0] 0.75 | G2 006 |— .o24+11 
+17°4708 22 06.7 | +17 36] 9.3 | 0.54] Fa2 cn ea ee eS 
Ross 272 22 06.8 | +17 34 | 11.9 | 0.37 | G2 Cie Lerevnc sta 
Wolf 1216 22 16.2 | +29 17 | 4.316.975 | Hg: CGeF hos. cee 
Ross 666 22 22.6 | +54 02 | 12.2 | 0.42 A2 ool Pe ewan A Tiel 
Wolf 1037 22 23.7 + 5 20] 14.0] 1.57 | Ka: 003 052+8 
Ross 288 22 39.7 | — 2 52 | 25.4 | 0.85 G8 CGR is .sisoeeee 

-O 4470 23 04.4 | + 0 12 Io. 2 | 3:24 G2 a ee a 
+ 38°4055 | 23 08.9 | +38 53 | 11.0 | 0.57 F8 004 Re no 
— 14°6437AB i. 2% 356 ~—14 22 | &.3 | 1.37 F8 os o19+7 
+59°2723 23 21.9 | +60 04 | 10.4 | 0.46 F4 004 .o19g+9 
—42°16457 | 23 41.2 —42 07] 7.5 | 0.88 A3 0.008 0.026+7 

| 
$ Wolf found » = 061; Luyten states that uw is smaller but does not give its value. 


average main sequence. It is of interest to note that Ross 578 is 
shown also spectroscopically to be different from ordinary M dwarfs. 

The most extreme subdwarfs thus far found are Wolf 134 and 
Wolf 1037. From the proper motion we find that they must be more 
than 3.5 mag. below the main sequence; the trigonometric parallaxes 
suggest that the deviation may be as much as about 6 mag. Addi- 
tional parallax determinations will be of interest, particularly for 
these two stars. 

There are twenty-six entries in Table 1 for which independent 
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trigonometric parallaxes are available and for which 7(Sp) < o.o1y. 
These stars would have been classified as subdwarfs even if the trigo- 
nometric parallax had not been known. It is of interest to obtain 
from these stars the most probable ratio of 7/a(Sp). Taking into 
account the greatly different weights of the twenty-six determina- 
tions, we find for the weighted mean value 2.2 + 0.2 (p.e.). This 


TABLE 2 


| | ces 
Star | a | 6 | Bu | w (tr. | My ; op. | 


) ® (Sp.) | d 
—13°249A 1:18.0 | —13:29 | 0.46 | o”oggt10] 7.5 | G8 0"042 
—13°240B. ..| | | 10.5 |}dK6 | .042 | 40 
+3°275A....| 1 56.9 iF. 3°2e) 491] 10.0 | K3+) .025 
+3°275B | 13.0:| M2 026:| 15 
L-Ebb. 7A...| 2 30.8 | +31 37 40 13.5 | M2 021 
L-Ebb. 7B 14.3 | 3 
Ross 45A....| 5 28.7 | +10 16 36 12.4} Ms 096 
Ross 45B 14.4 | M6 060 5 
— 3°2001 7 35.0 | — 3 22 20 ogo+7 7.2 | K2+ 074 
— 3°2002 8.7 |dK5 072 | 50 
Ross 825A 21 07.1 | +33 07 48 11.9\] ¢, - 
Ross 825B.. .| | 13.1f| 3+ | -— r 
Ross 773A...| 21 12.8 | +20 28 | .55 12.2 | M3 05 
Ross 773B...| 13.6 | M4 04 22 
Ross 223A...| 22 41.2 | +44 02 | 0.34 | 0.033+9 | 10.9 | Ky 0.021 
Ross 223B.. .| | 12.9 2 





means that the harmonic mean value of the transverse motion for 
these stars is 215 km/sec. This mean ratio corresponds to a devia- 
tion of 1.7 + 0.2 mag. from the average main sequence. The range 
seems to be from about 1 mag. (a somewhat arbitrary limit) to per- 
haps 5 mag. below the main sequence. 

Table 2 contains some new binaries found recently. For the two 
wide pairs proper-motion data will be required to establish the physi- 
cal relationship definitely. 

G. P. KULPER 
McDONALD AND YERKES OBSERVATORIES 
January 1940 


